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experimentally .  An  Exp&Umentat  HicAXMive.  Imaging  fiaci&ity  is  used  to  carry 
out  scattering  measurements  of  scale-models  of  targets  of  interest  in  an 
anechoic  chamber  environment.  Automated  scattered  field  measurements  with  un-. 
precedented  accuracy  are  possible  for  any  state  of  polarization  in  the  (2-26.5) 
GHz  range  with  extension  to  40  GHz  nearly  completed.  Smart  data  acquisition 
with  which  one  cam  efficiently  access  the  3-D  Fourier  space  of  a  scattering 
object  over  any  range  of  aspect  angles  are  made  possible  with  the  facility. 

Two  equivalent  image  reconstruction  algorithms:  TouJiizA.  invcuion  of  inter¬ 
polated  polar  formated  data  and  ^ittziicd  back  pfuojcctiOYi  have  been  developed 
and  refined.  The  filtered  back -projection  algorithm  is  found  to  be  more 
directly  related  to  the  physical  scattering  mechanism  of  plane  wave  illumina¬ 
tion  at  a  conducting  target.  Our  research  has  led  to  the  development  of  a 
tOAgct  dzAivcd  KC\ (zAencc  (TP R)  technique  in  which  aspect  dependent  fuxngc-pKofaiJLi 
information  is  utilized  to  form  a  iinogham  HCptiZA entaXion  of  the  scattering 
centers  on  the  target.  When  the  multiaspect  range  profiles  of  the  target  are 
aligned  properly,  a  clear  sinogram  pattern  emerges.  Point-like  scattering 
centers  on  the  target  give  rise  to  a  sinusoidal  trace  or  segment  thereof  in 
the  pattern.  A  direct  consequence  of  TDR  is  that  a  common  phase  reference  is 
synthesized  on  the  target  and  that  the  reconstructed  image  is  always  obtained 
centered  within  the  image  plane.  In  TDR,  the  phase  of  data  acquired  by  arrays 
of  independent  pairs  of  broad-band  coherent  transmitter  receivers  that  do  not 
share  a  common  phase  reference  can  be  referred  to  a  common  central  point  on 
the  target.  TDR  also  makes  the  imaging  process  more  immune  to  atmospheric 
effects.  This  is  specially  useful  when  millimeter  wave  spectral  windows  are 
used  for  higher  resolution.  Because  of  the  nature  of  scattering  mechanisms  of 
microwave  (yw)  and  millimeter  wave  (mmw)  illumination  from  man-made  objects, 
images  of  such  targets  appear  to  be  edge  enhanced.  The  image  thus  contains 
only  information  about  characteristic  detail  on  the  target,  such  as  edges, 
outlines,  protrusions  etc.,  which  aids  recognitions  by  the  eye-brain  system 
and  has  attributes  of  a  pfuxnat  sketch  that  makes  it  also  very  useful  for  machim 
recognition.  Techniques  for  accessing  the  Fourier  space  of  a  moving  6cattZAeA 
with  the  same  accuracy  realized  with  stationary  objects  in  our  imaging 
facility  have  also  been  studied.  A  doppler  compensation  scheme  employing  high 
speed  frequency  measurement  and  synthesis  is  found  to  be  the  most  promising. 

Resolution  in  yw/mmw  diversity  imaging  is  proportional,  among  other  things, 
to  the  width  of  the  spectral  window  over  which  data  is  acquired.  Very  wide 
spectral  windows  can  in  principle  yield  images  of  remote  targets  e.g.,  aero¬ 
space  targets  with  near  optical  resolution.  The  use  of  such  extended  spectral 
range  is  not  possible  in  practice  because  of  equipment  limitation,  frequency 
band  allocation,  and  atmospheric  effects.  It  is  possible  however  to  conduct 
measurements  over  separate  segments  of  an  extended  spectral  range .  For  this 
reason  a  new  iterative  algorithm  for  extrapolation  of  data  available  in 
multiple  segmented  regions  covering  an  extended  spectral  range  has  been 
developed.  For  objects  with  discrete  scattering  centers,  the  algorithm  is 
shown  to  furnish  images  with  quality  equal  to  the  image  obtained  had  the  data 
been  collected  continuously  over  the  entire  extended  spectral  range.  This 
result  has  obvious  practical  implications  in  that  hardware  complexity  is 
traded  by  cheaper  and  more  readily  available  computational  expenditure.  At 
the  same  time  limitations  imposed  by  frequency  allocation  and  atmospheric 
effects  are  cleverly  bypassed. 

Atmospheric  effects  limit  the  spectral  range  useful  for  microwave  diversity 
imaging  to  the  lower  mmw  region.  As  a  result  the  aperture  (physical  or 
synthetic)  needed  to  realize  high  resolution  on  remote  aerospace  targets  turns 
out  to  be  very  large  as  it  must  collect  data  for  as  wide  a  range  of  aspect 
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angles  of  the  target  as  possible.  Zt  is  worth  recalling  that  different  aspect 
"looks"  at  the  target  is  what  is  needed  for  image  formation.  Although  feasible, 
highly  thinned,  extended  broad-band  aperture  of  the  type  needed  can  be  very  costly. 

A  single  coherent  broad-band  transmitter/receiver  can  be  used  then  in  an  ISAR 
(inverse  SAR)  mode  to  collect  the  data  as  the  target  is  moving.  This  however 
is  a  non-real-time  approach  and  can  be  time-consuming  and  even  unacceptable  in 
some  situations  when  a  large  aperture  is  to  be  synthesized  in  order  to  achieve 
high  resolution.  One  is  faced  then  with  the  classical  problem  of  how  to  recognize 
the  target  in  microwave  diversity  imaging  from  a  few  looks  i.e.,  the  problem  of 
■  super-resolution  or  target  recognition  from  partial  information. 

For  this  reason  a  second  approach  to  ym/mow  target  identification  based  on 
machine  recognition  where  no  explicit  image  formation  takes  place  was  pursued. 

, This  approach  developed  more  recently  in  our  research  as  a  result  of  research  in 
neural  net  models  and  their  analogs  as  a  new  robust,  high-speed,  and  fault 
i  tolerant  approach  to  signal  processing.  A  research  effort  in  radar  target  recogni- 
\  tion  from  partial  information  based  on  neural  net  models  and  their  opto-electronic 
implementation  has  been  initiated.  Sinogram  representations  of  targets  of  interest 
are  considered  as  learning  sets.  The  basic  assumption  here  is  that  sinogram 
i  representations  and  other  kinds  of  signatures,  e.g.,  polarization  maps,  of  scale 
models  of  targets  of  interest  can  be  generated  cost-effectively  in  our  an echoic 
{chandler  facility  to  serve  as  learning  sets  for  a  neural  processor  that  can  be 
i subsequently  used  in  the  recognition  of  the  actual  radar  targets  from  few  echoes 
collected  by  broad-band  tracking,  radar  systems.  Careful  attention  is  being  given  to 
scaling  questions  and  to  the  principle  of  electromagnetic  similitude  in  order  to 
make  the  "laboratory"  generated  representations  be  as  realistic  as  possible  and 
therefore  useful  in  the  recognition  of  echoes  obtained  by  any  broad-band  radar 
installation.  Preliminary  results  show  that  a  neural  net  processor  is  capable  of 
distinguishing  between  models  of  three  aero-space  objects  from  a  few  looks  on  each 
corresponding  to  10  to  20  percent  of  their  corresponding  sinogram  representations. 
‘The  neural  net  processor  in  this  work  performed  the  function  of  information  storage, 
processing,  and  labeling  of  the  target  recognized  simultaneously.  Current  work  is 
considering  the  question  of  self-organization  and  learning  in  neural  nets  as  a 
swans  for  the  net  to  generate  its  own  representations  of  the  items  it  is  supposed 
to  become  familiar  with  rather  than  being  deliberately  programmed  to  do  so  as  was 
the  case  in  our  preceeding  work.  Stochastic  learning  algorithms  and  fast  opto¬ 
electronic  analog  architectures  and  implementations  based  on  simulated  annealing 
in  the  context  of  a  Boltzmann  machine  formalism  are  being  studied.  Work  in  this 
area  is  being  vigorously  pursued  as  it  might  have  far  reaching  implications  in 
smart  sensing  and  recognition  and  artificial  intelligence  systems.  Simulated 
annealing  as  an  optimization  process  has  also  been  found  in  our  work  to  be  very 
useful  in  phased -array  antenna  pattern  synthesis. 
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RESEARCH  IN  IMAGE  UNDERSTANDING  AS  APPLIED  TO  3-D  TOMOGRAPHIC 
IMAGING  WITH  NEAR  OPTICAL  RESOLUTION 


1 .  Introduction 

^  Our  research  in  high  resolution  microwave  imaging  shows  that  massive 
amounts  of  information  about  a  scattering  object  can  be  gathered  by 
measurements  that  combine  angular,  spectral  and  polarization  diversity  and 
trade-offs  between  spectral  and  angular  degrees  of  freedom  can  be  used  to 
enhance  cost  effectiveness  of  microwave  diversity  imaging  systems.  The  data 
collected  can  be  used  to  access  the  3-D  Fourier  space  or  p-space  of  the 
scatterer  which  can  be  Fourier  inverted  to  yield  3-D  image  information 
either  projectively  or  tomographically  (in  slices  or  cross-sectional 
outlines  when  the  object  is  metallic).  This  indicates  that  three- 
dimensional  microwave  imaging  of  distant  aerospace  targets  at  resolutions 
comparable  or  exceeding  the  capabilities  of  optical  imagers  (whose  operation 
is  severely  degraded  by  atmospheric  effects)  is  feasible.  The  aim  of 
research  described  in  this  report  was  to  attain  new  levels  of  image 
understanding  that  are  useful  in  the  study  and  demonstration  of  3-D 
microwave  remote  imaging  with  near  optical  resolution  or  better.  Real-time 
operation  and  cost-effectiveness  were  of  primary  concern.  '  Also  the 
presentation  of  3-D  image  detail  to  the  observer  in  such  a  fashion  as  to  aid 
recognition  by  the  eye-brain  system  was  an  objective,  although  the  study  of 
automated  identification  and  classification  technique  was  not  excluded  and 
encouraging  initial  results  in  this  regard  have  recently  been  obtained  and 
are  included  here. 
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2.  Summa 

Achievement  of  the  above  goals  required  a  new  level  of  insight  and 
understanding  of  the  dynamics  of  electro-magnetic  scattering  and  inverse 
scattering  in  general  and  those  of  microwave/millimeter  wave  diversity 
imaging  in  particular.  This  understanding  has  lead  to  the  identification  of 
imaging  algorithms  capable  of  judiciously  combining  smart  data  acquisition, 
signal  processing  including  image  enhancement  and  restoration,  and  efficient 
3-D  image  presentation  and  display.  New  levels  of  insight  and  understanding 
developed  by  the  program  have  also  contributed  to  the  development  of 
generalized  imaging  theory  applicable  to  all  modes  of  coherent  or 
incoherent,  synthetic  or  nonsynthetic  multidimensional  (2-D  and  3-D)  imaging 
systems  and  to  the  development  of  new  indirect  imaging  modalities  where  the 
object  transform  is  the  data  collected  and  from  which  an  image  is  formed. 
The  availability  of  the  object  transform  (e.g.,  Fourier  transform)  permits 
the  application  of  sophisticated  analysis  and  filtering  methods  that  play  a 
major  role  in  facilitating  the  study  of  imaging  processes.  Progress  in  all 
these  aspects  has  been  reported  in  earlier  annual  reports,  [1]-[4]  and  in 
journal  articles  and  proceedings  of  technical  meetings  (see  list  of 
publications  in  section  3  of  this  report). 

Real-time,  cost-effective  imaging  of  remote  scattering  objects  with 
near  optical  resolution  using  microwaves  (or  millimeter  waves)  provides  a 
valuable  tool  in  several  applications.  These  include  damage  assessment  of 
satellites  and  other  aero-space  objects,  identification  and  classification 
of  ground  vehicles  and  aero-space  objects  and  debris  with  low  proability  of 
error,  radar  cross-section  reduction  (diaphanization)  to  reduce  probability 
of  detection,  non-destructive  evaluation  (NDE)  of  critical  dielectric  and 
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composite  material  components  that  do  not  lend  themselves  to  other  NDE 
techniques,  and  3-D  analysis  and  study  of  storm  centers. 

Specifically,  the  approach  followed  in  this  research  is  based  on 
intelligent  combination  of  thorough  understanding  of  scattering  mechanisms, 
image  enhancement  and  restoration,  with  smart  data  acquisition  over  extended 
spectral  windows  that  start  in  the  lower  microwave  range  and  extend  into  the 
millimeter  wave  region.  The  argument  for  using  extended  spectral  windows  is 
found  in  the  results  of  our  X  and  polarization  diversity  imaging  work  in  the 
(6-18)  GHz  range  described  in  the  reports  and  publications  cited  earlier. 
These  results  made  it  abundantly  clear  that  the  utility  of  our  Experimental 
Microwave  Measurement  and  Imaging  Facility  described  in  earlier  reports  [1 ]— 
[2]  can  be  drastically  enhanced  by  extending  its  spectral  range  of  operation 
down  to  2  GHz  and  upto  40  GHz.  In  particular  the  reduction  in  the  side-lobe 
level  of  the  spatially  variant  point  spread  function  (PSF)  provided  by 
inclusion  of  lower  spectral  frequencies  and  the  ensuing  reduction  of 
speckle,  indicated  by  an  idealized  analysis  given  in  [1]  are  two  obvious 
motives.  Furthermore,  extending  operation  down  to  2  GHz  provides  probing 
wavelength  ranges  up  to  15  cm  which  exceed  the  characteristic  sizes  of  a 
number  of  sub-structures  on  the  test  objects  utilized  in  our  studies.  For 
example  the  engines,  fuel  tanks,  and  tail  section  of  a  100:1  scale  model  of 
a  B-52  we  frequently  use  as  test  target  are  all  of  the  order  of  or  less  than 
15  cm.  This  places  scattering  from  these  components  in  their  Rayleigh  or 
resonance  regime  and  enhances  image  quality  by  causing  a  "filling-up”  of  the 
image.  This  leads  for  example  to  a  desirable  alteration  of  the  image  of  the 
engines,  appearing  usually  as  two  dots  located  at  the  ends  of  their  nearly 
cylindrical  shapes,  to  a  more  continuous  image  resembling  more  closely  their 
optical  image.  It  is  possible  of  course  to  investigate  Rayleigh  and 


resonance  scattering  with  a  ( 6— 1 8 )  GHz  system  by  reducing  the  size  of  the 
test  models  employed.  This  however  will  result  in  an  unacceptable  reduction 
of  the  radar  return  necessitating  the  use  of  higher  levels  of  irradiance  at 
the  target.  The  argument  for  extending  the  operational  range  up  to  40  GHz 
stems  from  the  desire  to  collect  information  about  finer  detail  on  the 
scaled  target  representing  for  example  damaged  structure  on  an  actual 
target . 

Figure  1  shows  examples  of  two  microwave  diversity  projection  images  of 
the  scale  model  of  a  B-52  used  as  a  test  target  in  our  experimental 
Microwave  Imaging  Facility.  The  image  in  Fig.  1(a)  is  obtained  using  a 
spectral  window  of  (6-18)  GHz  while  that  in  Fig.  1(b)  is  obtained  with  a 
spectral  window  of  (2-1 8)  GHz.  The  angular  window  of  data  acquisition  in 

both  cases  was  about  90°  from  head-on  to  broad-side  in  128  looks.  The 
spectral  window  were  covered  by  incrementing  the  frequency  in  steps  of  “86 

MHz.  The  elevation  angle  in  both  cases  was  30°.  It  is  easily  verified  that 
the  (2-18)  GHz  image  exhibits  improved  delineation  of  detail  such  as  the 
engines  and  fuel  tanks  as  compared  to  the  (6-18)  GHz  image.  Also  the 
fuselage  d  sail  (vertical  tail  sections)  are  better  delineated.  The 
outlines  of  the  wings  and  the  tail  are  however  de-emphasized  somewhat  in  the 
(2-1 8)  GHz  image  because  of  the  effect  of  low  frequency  scattering  data 
included  in  this  image.  Extension  of  the  upper  limit  spectral  window  for 
example  to  40  GHz  while  maintaining  the  2  GHz  lower  limit  is  expected  to 
restore  the  sharpness  of  the  wings  and  tail  section  and  to  result  in  further 
image  improvement. 

The  images  shown  are  symmetry  enhanced  and  polarization  enhanced  as 
described  in  earlier  reports  and  publications  and  have  been  obtained  by 
means  of  an  elaborate  target  derived  reference  technique  (TDR)  developed  and 


refined  in  our  work  which  removes  the  range-phase  caused  by  the  unknown 
range  to  a  specific  reference  phase  center  on  the  target  in  order  to  be  able 
to  form  an  intelligible  image. 

Resolution  in  yw/mmw  diversity  imaging  depends  among  other  things  on 
the  width  of  the  available  spectral  window.  Extended  spectral  windows 
covering  the  yw  and  portions  of  the  mmw  range  can  lead  to  resolution  on 
remote  targets  (e.g.,  aerospace  targets)  exceeding  those  of  optical  systems 
which  are  more  susceptible  to  atmospheric  effects  than  yw/mmw  imaging 
systems.  Hardware  availability  and  frequency  band  allocation  may  limit 
however  data  acquisition  in  practice  to  nonoverlapping  bands  that  can  cover 
a  wide  spectral  range.  The  question  then  is  how  to  employ  information 
collected  in  such  multiple  restricted  frequency  bands  to  retrieve  by 
extrapolation  into  the  missing  bands  an  image  of  the  target  equivalent  or 
approaching  the  quality  of  that  obtained  had  the  entire  un-interrupted 
spectral  range  been  available  . 

A  new  iterative  method  for  extrapolation  of  incomplete  segmented  data 
available  in  multiple  separated  bands  has  therefore  been  proposed  and 
tested.  The  method  uses  the  Burg  algorithm  to  find  the  linear  prediction 
parameters  and  an  iterative  procedure  to  improve  the  estimation  of  the 
linear  prediction  parameters  and  extrapolation  of  the  data.  This  method  is 
especially  effective  when  the  spectra  (Fourier  transform  of  the  observed 
data)  are  in  discrete  forms.  In  the  context  of  radar  imaging  represented 
here,  this  means  the  object  consist  of  discrete  or  distinctly  spaced 
scattering  centers.  The  advantage  of  this  algorithm  were  demonstrated  using 
both  numerically  generated  and  realistic  experimental  data  pertaining  to 
high  resolution  radar  imaging.  Detail  of  this  work  is  given  in  Appendix  V. 


The  image  formation  principle  for  microwave  diversity  imaging  systems 
can  be  understood  by  investigating  the  scattering  mechanisms  of  the  object 
and  the  image  reconstruction  algorithms.  After  looking  into  these  two 
factors,  several  numerical  and  experimental  results  can  be  satisfactorily 
explained.  The  image  of  a  metallic  object  obtained  using  prescribed 
spectral  and  angular  windows  can  be  predicted  accordingly.  This  successful 
interpretation  is  fundamental  to  research  in  radar  cross  section  management 
studies  and  target  imaging  and  recognition.  Recent  research  in  this  regard 
and  the  results  obtained  are  summarized  in  some  detail  in  Appendix  I. 

The  ultimate  aim  of  all  imaging  systems  is  recognition  and 
identification.  In  the  microwave  diversity  imaging  work  described  above, 
the  recognition  and  identification  is  performed  by  an  observer  i.e. 
identification  is  done  by  the  eye-brain  system. 

Another  approach,  to  target  identification  which  emerged  more  recently 
in  our  research,  is  less  traditional  than  mere  formation  of  a  high 
resolution  image  suitable  for  human  analysis.  It  involves  automated  machine 
recognition.  Here  we  are  concerned  with  Issues  of  correct  recognition  and 
labeling  by  a  machine  from  partial  or  sketchy  information  and  with  systems 
that  can  do  this  in  a  robust  and  fault  tolerant  manner.  In  this  new 
approach,  the  role  of  the  eye-brain  system  in  recognizing  an  image  is 
mimicked  by  the  machine  which  makes  the  recognition  process  "brain-like"  in 
nature.  The  approach  is  based  on  neural  net  models  and  their  opto¬ 
electronic  analogs  which  exploits  the  fit  between  what  optics  can  offer 
(parallelism  and  massive  interconnectivity)  and  what  neural  net  models  can 
offer  (new  approach  to  signal  processing  that  is  nonlinear  and  collective 
and  therefore  f ast ,. robust ,  and  fault  tolerant).  The  two  approaches  are 
coupled  in  that  automated  machine  recognition  with  artificial  neural 


networks  relies  on  the  generation  of  target  or  object  representations  or 
feature  spaces  that  can  lead  to  "distortion  tolerant"  recognition  i.e., 
recognition  irrespective  of  target  range,  orientation,  or  location  within 
the  field  of  view  i.e.,  what  is  traditionally  referred  to  as  scale, 
rotation,  and  shift  invariant  recognition.  The  generation  of  such 
representations  usually  involves  the  same  gear  employed  in  traditional 
iiw/mmw  diversity  imaging  where  spectral,  angular,  and  polarization  degrees 
of  freedom  are  combined  to  realize  images  of  the  scattering  targets  with 
near  optical  resolution.  In  fact  the  representations  or  feature  spaces 
often  contain  exactly  the  same  information  contained  in  a  yw/mmw  image  of 
the  target  except  it  is  arranged  in  a  different  format  that  is  more  amenable 
for  use  in  automated  recognition  schemes.  This  approach  has  lead  to  the 
emergence  of  a  new  concept  of  achieving  the  ultimate  resolution  thrugh 
correct  recognition  of  the  target  by  a  machine.  In  other  words  correct 
recognition  of  the  target  by  a  machine  is  synonymous  with  producing  a  very 
high  resolution  image  with  detail  sufficient  for  reliable  recognition  by  the 
eye-brain  system  and  for  discrimination  between  closely  similar  targets. 

Among  its  many  fascinating  capabilities  such  as  robustness  and  fault 
tolerance,  the  brain  is  also  able  to  recognize  objects  from  partial 
information.  We  can  recognize  a  partially  obscured  or  shadowed  face  of  an 
acquaintance  or  a  multilated  photograph  of  someone  we  know  with  little 
difficulty.  The  brain  has  a  knack  for  supplementing  missing  information. 
Capitalizing  on  this  feature  and  on  our  knowledge  of  neural  models  and  their 
collective  computational  properties,  (see  15  and  16  of  list  of  publications) 
a  study  of  "neural  processing"  for  recognizing  objects  from  partial 
information  has  been  initiated.  An  example  of  the  capabilities  of  this  new 


and  intriguing  approach  in  automated  recognition  and  labeling  of  radar 
targets  from  partial  information  is  given  in  Fig.  2.  Details  of  this 
preliminary  work  is  given  in  Appendix  IV. 

Object  recognition  employing  hetero-associative  memory  based  on  models 
of  neural  networks  may  obviate  the  need  for  expensive  or  non  real-time  image 
formation  in  the  field.  Images  of  scale  models  of  objects  of  interest  and 
their  sinogram  representations  are  produced  indoors  in  our  microwave  imaging 
facility  at  a  fraction  of  the  cost  of  an  actual  radar  system  operating  to 
form  images  of  actual  targets  intended  for  recognition  by  the  eye-brain 
system.  (See  Appendix  IV).  The  sinogram  representations  which  can  be 
viewed  also  as  signatures  or  feature-spaces  of  targets  of  interest  are 
taught  to  a  hetero-associative  memory  by  associating  them  with  identifying 
letter  labels.  When  such  a  hetero-associative  memory  is  presented  with 
partial  sinogram  information,  that  can  be  as  low  as  1 0%  of  the  full 
sinogram,  the  memory  will  supplement  the  missing  information  and  produce  the 
correct  label.  This  ability  of  supplementing  missing  information  is 
synonymous  with  "super  resolution".  It  is  akin  to  our  ability  for  example 
to  recognize  a  photograph  of  someone  we  know  even  when  part  of  it  is  missing 
or  to  understand  a  sentence  in  a  handwritten  letter  even  when  some  of  the 
words  are  not  legible.  Partial  sinogram  information  can  be  produced  much 
more  rapidly,  easily,  and  economically  in  the  field  by  actual  radar  systems 
than  producing  full  sinogram  or  image.  Thus  a  few  looks  at  the  target  from 
slightly  different  aspects  may  be  sufficient  for  recognition.  The  real 
challenge  is  to  recognize  the  target  from  one  look  i.e.,  from  one  broadband 
radar  echo  and  also  for  non-aerospace  targets  in  the  presence  of  clutter. 
The  achievement  of  this  aim  employing  artificial  neural  processors  is  a 
major  objective  of  our  current  and  future  research.  It  is  important  to  note 
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that,  unlike  a  conventional  memory  in  a  computer,  the  hetero-asaociative 
memory  here  performs  the  functions  of  storage,  processing,  and  recognition 
simultaneously.  As  such  it  is  more  correct  to  regard  it  as  a  "neural 
processor".  A  primary  objective  of  our  future  research  in  this  area  is  the 
study  of  generating  realistic  target  representation  employing  scale  models 

in  an  anechoic  chamber  environment.  Careful  consideration  of  scaling  issues 

w 

and  of  the  principle  of  electromagnetic  similitude  are  called  for  in  order 
to  be  able  to  generate  from  scale  models  feature  spaces  (e.g.,  sinogram  and 
polarization  representations)  that  resemble  those  obtained  had  the  actual 
object  been  interrogated  in  the  field  by  a  realistic  broad-band  radar 
system. 

Our  research  during  this  period  examined  also  several  approaches  to  the 
problem  of  accessing  the  3-D  Fourier  space  of  a  scattering  body  undergoing 
irrotational  motion.  The  criterion  for  selecting  the  best  approach  were 
speed,  accuracy  and  economy.  The  method  chosen,  shown  in  Fig.  3.  exploits 
recent  advances  in  the  state-of-the-art  of  Microwave  instrumentation  to 
enable  the  study  of  innovative  and  efficient  approaches  to  data  acquisition 
that  were  generally  thought  not  to  be  feasible  a  3hort  time  ago.  The  scheme 
utilizes  two  capabilities  of  current  microwave  instrumentation.  High  speed 
computer  controlled  frequency  synthesis  and  high-speed  frequency 
measurement.  These  capabilities  are  combined  in  the  arrangement  shown  to 
realize  rapid  determination  of  the  doppler  shift  FD  in  the  signal  scattered 

from  the  moving  scatterer  (target)  which  can  then  be  compensated  for  in  the 


phase  and  amplitude  measurement  of  the  scattered  field  carried  out  by  the 
coherent  receiver  (hp  841  OB  Network  Analyzer)  whose  intermediate 

*  J.D.  Stratton,  Electromagnetic  Theory,  McGraw  Hill,  New  York  (1941).  pp. 
488-490. 


Fig.  Doppler  compensated  arrangement  for  measuring  the  stepped 
frequency  response  of  a  moving  scatterer. 


local  oscillator  Is  effectively  made  to  equal  FRp  +  FD  using  an  externally 

synthesized  signal  at  frequency  FIF  *  FD  for  the  reference  IF  channel  of  the 

network  analyzer.  This  eliminates  the  effect  of  Doppler  shift  on  the  phase 
measurement  and  provides  data  very  much  similar  to  that  generated  had  the 
target  been  stationary  and  it  merely  changes  its  aspect  relative  to  the  line 
of  sight  of  the  Interrogating  radar  system.  Note  the  low  frequency 
synthesizer  source  of  the  fif  *  fd  signal,  the  hp  3325A  synthesizer,  is 

phase-locked  to  the  RF  source  of  illumination  of  the  target.  This 
synthesized  source  has  been  carefully  chosen  because  its  high  frequency 
resolution  which  enables  setting  the  doppler  shift  of  the  IF  reference 
precisely  to  the  measured  value  of  F^.  This  frequency  resolution  is 

necessary  because  any  discrepancy  between  the  Intermediate  frequency  (IF) 
reference  signal  and  the  down  converted  IF  version  of  the  target  echo  signal 
in  the  network  analyzer  results  in  an  undesirable  time  varying  phase  error. 
In  this  arrangement  the  complex  frequency  response  of  the  moving  scatterer 
is  measured  rapidly  in  discrete  frequency  steps  following  the  determination 
of  the  doppler  shift  and  the  setting  of  the  low  frequency  intermediate 
frequency  reference  synthesizer  to  (Fip  ♦  Fp).  If  we  assume  that  100 

frequency  stepped  measurements  per  second,  as  limited  by  the  bandwidth  of 
the  phase/amplitude  measuring  segment  of  the  network  analyzer,  can  be 
performed,  then  the  Doppler  shift  during  this  interval  is  expected  to  change 
negligibly  for  most  practical  targets  including  aerospace  targets. 
Therefore  the  Doppler  measurement  and  doppler  compensation  operation  need 
not  be  performed  frequently  by  the  system  but  only  occasionally  or  at  most 
once  at  the  beginning  of  each  digital  frequency  sweep. 


All  the  instrumentation  to  verify  this  scheme  is  now  available  in  our 
laboratory  with  the  exception  of  the  hp  3325A  low  frequency  synthesizer.  In 
the  next  phase  of  this  research  we  plan  to  either  purchase  or  lease  this 
instrument  in  order  to  carry  out  this  experiment. 

As  described  earlier,  the  smart  sensing  capabilities  of  our  yw/mmw 
diversity  imaging  approach  were  exploited  to  generate  target  representations 
(sinogram  representations  or  feature  spaces)  which  were  used  in  turn  to 
compute  the  synaptic  connectivity  matrix  of  the  neural  net  processor.  This 
operation  is  equivalent  to  deliberate  programming  of  the  neural  net 
processor  to  perform  certain  pattern  recognition  task.  One  of  the  most 
intriguing  and  potentially  useful  properties  of  neural  nets  is  their  ability 
to  self-organize  and  learn  where  the  net  forms  its  own  internal 
representations  of  the  associations  it  is  presented  with  by  internally 
altering  the  weights  of  interconnections  between  its  neurons  in  accordance 
to  sane  learning  algorithm.  Self-organization  and  learning  in  neural  nets 
by  simulated  annealing  in  the  context  of  a  Boltzmann  machine  formalism  is 
one  of  their  attributes  that  sets  them  apart  from  other  approaches  to  signal 
processing.  Self  organization  and  learning  in  neural  nets  can  play  an 
important  role  in  the  development  of  smart  sensing,  artificial  intelligence, 
and  high-speed  knowledge  processing  systems. 

Simulated  annealing  is  known  to  be  an  effective  tool  in  the  solution  of 
combinatorial  optimization  problems  where  an  "energy”  or  "cost"  function  is 
to  be  minimized  as  a  means  of  finding  a  solution.  It  has  been  applied  in 
the  solution  of  a  wide  range  of  combinatorial  optimization  problems.  We 
have  found  recently  that  it  can  also  be  an  effective  tool  in  phased  array 
antenna  design  where  simulated  annealing  is  applied  to  determine  the  binary 
(♦1,-1)  weights  of  the  phased  array  elements  that  yield  an  optimal  radiation 


pattern  as  far  as  beamwidth  and  side-lobe  level  are  concerned.  Detail  of 
this  work  is  given  in  Appendix  V. 
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IMAGE  UNDERSTANDING  AND  INTERPRETATION  IN 
MICROWAVE  DIVERSITY  IMAGING 

HsmkJyh  Li,  Nabit  H.  Farhat,  Yuhsyen  Shen,  Charles  Werner 

The  Moore  School  of  Electrical  Engineering 
University  of  Pennsylvania 
Philadelphia,  PA.  19104 


ABSTRACT 

The  image  formation  principle  for  microwave  diversity  imaging  systems 
can  be  understood  from  investigating  the  scattering  mechanisms  of  the  object 
and  the  image  reconstruction  algorithms.  After  looking  into  these  two  factors, 
several  numerical  and  experimental  results  can  be  satisfactorily  explained.  The 
image  of  a  metallic  object  obtained  by  prescribed  spectral  and  angular  win¬ 
dows  can  be  predicted  accordingly.  This  successful  interpretation  is  fundamen¬ 
tal  to  research  in  radar  cross  section  management  studies  and  target  recogni- 
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IMAGE  UNDERSTANDING  AND  INTERPRETATION  IN 
MICROWAVE  DIVERSITY  IMAGING 

Hsueh-Jyh  Li,  Nabil  H.  Far  hat,  Yuhsyen  Shen,  Charles  Werner 

The  Moore  School  of  Electrical  Engineering 
University  of  Pennsylvania 
Philadelphia,  PA.  19104 


L  INTRODUCTION 

The  image  properties  obtained  by  microwaves  are  different  from  those  obtained  by  opti¬ 
cal  waves.  The  relative  long  wavelength  of  microwaves  is  the  main  reason.  Optical 
wavelengths  are  usually  smaller  than  the  surface  roughness  of  most  objects.  The  surface  can 
be  seen  through  the  waves  scattered  from  the  rough  surface.  The  same  surface  might  be 
smooth  compared  to  microwave  wavelengths,  and  therefore  its  microwave  image  can  be 
entirely  different  from  its  optical  counterpart 

Determination  of  the  shape  of  a  perfectly  conducting  body  from  a  finite  number  of  meas¬ 
urement  of  the  scattered  fields  has  been  considered  as  an  inverse  scattering  problem  [1,2,3], 
where  the  input  (the  illuminating  plane  wave)  is  known  and  the  outputs  (the  scattered  far 
fields)  are  measured,  and  the  object  function  to  be  determined  is  the  geometrical  shape  of  the 
body.  It  has  been  proven  that  under  the  physical  optics  (PO)  approximation,  the  object  func¬ 
tion  could  be  determined  if  the  back- scattered  fields  could  be  measured  for  all  frequencies  and 
all  directions[2,3].  In  reality,  the  spectral  and  angular  windows  for  the  data  are  limited.  Furth¬ 
ermore,  the  PO  approximation  will  fail  when  the  wavelength  is  not  small  compared  to  the 
characteristic  dimension  of  the  object.  Therefore,  the  microwave  image  reconstructed  over 
finite  spectral  and  angular  windows  are  not  exactly  the  geometrical  shape  of  the  object 

Microwave  images  have  been  obtained  by  several  different  schemes,  for  example,  side¬ 
looking  radar,  synthetic  aperture  radar  (SAR),  inverse  synthetic  aperture  radar  (ISAR)[4], 
radio  camera  [5],  microwave  holography  [6],  microwave  diversity  imaging  [7,8],  etc.  Different 
schemes  may  produce  different  images  and  the  images  were  interpreted  in  different  ways. 
Radar  workers  modeled  the  object  with  a  reflectivity  distribution  function  [9],  while  elec¬ 
tromagnetics  workers  modeled  the  object  from  more  rigorous  expressions  [1,2,3].  The  nature 
of  the  reconstructed  images  depends  on  the  data  acquisition  systems  and  the  reconstruction 


algorithm.  For  example,  the  SAR  image  is  derived  using  a  range  measuring  technique  in  one 
dimension  and  the  cross  range  coordinate  is  then  introduced  by  processing  many  range  returns 
along  the  synthesized  aperture  [9].  The  effect  of  finite  range  resolution  produces  a  "breaking 
upH  of  the  image  in  the  range  direction  that  accentuates  edges  of  the  image  and  eliminates 
portions  of  the  image  which  have  constant  intensity  over  long  distances.  The  radio  camera 
method  treats  the  image  problem  from  the  antenna  design  point  of  view.  Beamforming  and 
scanning  are  the  two  main  processes  of  the  image  reconstruction  algorithm,  and  the  element 
separation  is  in  terms  of  wavelength  [5].  Microwave  diversity  imaging  systems  make  full  use 
of  angular,  wavelength,  and  polarization  diversity  and  to  obtain  images  of  metallic  objects 
with  nearly  optical  resolution  [8].  2-D  Fourier  transform  method  or  filtered  back-projection 
method  are  the  image  reconstruction  algorithms  of  that  imaging  scheme.  The  image  properties 
are  close  to  those  of  SAR  images. 

The  principle  of  microwave  diversity  imaging  is  based  on  the  PO  approximation  and 
Bom  approximation.  An  object  scattering  function  was  first  defined,  and  then  it  was  shown 
that  the  object  scattering  function  and  the  corrected  scattered  fields  are  3-d  Fourier  transform 
pair.  It  was  thought  that  the  reconstructed  image  represents  the  object  scattering  function. 
However,  if  the  position  of  the  transmitter  is  changed,  or  if  the  object  has  been  rotated,  then 
the  object  scattering  function  cannot  be  defined  clearly.  This  situation  also  occurs  in  SAR  or 
ISAR  images.  What  the  image  represents  and  how  the  image  will  look  like  over  a  given  finite 
spectral  and  angular  windows  were  not  fully  explained  and  discussed. 

The  principle  of  image  formation  for  microwave  diversity  imaging  system  can  be  under¬ 
stood  from  investigating  the  scattering  mechanism  of  the  object  and  the  procedures  of  the 
reconstruction  algorithm.  After  looking  into  these  two  factors,  several  numerical  and  experi¬ 
mental  results  can  be  satisfactorily  explained.  Also,  the  image  of  a  metallic  object  obtained  by 
prescribed  spectral  and  angular  windows  can  be  predicted  accordingly. 

In  this  paper,  the  microwave  diversity  imaging  system  will  be  briefly  described;  the 
scattering  mechanism  of  the  object  and  the  procedures  of  the  reconstruction  algorithm  will  be 
studied;  what  the  images  represent  and  what  the  image  will  look  like  over  prescribed  angular 
and  spectral  windows  will  be  interpreted  and  predicted;  some  numerical  and  experimental 
examples  of  reconstructed  images  will  be  given  and  discussed;  an  alternative  interpretation  of 
the  edge  enhancement  property  of  microwave  diversity  image  will  be  given  and  discussed. 


IL  DESCRIPTION  OF  MICROWAVE  DIVERSITY  IMAGING  SYSTEMS 


The  simplified  microwave  diversity  imaging  arrangement  is  shown  in  Fig.l.  A  metallic 
object  seated  on  a  rotating  pedestal  is  illuminated  by  a  plane  wave  with  unit  amplitude.  The 
distance  between  the  rotation  center  O  and  the  transmitter  and  receiver  are  R,  and  Rr  respec¬ 
tively,  the  unit  vectors  in  the  direction  of  transmitter  and  receiver  are  lt  and  lr  respectively, 
and  both  are  perpendicular  to  the  rotation  axis.  Using  the  PO  approximation  and  far  field  con¬ 
dition,  the  vector  potential  at  the  receiver  can  be  expressed  as  [8] 
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where  k  is  the  wave  number,  Sa/  the  illuminated  region,  h  ( 7* ' )  the  unit  normal  vector  at 
7* ',  Iff  the  unit  vector  of  the  incident  magnetic  field.  The  scattered  field  is  related  to  the  vec¬ 
tor  potential  by 

?s(kJtJr)  =  jMCr(k,itjr)  (2) 

where  XT  is  the  transverse  component  of  A  along  the  direction  lr  .  Express  k(lr  - 1,  )  as 
the  3-d  position  vector  P  in  Fourier  space  and  n(F*'  )x  lH  as  the  object  scattering  function 
T*  ( 7* ' ),  the  vector  potential  becomes 

A*  (?)  =  tfp'  j  **(**')«  (3) 
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If  the  term  in  front  the  integral  is  removed,  the  integral  in  Eq.(3)  shows  that  measurement  of 
the  vector  potential  carried  out  for  a  sufficient  range  of  values  of  P  that  are  realized  by  vary¬ 
ing  lr  (angular  diversity)  and  by  varying  k  (frequency  or  wavelength  diversity)  permits 
accessing  a  finite  region  of  the  3-d  Fourier  space  data  [8]. 

It  is  noted  that  the  object  scattering  function  is  uniquely  defined  only  if  both  the  object 
and  1^  are  fixed  during  the  measurement  interval.  This  condition  holds  only  if  the  object  and 
the  transmitting  antenna  are  kept  stationary  relative  to  each  other  while  the  receiving  antenna 
is  moved.  If  the  object  is  rotated  (as  in  the  case  of  Fig.l)  or  the  transmitter  is  moved  (as  in 
SAR),  then  the  object  scattering  function  can  no  longer  be  defined  uniquely,  and  the  3-d 


Fourier  transform  relation  between  the  corrected  far  field  and  the  "object  scattering  function" 
will  no  longer  be  exactly  satisfied 

In  spite  of  the  ambiguity  of  the  relationship  between  the  scattered  field  and  "object 
scattering  function"  as  stated  above,  microwave  images  were  reconstructed  by  Fourier 
transform  of  the  corrected  scattered  fields.  What  the  image  represented  and  what  the  image 
will  look  like  for  given  spectral  and  angular  windows  needed  more  investigation. 

In  the  following,  we  will  investigate  the  scattering  mechanism  of  the  object  and  the  pro¬ 
cedures  of  the  reconstruction  algorithm  to  interpret  what  the  image  represents  and  to  predict 
what  the  image  will  look  like  over  prescribed  angular  and  spectral  windows. 

ffl.  SCATTERING  PROPERTIES  OF  A  METALLIC  OBJECT 

For  an  object  large  compared  to  a  wavelength  of  the  incident  wave,  the  major  contribu¬ 
tions  to  the  scattered  field  are  specular  scattering  points,  edge  diffraction,  and  multiple  scatter¬ 
ing  points  [10].  The  properties  of  the  above  three  scattering  mechanisms  will  be  briefly  dis¬ 
cussed. 


1.  Specular  scattering  points 

Under  the  physical  optics  approximation,  the  asymptotic  expression  of  Eq.(l)  as  k  goes 
to  infinity  can  be  obtained  by  applying  the  stationary  phase  method  [11]  and  the  result  is 

A*  (*/„/,)  =  -rj|-  e~jkR'  £  -pj=  2n  (r/)xlHe  *  (5) 
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where  7*:'  are  the  vectors  such  that 


da' 


I 


r  ’  =  ts‘  ~  0  * 


(6) 


and 


5;  = 


d2[Or-it)-r']\ 


da' 


r '  =  r; 


(7) 


where  *77  is  derivative  with  respect  to  the  surface  curvature.  The  points  corresponding  to 
da 

the  solutions  of  Eq.(6)  are  called  stationary  points  or  equiphase  points  or  scattering  centers. 
The  absolute  value  of  the  term  n  ( T*  ')xlH  /  yjs~  is  called  the  scattering  strength  for  those 


scattering  centers.  It  is  seen  that  locations  of  the  scattering  centers  depend  on  directions  of  lr 
and  1,  as  well  as  the  shape  of  the  metallic  surface.  The  scattering  strength  for  each  scattering 
center  depends  on  the  local  properties  of  the  scattering  centers.  The  second  derivative  at  the 
stationary  point  is  proportional  to  the  curvature  of  that  point  The  above  analysis,  illustrates 
that  the  objects  we  would  be  dealing  with  in  high  frequency  radar  imaging  are  of  discrete 
form,  consisting  of  point  scattering  centers. 


2.  Conducting  wedge  with  finite  length 

The  diffracted  field  of  a  finite  edge  for  arbitrary  incidence  and  observation  or  diffraction 
angles  has  been  derived  [12],  where  the  concept  of  equivalent  electrical  current  and  equivalent 
magnetic  current  has  been  applied. 


Without  loss  of  generality,  consider  a  z-directed  finite  wedge  with  included  angle  (2-n  )n 
and  edge  length  L  as  shown  in  Fig.2.  The  origin  is  chosen  at  the  center  of  the  edge.  Let  the 
x-axis  be  directed  along  the  normal  to  the  edge  lying  on  the  front  face  of  the  wedge,  and  the 
y-axis  be  coincident  with  the  normal  to  this  face,  with  s  '  the  direction  of  incident  wave,  and 
S  the  direction  of  diffracted  wave.  The  angles  between  the  s  and  s  '  and  the  edge  are  9  and 
0  '  respectively.  The  angle  between  the  front  face  and  the  edge-fixed  plane  of  incidence  con¬ 
taining  the  vectors  S  '  and  f  is  0  ' .  The  corresponding  angle  for  the  edge-fixed  plane  of 
observation  is  0. 

The  incident  field  can  be  decomposed  into  TE  and  TM  components(or  parallel  and  per¬ 
pendicular  polarization)  with  respect  to  the  plane  of  incidence  while  the  diffracted  wave  can 
also  be  decomposed  into  TE  and  TM  components  with  respect  to  the  plane  of  observation. 
Denote  the  equivalent  electric  current  and  equivalent  magnetic  current  on  the  edge  as  I  and 
M.  The  expression  of  I  and  M  for  arbitrary  incidence  and  observation  angles  has  been  found 
in  [12].  For  the  special  case  of  backscattering  and  plane  conductor  (i.,e,  n=2),  the  expression 
of  I  and  M  can  be  simplified  as: 
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where  Hj  and  £,'  are  the  z  component  of  the  incident  fields  along  the  edge.  The  scattered 
fields  can  be  shown  to  be 
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3.  Multiple  reflection  points 

In  some  situations,  multiple  reflections  become  the  major  contributors  to  the  scattered 
field.  To  determine  the  positions  of  the  multiple  reflection  points,  the  ray  optics  approxima¬ 
tion  can  be  used.  The  ray  should  satisfy  Snell’s  law,  i.e.,  the  angle  of  incidence  should  be 
equal  to  angle  of  reflection.  Assume  now  that  the  multiple  reflection  points  have  been  found, 
then  the  scattered  field  due  to  a  set  of  two  multiple  reflection  points  may  be  approximated  by 
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where  /?2i  is  the  distance  between  the  two  scattering  points,  r2l  the  unit  vector  connecting 
from  scattering  point  1  to  point  2,  r*2 '  and  7*]  '  the  position  vectors  of  the  two  multiple 
scattering  points,  Sx  and  S2  the  second  derivative  at  points  1  and  2  as  defined  in  Eq.(7),  and 
n  (7*!  Oand  n  (7*2  ' )  the  normal  vectors  at  position  7*!  '  and  7*2  '  respectively. 


IV.  INTERPRETATION  OF  THE  RECONSTRUCTED  IMAGES 

The  reconstruction  method  used  in  microwave  diversity  imaging  systems  is  either  the 
back-projection  method  or  the  2-d  FT  method.  These  two  methods  have  been  proven  to  yield 
equivalent  results  [13].  The  back-projection  algorithm,  however,  provides  more  physical 
insight  to  the  image  formation  process  of  microwave  diversity  imaging.  The  two  basic  pro¬ 
cedures  in  the  back-projection  algorithm  are  first  to  obtain  filtered  range  profiles  and  second 
the  back-projection.  The  Fourier  transform  of  the  frequency  response  with  respect  to 


frequency  will  give  the  range  profiles  at  that  view  angle.  The  filtered  range  profile  is  the  FT 
of  the  product  of  the  frequency  response  and  a  frequency  weighting  function. 

The  range  profiles  due  to  different  scattering  mechanisms  can  be  obtained  by  taking  the 
Fourier  transform  of  the  respective  corrected  scattered  field  frequency  response.  Because  the 
scattered  fields  are  measured  only  over  a  finite  bandwidth  and  low  frequency  data  is  not  avail¬ 
able,  the  range  profile  obtained  by  FT  of  the  passband  data  will  be  different  from  that 
obtained  by  FT  of  the  infinite  bandwidth  data.  The  effect  of  passband  filter  will  be  discussed 
in  more  detail  in  the  later  section.  Where  the  scattering  centers  are  the  major  contributors,  the 
Fourier  transform  of  Eq.(5)  with  respect  to  frequency  will  give  peaks  at  the  differential  ranges 
of  the  scattering  centers  with  amplitude  equal  to  the  scattering  strength.  Where  an  edge  is  the 
major  contributor,  the  Fourier  transform  of  Eq.(10)  and  (11)  will  give  two  peaks  which  are  at 
the  differential  ranges  of  the  end  points  of  the  edge  with  an  amplitude  roughly  proportional  to 
sine  (kL  cos  0  ).  Where  the  multiple  reflection  points  are  the  major  contributors,  the  Fourier 
transform  of  Eq.(12)  will  give  peaks  at  a  range  equal  to  the  differential  range  of  the  multiple 
scattering  path.  To  a  first  order  approximation,  the  scattered  field  of  a  complex  object  can  be 
expressed  as  the  linear  superposition  of  the  scattered  field  from  discrete  scattering  centers, 
diffraction  from  the  "visible"  edges,  and  the  scattered  fields  from  any  present  multiple 
reflection  points.  Peaks  of  the  range  profiles,  therefore,  correspond  to  the  differential  ranges 
of  the  specular  points,  the  end  points  of  the  "visible"  edges,  and  the  differential  ray  path 
length  of  the  multiple  scattering  points. 

Once  the  range  profiles  are  known  for  sufficient  number  of  object  aspects,  the  next  step 
is  to  back-project  the  range  profiles  onto  each  pixel  of  the  image  plane.  When  the  object  is 
rotated  or  the  receiver  is  moved,  the  scattering  centers  will  move  to  those  points  which  satisfy 
Eq.(6)  with  scattering  strengths  given  by  Eq.(7).  The  differential  ranges  and  the  equivalent 
scattering  strength  of  the  end  points  of  edges  will  also  be  changed. 

The  implemention  of  the  back-projection  algorithm  should  satisfy  the  following  require¬ 
ments  [14]: 

(1) .  The  back-projection  line  should  contribute  to  those  pixels  which  it  intersects  and  to  no 

others. 

(2) .  The  contribution  to  the  backprojection  line  must  be  proportional  to  the  scattering 

strength. 


(3).  The  contributions  to  a  pixel  from  different  filtered  range  profiles  must  add  ’coherendy’ 

such  that  the  bipolar  nature  of  the  filtered  range  profiles  is  preserved. 

If  only  the  scattering  centers  are  considered,  after  back-projection,  the  image  recon¬ 
structed  can  be  interpreted  as  "collections  of  the  projected  scattering  centers  over  the  angular 
window".  To  predict  what  the  image  will  look  like  over  the  angular  window,  one  can  formu¬ 
late  the  surface  equations  of  the  target,  use  Eq.(6)  and  Eq.(7)  to  locate  the  positions  and  cal¬ 
culate  the  scattering  strength  of  the  scattering  centers  for  each  view  angle,  and  then  use  back- 
projection  to  reconstruct  the  images.  If  the  surface  of  the  object  is  smooth,  and  the  incremen¬ 
tal  rotation  angle  is  small,  then  the  image  will  be  "the  projected  outline  of  the  scattering 
centers  of  the  object  over  the  angular  window". 

For  an  edge  with  finite  length,  its  image  can  also  be  formulated  through  back-projection. 
To  simplify  the  discussion,  only  the  back-scattering  case  (i.e.,  when  S  -  -S'  in  Fig.  2)  is  con¬ 
sidered.  Let  L  be  the  length  of  the  wedge.  For  a  specific  aspect  angle,  the  positions  of  the 
peaks  in  the  range  profile  are  located  at  the  projections  of  the  end  points  along  the  line  of 
sight.  Amplitudes  of  the  peaks  arc  functions  of  6  and  $  and  are  proportional  to  sine 
(  kL  cos  0  ),  where  6,  <j>  are  defined  in  the  edge-fixed  coordinate  as  described  in  Fig.2.  After 
back-projection,  the  contribution  to  the  reconstructed  image  due  to  this  view  angle  will  be  two 
tines.  Each  of  these  tines  passes  through  one  of  the  end  points.  When  the  object  is  continu¬ 
ously  rotated,  0  and  <j>  will  also  be  changed.  The  contributions  to  the  reconstructed  image  due 
to  the  data  of  each  view  is  two  tines,  except  when  9  =  90°.  All  the  back-projection  tines  pass 
through  one  of  the  two  points,  which  means  the  brightness  of  the  end  points  is  intensified.  If 
the  angular  window  containes  a  view  so  that  9  =  90°  (the  incident  wave  is  normal  to  the 
edge),  then  the  back-projection  due  to  this  view  will  produce  a  bright  line  because  L  cos0  =0, 
and  sine  (  kL  cos  0  )  =  1,  and  this  bright  tine  will  overlap  with  the  edge.  The  edge,  then,  will 
appear  in  the  reconstructed  image  only  if  the  angular  window  contains  this  view;  otherwise, 
the  reconstructed  image  of  an  edge  will  be  only  two  bright  points. 

V.  EXAMPLES 

Some  examples  of  the  reconstructed  images  and  their  interpretations  are  given  next  In 
the  following  numerical  examples,  the  frequency  coverage  is  from  6  GHz  to  16  GHz.  In  the 
first  example,  we  show  the  dependence  of  the  image  on  the  data  acquisition  geometry  (mono¬ 
chromatic  and  bistatic  schemes).  Consider  an  infinite  metallic  cylinder  with  radius  20  cm.  The 


cylinder  axis  is  parallel  to  the  rotation  axis.  The  transmitting  and  receiving  antennas  have  the 
same  linear  polarization  which  is  parallel  to  the  axis  of  cylinder.  The  reconstructed  images 
are  obtained  by  rotating  the  target  over  360°  and  by  changing  the  bistatic  angle  of  the 
transmitter  and  receiver  over  360°  .  The  scattering  center  of  a  cylinder  in  the  monostatic  case 
is  at  a  point  with  its  surface  normal  vector  parallel  to  the  incident  wave,  while  the  scattering 
center  in  the  bistatic  case  is  at  a  point  with  its  surface  normal  vector  parallel  to  the  line 
bisecting  the  directions  of  incidence  and  observation.  The  scattering  strength  of  the  scattering 
center  for  the  bistatic  case  is  proportional  to  cos  a,  where  a  is  the  bistatic  angle.  The  collec¬ 
tion  of  the  scattering  centers  is  a  circle  for  the  monostatic  case  and  is  a  semicircle  for  the  bis¬ 
tatic  case.  The  positions  of  the  scattering  centers  on  the  cylinder  by  rotating  the  object  (or 
equivalently  by  rotating  the  transmitter/receiver  around  the  object)  and  by  rotating  the  receiver 
only  and  the  sketchs  of  the  back-projection  for  these  two  cases  are  shown  in  Fig.3(a)  and  3(b) 
respectively.  The  numerical  simulated  images  are  shown  in  Fig.3(c)  and  3(d). 

In  another  example,  we  show  the  effect  of  the  angular  window  on  the  reconstructed 
image.  The  geometry  of  a  plate  is  shown  in  Fig.4(a).  The  inclination  angle  between  the  plate 
and  the  rotational  axis  is  a.  The  polarizations  of  the  transmitting  and  receiving  antennas  are 
right-hand  circularly  polarized  and  left-hand  circularly  polarized  respectively,  which  means 
that  the  image  is  formed  from  the  co- polarized  scattered  field.  The  scattered  fields  of  the  hex¬ 
agonal  plate  are  assumed  to  be  contributed  by  the  six  edges.  The  images  reconstructed  from 
angular  window  1  (  0°,  180°  ),  window  2  (  0°,  110°  ),  window  3  (  70°,  180°  ),  and  window  4 
(  0°,  87  °)  are  shown  in  Fig.4(b)  to  4(e)  respectively.  In  window  1,  each  edge  becomes  nor¬ 
mal  during  rotation  of  the  object  to  a  particular  incident  vector,  therefore,  all  edges  appear  in 
the  image  of  Fig.4(b).  In  window  2,  there  is  no  incident  vector  normal  to  edge  1  and  4,  hence 
these  two  edges  disappear  in  Fig.4(c);  but  the  end  points  have  been  intensified.  The  image  of 
Fig.4(d)  and  4(e)  can  be  interpreted  in  the  same  manner. 

In  another  example,  we  show  the  effect  of  multiple  scattering  on  the  image  and  the  effect 
of  different  bistatic  angles  of  data  acquisition.  The  object  consists  of  two  parallel  separated 
long  cylinders  as  shown  in  Fig.5(a).  The  radius  of  each  cylinder  is  4  cm,  and  the  distance 
between  the  two  cylinders  is  60  cm.  Both  the  transmitting  and  receiving  antennas  are  right- 
hand  circularly  polarized  which  mean  that  the  image  is  formed  from  the  cross-polarized  scat¬ 
tered  field.  In  this  transmitting  and  receiving  antenna  system,  multiple  scattering  phenomenon 
is  more  pronounced,  because  the  specular  reflection  will  be  rejected  by  the  antenna  arrange¬ 
ment  Experimental  and  numerical  results  concerning  multiple  scattering  on  this  structure  has 


been  discussed  in  [IS].  Series  expansion  solutions  were  used  to  calculate  the  scattered  fields, 
and  lengthy  computing  time  is  needed  to  generate  the  field  data.  However,  under  high  fre¬ 
quency  approximation,  the  multiple  reflection  points  for  a  given  set  of  aspect  angle  and  bis¬ 
tatic  angle  can  be  found  by  ray  optics.  The  detail  of  the  derivation  is  given  in  the  appendix. 
The  differential  path  length  versus  rotation  angle  $  and  the  numerically  reconstructed  images 
with  bistatic  angle  a  =  0°  ,  a  =  16°,  and  a  =  40°  are  shown  in  Fig.  5(c),(d)  and  (e) 
respectively.  For  a  bistatic  case,  there  are  two  sets  of  multiple  reflection  points.  The 
differential  path  lengths  of  these  two  sets  are  different  and  is  a  function  of  the  bistatic  angle, 
which  explains  the  difference  between  these  images.  It  is  seen  that  multiple  reflections  has 
produced  artifact  in  the  image,  that  does  not  correspond  to  the  physical  shape  of  the  object. 

An  experimental  example  of  the  reconstructed  image  of  a  complex  object  is  given  next. 
The  test  object,  a  metalized  100:1  scale  model  of  a  B-S2  aircraft  with  79  cm  wing  span  and 
68  cm  long  fuselage  was  mounted  on  a  computer-controlled  positioner  situated  in  an  anechoic 
chamber  enviroment.  201  equal  frequency  steps  covering  the  6.1  to  17.5  GHz  range  were  used 
to  obtain  the  frequency  response  of  the  object.  The  transmitting  and  receiving  antennas  are 
both  right-hand  circularly  polarized.  The  reconstructed  image  of  the  test  object  using  data  col¬ 
lected  in  an  angular  window  of  90°  extending  from  head-on  to  broadside  in  128  looks  is 
shown  in  Fig.6.  The  intensity  or  brightness  of  each  pixel  in  the  image  is  proportional  to  the 
average  scattering  strength  of  that  portion  over  this  angular  window.  This  fact  indicates  that 
microwave  diversity  imaging  can  be  applied  to  the  study  of  radar  cross  section  (RCS) 
management 

Over  this  angular  window,  there  exist  certain  rotation  angles  such  that  the  edges  of  the 

right  wing  and  right  tail  are  perpendicular  to  the  bisection  line  or  of  the  angle  between  the 

* 

lines  joining  the  transmitter  and  receiver  to  the  rotation  center  of  the  target.  However,  there  is 
no  such  rotation  angle  that  the  bisecting  line  is  perpendicular  to  the  edges  of  the  left  wing  and 
tail.  Therefore,  the  edges  of  the  right  wing  and  right  tail  can  be  seen  in  the  reconstructed 
image,  while  the  appearance  of  the  left  tail  in  the  image  shows  only  the  end  points  of  the 
edges  as  shown  in  Fig.  6(b). 

If  the  angular  window  is  reduced  to  72°  from  head-on  as  shown  in  Fig.6(c),  the  recon¬ 
structed  image  becomes  as  shown  in  Fig.6(d).  Over  this  reduced  angular  window,  there  is  no 
specular  reflection  from  the  fuselage  and  engines  and  fuel  tank  on  the  right  side.  This  explain 
the  difference  between  Fig.6(b)  and  6(d).  From  the  above  experimental  examples,  we  can 


conclude  that  the  microwave  images  represent  the  "collection  of  the  projected  scattering 
centers,  those  visible  edges  which  are  perpendicular  to  the  bisection  line,  and  end  points  of 
the  visible  edges,  over  the  precribed  angular  window". 

The  above  examples  show  that  the  reconstructed  image  is  dependent  on  the  data  acquisi¬ 
tion  parameters  (monostatic,  bistatic,  polarization,  frequency  range,  angular  window)  and 
hence  care  must  be  taken  for  proper  interpretation  of  the  image. 

VL  ALTERNATIVE  INTERPRETATION  OF  EDGE  ENCHANCEMENT  PROPERTY 
OF  MICROWAVE  IMAGE 

As  mentioned  in  Sec.n  and  ref  [2,3],  the  scattering  object  function  defined  is  related  to 
the  geometrical  shape  of  the  object  SAR  images  or  the  images  obtained  by  microwave  diver¬ 
sity  imaging,  however,  are  not  solid  filled  images,  but  are  contoured  or  edge-enchanced 
images  as  illustrated  in  the  previous  examples.  This  property  can  be  explained  from  the  image 
formation  principles  described  in  the  previous  sections,  however,  it  can  also  be  explained  from 
another  point  of  view. 

In  SAR,  in  stead  of  the  object  scattering  function,  a  target  reflectivity  p(x,  R)  is  defined 
in  terms  of  the  range  coordinate  R  and  cross-range  coordinate  x  measured  from  the  center  of 
the  synthetic  aperture  [9].  The  effect  of  the  finite  range  resolution  on  the  nature  of  SAR 
image  was  discussed  [9].  The  range  processed  retur.  with  finite  resolution  was  derived  using 
successive  integration  by  parts,  and  a  conclusion  was  made  that  the  finite  range  resolution 
accentuates  edges  of  the  image  and  eliminates  portions  of  the  image  which  have  constant 
intensity  over  long  distances.  This  conclusion  can  also  be  easily  obtained  by  analyzing  the 
scattered  fields  from  the  frequency  domain  over  a  finite  bandwidth.  By  using  the  same  exam¬ 
ple  as  in  [9],  the  range  reflectivity  is  shown  in  Fig.  7.  The  scattered  field  which  is  the  Fourier 
transform  of  the  range  reflectivity  can  be  obtained  by  successive  differentiation  of  the  range 
reflectivity  with  respect  to  range.  From  the  fundamental  Fourier  transform  theorem,  we  have 
the  following  FT  pairs: 

/(*)< — *F(o>) 
juF((i>) 

/"(R)A— 5-m2  F(  (0  ) 

The  FT  of  $(R  -  R0  )  is  The  FT  of  Fig.  7  can  be  easily  expressed  as 
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F(g>)  =  flj- 


,-J*  1“ 


,-jR  i<o 


where  a2  and  <23  are  the  slope  of  the  segment  2  and  3  respectively.  If  the  ratio  of  the  highest 
frequency  to  the  lowest  frequency  of  the  measurement  system  is  not  much  larger  than  1,  then 
the  denominator  can  be  approximated  by 

J_  =  1  =  JL  (16) 

to  0.5(aw  +  gw  )  <o0 

-4  =  - T  <H> 

or  (o; 

The  estimated  target  reflectivity  (or  range  profile)  will  be 

f(R)  =  [  a&R-RO  +  a26(R-R2)  ~  a&R-Rj)  -  a&R-RJ]  *  sine  (  ~  )  (18) 

The  above  simple  example  and  analysis  show  the  range  profile  estimated  from  pass  band 
data  accentuate  the  discontinuity  of  the  reflectivity  function  or  object  scattering  function.  From 
this  point  view,  we  can  also  conclude  that  the  reconstructed  image  represents  the  collection  of 
the  discontinuities  of  the  reflectivity  function  over  the  prescribed  window. 

Although  discontinuity  of  the  reflectivity  distribution  function  can  be  used  to  understand 
and  interpret  the  image  formation  of  the  microwave  diversity  imaging,  high  frequency  asymp¬ 
totic  approximation  can  give  better  understanding  of  the  scattering  mechanism.  Range  resolu¬ 
tion  of  SAR  or  ISAR  is  usually  limited  by  the  pulse  width  of  the  transmitted  signal  or  the 
bandwidth  of  the  chirp  signal.  The  range  resolution  of  microwave  diversity  imaging  is  also 
limited  by  the  bandwidth  of  the  measuring  system.  However,  under  high  frequency  assump¬ 
tion,  the  scattered  fields  outside  the  available  bands  can  be  extrapolated  and  better  resolutions 
can  be  achieved  [16].  The  applicability  of  extrapolations  is  based  on  the  fact  that  the  objects 
dealt  with  in  high  frequency  radar  imaging  are  of  finite  extent  and  discrete  form,  consisting  of 
scattering  centers  and  visible  edges.  Therefore,  the  understanding  of  the  scattering  mechanism 
and  reconstruction  algorithm  is  more  helpful  to  interpret  and  to  predict  the  images. 

m  CONCLUSION 

The  principle  of  image  formation  of  microwave  diversity  imaging  systems  has  been 
understood  through  investigating  the  scattering  mechanism  of  the  object  and  the  procedures  of 
the  reconstruction  algorithm.  Some  numerical  and  experimental  results  concerning  the  images 
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of  a  metallic  object  obtained  by  prescribed  spectral  and  angular  windows  have  been  explained 
satisfactorily  and  predicted  correctly.  The  edge-enhanced  nature  of  microwave  diversity  imag¬ 
ing  has  been  interpreted  by  a  rather  simple  analysis.  These  successful  interpretations  are  fun¬ 
damental  to  research  in  several  areas,  for  example,  radar  cross  section  management  studies 
[171,  target  recognition.  Furthermore,  it  gives  directions  to  design  some  techniques  to  obscure 
the  reconstructed  images,  which  are  the  research  topics  under  investigation. 

APPENDIX 

Consider  two  parallel  identical  cylinders  with  radius  a,  separated  by  a  distance  2r0  seated 
on  a  rotating  pedestal  and  illuminated  by  a  plane  wave  as  shown  in  Fig.  5(a).  The  cross  sec¬ 
tion  of  that  arrangement  is  shown  in  Fig.Al.  The  line  connecting  the  centers  of  the  cylinders 
passes  through  the  rotating  center  and  makes  an  angle  8  with  the  respect  to  the  rotation  axis. 
The  polar  coordinates  and  rectangular  coordinates  of  the  points  on  the  surfaces  of  these  two 
cylinders  with  respect  to  each  cylinder  center  are  ( a ,  y  i0.  and  (a,  fcO.  C*2\  yz) 

respectively.  Let  the  bistatic  angle  between  the  transmitter  and  receiver  be  a,  the  unit  vectors 
in  the  direction  of  transmitter  and  receiver  be  /,  and  lr  respectively.To  find  the  stationary 
points  of  the  multiple  reflection,  geometrical  optics  will  be  applied.  The  reflection  law  requires 
that  the  angle  between  the  incident  ray  and  normal  line  must  be  equal  to  the  angle  between 
the  reflected  ray  and  the  normal  line. 

The  stationary  points  of  multiple  reflection  on  the  respective  cylindrical  surfaces  for  a 
given  set  of  (r0,  a,  0,  a)  are  to  be  determined.  Assume  the  incident  ray  first  hit  the  cylinder 
2,  reflects  to  cylinder  1,  and  then  bounces  back  to  the  receiver  as  shown  in  Fig.Al  (a).  The 
notation  of  the  angles  Vi*  V2>  a*  9.  0i'.  W  are  defined  in  Fig.  Al(a).  It  can  be  shown  the 
following  relationship  must  be  satisfied  in  order  to  satisfy  the  reflection  law. 

2r0sin6  -  a  [cosCfo'  -  a/2)  +  sinfo'] 

— 2 — - — - -  =  -cot(2<|>2'  -  a/2  )  (A.l) 

2rgcos0  +  a  [sm(<{>2  -  a/2)  -cost^T  2 

4>r  =  -  90°  -  a/2  (A.2) 

where  $2  is  restricted  to  be 

270°  £  S  360°  when  -90°  £  9  £  90° 

180°  £  02  *  270°  when  90°  £  0  £  270° 


For  a  given  set  of  parameters  {tq,  a,  9,  a  },  the  azimuth  angle  $2  anti  therefore  0j'  can  be 
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detcmrincd  by  solving  Eq.(A.l)  and  (A.2). 

The  other  case  is  that  the  incident  ray  first  hit  cylinder  1  .reflect  to  cylinder  2,  and  then 
bounce  back  to  the  receiver.  The  geometry  of  this  situation  is  shown  in  Fig.Al(b).  The  sta¬ 
tionary  points  must  satisfy  the  following  relationship 


<h"  *  -  90°  +  ct/2 

2rosinO  -  a  [cos(  +  a/2)  +  sin^'T 
2rocos0  +  a  [sin(  fa"  +  a/2  )  -  cos^"] 


-cot(202,#  +  ot/2) 


(A.3) 


(A.4) 


From  the  above  derivation,  it  can  be  seen  that  there  are  two  pairs  of  multiple  reflection  points 
for  bistatic  system  and  only  one  multiple  reflection  points  for  monostatic  case. 
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FIGURE  CAPTION 

Fig  1.  A  simplified  microwave  diversity  imaging  arrangement 

Fig.2  Geometry  of  a  finite  wedge  diffraction  problem 

Fig.3  Positions  of  the  scattering  centers  of  a  cylinder  and  the  sketch  of 
back-projection  for  the  case  of  (a),  rotating  the  transmitter/receiver, 
and  (b).  rotating  the  receiver  only;  over  360°. 

(c) .  the  numerically  reconstructed  image  of  (a); 

(d) .  the  numerically  reconstructed  image  of  (b). 

Fig.4  (a).  Geometry  of  a  hexagonal  plate  and  the  measurement  arrangement 
Image  reconstructed  from  different  angular  windows: 

(b).  window  1  (0°  ,  180°),  (c).  window  2  (0°,  110°), 

(d).  window  3  (70°,  180°),  (e).  window  5  (0°,  87°). 

Fig.5  (a).  A  dual-cylinder  object, 

and  plots  of  the  its  numerically  calculated  differential  range  of 
multiple  reflection  path  versus  rotation  angle  and  simulated  images 
obtained  by  using  different  bi-static  angle  equal  to 
(b).  0°,  (c).  16°,  (d).  40°. 

Fig.6  (a).  Sketch  of  a  B-52  test  object  aiiplane  and  the  angular  window. 

(b) .  Reconstructed  image  of  (a). 

(c) .  Reconstructed  image  of  (c). 

Fig.7  An  example  of  target  range  reflectivity  together  with  its  successive 
differentials. 

Fig.Al  Multiple  reflections  of  a  two-cylinder  object. 

(a) .  The  incident  wave  reaches  the  left  cylinder  first,  and  is  then 

reflected  back  by  the  right  cylinder. 

(b) .  The  incident  wave  reaches  the  right  cylinder  first  and  is  then 
reflected  back  by  the  left  cylinder. 
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Fig.4  (a).  Geometry  of  a  hexagonal  plate  and  the  measurement  arrangement. 
Image  reconstructed  from  different  angular  windows: 

(b).  window  1  (0°  ,  180°),  (c).  window  2  (0°,  110°), 

(d).  window  3  (70°,  180°),  (e).  window  5  (0°,  87°). 
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Fig.5  (a).  Imaging  geometry  of  a  dual-cylinder  object  Plots  of  the  numerically  cal¬ 
culated  differential  range  of  multiple  reflection  path  versus  rotation  angle  and  simu¬ 
lated  images  obtained  by  using  different  bistatic  angles  equal  to  (b).  0°,  (c).  16°, 
(d).  40°. 
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RADAR  CROSS  SECTION  MANAGEMENT  STUDIES 
EMPLOYING  MICROWAVE  DIVERSITY  IMAGING 

HsuehrJyh  Li.  NabU  H.  Farhat,  Yuhsyen  Shot 

The  Moore  School  of  Electrical  Engineering 
University  of  Pennsylvania 
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L  INTRODUCTION 

Radar  cross  section  (RCS)  is  a  measure  of  the  equivalent  size  of  a  target  as  seen  by  the 
radar  [1].  RCS  is  a  function  of  aspect  angle,  frequency,  and  polarization.  For  certain  applica¬ 
tions,  it  is  desirable  to  enhance  the  RCS  over  some  specific  range  of  aspect  angles.  For  exam¬ 
ple,  strong  return  signals  are  required  to  track  a  missile  or  an  aircraft  in  flight  In  contrast  for 
other  applications  it  is  often  desirable  to  reduce  or  minimize  the  RCS  over  specified  spectral 
and  angular  windows  so  the  target  is  less  likely  to  be  detected. 

RCS  management  is  a  general  term  for  obtaining  the  RCS  of  a  scattering  object  by  mani¬ 
pulating  the  distribution  and  strengths  of  the  hot  spots  or  flare  spots  of  a  target  over 
prescribed  spectral  and  angular  windows  and  states  of  polarization.  Hot  spots  represent  all 
those  areas  that  have  major  contributions  to  the  received  scattered  fields.  There  are  several 
ways  to  determine  the  site  of  hot  spots  [2].  An  experienced  engineer  can  roughly  point  out 
the  possible  locations  of  the  hot  spots;  a  computer  code  can  be  used  to  estimate  the  locations 
of  the  scattering  centers  and  their  scattering  if  the  equations  that  describe  the  surface  of  the 
object  can  be  formulated.  However,  all  such  estimates  need  experimental  verification. 

Microwave  diversity  imaging  systems  using  angular,  wavelength,  and  polarization  diver¬ 
sity  have  been  used  to  obtain  images  of  metallic  objects  with  nearly  optical  resolution  [3]. 
The  pixel  intensity  or  brightness  of  the  reconstructed  image  is  proportional  to  the  mean 
scattering  strength  at  the  corresponding  position  averaged  over  the  specified  angular  window 
and  spectral  window.  This  fact  suggests  microwave  diversity  imaging  can  be  employed  to 
determine  the  hot  spots  and  to  study  RCS  management 


Two  steps  are  usually  involved  in  the  study  of  RCS  management  First  the  hot  spots 
must  be  identified.  Certain  techniques  are  then  utilized  to  manipulate  the  hot  spots.  Fine 
resolution  of  object  detail  in  the  reconstructed  image  is  essential  for  locating  the  hot  spots  of 
an  object  over  the  prescribed  spectral  and  angular  windows  and  the  given  states  of  polariza¬ 
tion  of  the  transmitter  and  receiver.  The  range  resolution  of  a  frequency  diversity  imaging 
system  is  proportional  to  the  bandwidth  of  the  measurement  system.  If  the  bandwidth  is  too 
narrow  to  give  acceptable  resolution,  the  algorithm  for  Ate  extrapolation  of  data  available  in 
multiple  restricted  bands  or  in  a  single  narrow  band  developed  in  [4]  can  be  applied  to  locate 
the  hot  spots. 

Previously,  radar  workers  usually  treated  the  RCS  from  the  detection  perspective.  There¬ 
fore,  techniques  of  reducing  RCS  were  developed  to  reduce  the  signal  to  noise  ratio  (S/N)  so 
that  the  target  can  more  easily  elude  detection  from  radar.  However,  as  will  be  shown  later  a 
low  S/N  is  not  a  sufficient  criterion  for  eluding  detection  when  microwave  diversity  imaging 
systems  are  employed  where  signals  from  many  broadband  sensors  are  combined  to  form  an 
image.  In  this  paper  we  will  consider  therefore  the  RCS  management  from  the  imaging  point 
of  view,  that  is,  we  consider  principles  and  techniques  for  making  the  image  unrecognizable 
by  radar. 

In  Sec.  II  image  formation  in  microwave  diversity  imaging  will  be  summarized.  A  new 
term,  “diaphanization”,  will  be  introduced  in  Sec.  m.  We  will  also  describe  how  to  use 
microwave  diversity  imaging  to  diaphanize  an  object  in  that  section.  Effects  of  Gaussian  noise 
on  the-  reconstructed  image  will  be  discussed  in  Sec.  IV.  Some  fundamental  concepts  for 
obscuring  the  image  will  be  proposed  in  Sec.  V.  Experimental  results  on  different  diaphaniza¬ 
tion  techniques  will  be  given  in  Sec.  VL 

II.  IMAGE  FORMATION  IN  MICROWAVE  DIVERSITY  IMAGING 

The  simplified  microwave  diversity  imaging  arrangement  is  shown  in  Fig.  1 .  A  metallic 
object  seated  on  a  rotated  pedestal  is  illuminated  by  a  plane  wave.  For  each  aspect  angle  a 
set  of  pulses  at  different  frequencies  are  transmitted  and  their  echoes  are  received.  The  object 
is  then  rotated  and  the  measurement  is  repeated  to  obtain  the  multiaspect  stepped  frequency 


response  of  the  scattering  object  Analysis  of  microwave  diversity  imaging  has  been  discussed 
in  [5]. 

The  reconstruction  method  used  in  microwave  diversity  imaging  systems  is  either  the 
back-projection  method  or  the  2-d  FT  method  These  two  methods  have  been  proven  to  yield 
equivalent  iesults[6].  The  back-projection  algorithm,  however,  provides  more  physical  insight 
into  the  image  formation  process  in  microwave  diversity  imaging  [12].  The  two  basic  pro¬ 
cedures  in  die  back-projection  algorithm  are  first  obtaining  filtered  range  profiles  and  then 
back-projecting  the  range  profiles  onto  each  pixel  of  an  image  plane.  The  range  profile  of  a 
given  aspect  is  obtained  by  Fourier  transforming  the  frequency  response  at  that  viewing  angle 
while  the  filtered  range  profile  is  the  FT  of  the  product  of  the  frequency  response  and  an 
appropriate  frequency  weighting  function. 

At  high  frequencies  the  major  contributions  to  the  scattered  field  of  a  complex  shaped 
metallic  object  are  specular  scattering  points,  edge  diffractions,  multiple  reflection  points, 
creeping  waves,  and  surface  traveling  waves  [7].  Properties  of  the  range  profile  for  each 
scattering  mechanism  has  been  discussed  in  [5].  Locations  of  peaks  in  a  range  profile  reflect 
the  differential  ranges  of  specular  points,  end  points  of  the  “visible”  edges,  differential  ray 
path  length  covering  the  multiple  reflection  points,  and  /  or  the  ray  path  length  of  the  creeping 
waves  and  surface  traveling  waves. 

Once  the  range  profile  far  each  aspect  angle  is  determined,  the  next  step  is  to  back- 
project  the  range  profile  onto  each  pixel.  If  we  only  consider  the  scattering  centers,  then  after 
back-projection,  the  reconstructed  image  can  be  interpreted  as  the  projected  outline  of  the 
specular  points  of  the  object  for  the  prescribed  angular  window,  which  usually  reflects  the 
projective  outline  of  the  visible  portions  of  the  object  shape.  If  an  edge  with  finite  length  is 
involved,  the  microwave  image  of  an  edge  will  be  two  bright  points  whose  locations 
correspond  to  the  end  points  of  the  edge  unless  the  normal  aspect  ( i.e.,  incident  wave  is  nor¬ 
mal  to  the  edge)  is  included  in  the  angular  window;  only  then  will  an  edge  appear  in  the 
image  [5]. 


HL  DIAPHANIZATION 


According  to  the  Webster  dictionary,  “diaphanous”  means: 

1.  so  fine  or  gauzy  in  texture  as  to  be  transparent  or  translucent 

2.  vague  or  indistinctive,  airy. 

Accordingly,  we  introduce  a  new  term  “diaphanization”,  which  can  be  used  to  describe: 

1.  All  techniques  of  reducing  RCS  in  order  to  render  a  target  invisible  to  radar  detection. 

2.  All  techniques  of  obscuring  the  image  so  that  the  image  is  vague,  indistinctive  or  not 
recognizable  by  an  imaging  radar 

Therefore,  diaphanization  is  a  technique  which  will  influence  either  detection  or  image 
formation  and  recognition. 

From  the  detection  perspective,  diaphanization  will  generally  lower  the  target’s  signal  to 
noise  ratio  (S/N),  thus  decreasing  the  likelihood  of  detection  [1].  Therefore,  it  is  desirable  to 
minimize  the  RCS.  From  the  imaging  point  of  view,  however,  low  S/N  is  not  a  fatal  problem 
for  microwave  diversity  imaging  systems.  It  will  be  shown  later  that  microwave  diversity 
imaging  is  very  robust  to  Gaussian  noise,  even  though  die  noise  is  stronger  than  the  signal. 
Therefore,  techniques  that  distort  the  reconstructed  image  to  make  it  unrecognizable  should  be 
considered.  In  the  following,  we  will  demonstrate  how  to  apply  microwave  diversity  imaging 
to  diaphanize  an  object  and  to  explain  the  motivation  for  studying  image  distortion. 

The  procedure  of  utilizing  microwave  diversity  imaging  to  diaphanize  a  target  for 
prescribed  spectral  and  angular  windows  is  as  follows: 

1.  Obtain  the  image  of  the  target  before  diaphanization. 

2.  Identify  the  hot  spots  which  are  the  bright  portions  in  the  image. 

3.  Use  certain  techniques  to  diaphanize  the  object,  i.e.  to  minimize  reflection  from  the  hot 
spots. 

4.  Obtain  the  image  of  the  target  following  diaphanization. 


5.  Compare  the  images  before  and  after  diaphanizadon. 

Techniques  of  reducing  the  RCS  can  be  grouped  into  three  classes  [7]:  covering  the 
object  with  absorbing  material  or  substituting  metallic  parts  with  resistive  or  composite 
material;  shaping  the  target;  and  impedance  loading.  The  first  method  is  based  on  absorption 
of  the  incident  radiation.  The  second  method  is  based  on  redirecting  the  incident  waves,  Le., 
scattering  the  incident  illumination  into  as  wide  a  solid  angle  as  possible  or  into  a  direction 
away  from  the  probing  system.  The  third  method,  also  known  as  passive  cancellation,  is 
based  upon  die  concept  to  introduce  an  echo  source  whose  amplitude  and  phase  can  be 
adjusted  to  cancel  another  echo  source  [7].  In  this  paper,  we  will  give  examples  showing  how 
microwave  diversity  imaging  is  used  in  RCS  management  and  diaphanization  studies  of  a 
complex  object 

In  the  following  experimental  example,  201  equal  frequency  steps  covering  the  6.0  to 
17.5  GHz  range  were  used  to  obtain  the  frequency  response  of  the  object  for  different  aspect 
angles.  The  magnitude  of  the  scattered  fields  are  measured  relative  to  a  reference  cylinder. 
The  RCS  defined  in  this  paper  is  the  square  of  die  measured  relative  field  and  therefore  the 
unit  of  RCS  is  a  relative  level  but  not  an  absolute  value.  Co-polarized  waves  refer  to  the 
waves  created  when  both  transmitting  and  receiving  antennas  have  the  same  sense  of  circular 
polarization;  while  cross-polarized  waves  are  the  waves  created  when  the  transmitting  antenna 
and  receiving  antenna  have  an  opposite  sense  of  circular  polarization. 

The  test  object,  a  metalized  1:100  scale  model  of  a  B-52  airplane  with  79  cm  wing  span 
and  68  cm  long  fuselage  was  mounted  on  a  computer-controlled  positioner  situated  in  an 
anechoic  chamber  environment  A  sketch  of  the  B-52  and  the  reconstructed  image  using  data 
collected  in  an  angular  window  of  90°  extending  from  head-on  to  broadside  in  128  looks  with 
cross-polarized  waves  are  shown  in  Figs.  2(a)  and  (b)  respectively.  The  intensity  or  brightness 
of  each  pixel  is  proportional  to  the  average  scattering  strength  of  that  portion  over  the  angular 
and  spectral  windows. 

The  bright  spots  are  the  hot  spots  of  the  target.  To  reduce  the  RCS,  we  cover  the  hot 
spots  with  broad-band  absorbers  (Emerson  and  Cumming  AN-72).  Portions  of  the  B-52  model 
covered  with  broad-band  absorber  are  indicated  by  the  dark  lines  in  the  sketch  of  Fig.  2(c). 


These  parts  consists  of  the  right  hand  sides  of  the  fuselage,  all  edges  of  the  right  wing,  and 
the  right  sides  of  the  two  right  engines  and  fuel  tank.  Hie  RCS  of  an  object  is  usually 
reduced  after  diaphanization,  but  it  is  not  reduced  for  every  frequency  because  the  scattered 
field  is  a  result  of  interference  and  is  a  function  of  frequency.  Therefore,  it  is  not  proper  to 
measure  the  effectiveness  of  a  diaphanization  on  RCS  at  a  specific  frequency  if  a  prescribed 
spectral  range  is  of  interest.  We  have  adopted  a  more  reasonable  measure  —  that  of  compar¬ 
ing  die  mean  scattered  fields  before  and  after  diaphanization.  The  mean  RCS  is  obtained  by 
averaging  the  RCS  at  each  frequency  over  the  prescribed  spectral  window  for  a  given  aspect 
angle.  The  mean  RCS  versus  aspect  angle  before  and  after  diaphanization  using  cross- 
polarized  waves  are  shown  in  Figs.  2(d).  From  these  plots,  we  find  that  the  RCS  has  been 
reduced  for  most  aspect  angles.  The  reconstructed  image  after  diaphanization  are  shown  in 
Figs.  2(e).  The  images  before  and  after  diaphanization  have  the  same  dynamic  display  range. 
One  can  see  that  diaphanization  reduces  the  brightness  of  the  hot  spots.  However,  if  we  lower 
the  dynamic  display  range  of  the  post  diaphanization  image,  the  image  becomes  as  shown  in 
Figs.  2(f).  Clearly,  the  nature  of  the  target  is  still  recognizable.  From  a  detection  point  of 
view,  covering  the  hot  spot  portions  did  reduce  the  RCS.  However,  from  the  imaging  perspec¬ 
tive,  the  image  is  still  recognizable. 

IV.  EFFECT  OF  GAUSSIAN  NOISE  ON  THE  RECONSTRUCTED  IMAGE 

It  is  known  that  low  signal  to  noise  ratio  (S/N)  will  decrease  detection  probability  and 
increase  the  false  alarm  rate  [1].  For  a  practical  radar  detection  system,  S/N  above  15  db  is 
usually  required.  It  is  also  known  that  coherent  summation  will  produce  coherent  gain  or 
superposition  gain.  That  is,  if  signals  contaminated  by  Gaussian  noise  are  coherently  summed 
n  times,  the  resultant  S/N  will  increase  n  times  [1].  The  reconstruction  algorithm  of 
microwave  diversity  imaging  involves  obtaining  the  filtered  range  profiles  and  back-projecting 
them  in  the  image  plane  [12].  Both  procedures  involve  coherent  summations  and  hence  result 
in  coherent  gain.  Therefore,  the  robustness  of  microwave  diversity  imaging  to  Gaussian  noise 
is  expected. 


Microwave  diversity  imaging  usually  uses  stepped  frequency  response  measurement  For 
each  frequency  step,  the  bandwidth  of  the  receiver  is  kept  the  same.  The  standard  deviation  of 
the  Gaussian  noise  at  each  frequency  step  is  therefore  the  same. 

Under  die  high  frequency  approximation,  at  a  specific  aspect  the  field  scattered  from  a 
metallic  object  can  be  expressed  as  the  supeiposition  of  the  scattered  fields  of  N  discrete 
scattering  centers  [4]  such  that  the  scalar  fields  at  frequency  fm  is  given  by 

Es(fm)  -  exp(/2* /mr<)  fl  Zfm  Zf2  (1) 

i- rl 

where  a,-  is  the  scattering  strength  of  the  i-th  scattering  center,  rt  is  the  differential  range  of 
the  i-th  scattering  center  with  respect  to  the  reference  point,  /t  and  /2  are  the  lowest  and 
highest  frequencies  of  the  measurement  system  respectively,  Np  is  the  number  of  frequency 
points.  Assume  the  noise  power  level  of  the  receiver  at  each  frequency  step  is  N,  then  the 
S/N  due  to  the  i-th  point  scatter  is  la;  I2  /  N.  Assume  the  point  scatterers  are  well  separated 
apart  so  that  at  range  r(-  the  contributions  due  to  die  side-lobe  of  other  scatterers  can  be 
neglected.  The  signal  (sampled  value  at  rt  in  the  FT  domain  or  range  profile)  to  noise  ratio 
will  become  Np  \at\2  IN  since  the  Fourier  transform  is  a  process  of  coherent  summation. 
Therefore  Np  is  the  coherent  gain  or  superposition  gain  for  this  example. 

Back-projection  is  a  process  of  coherent  summation  of  the  contributions  of  the  back- 
projection  lines  obtained  from  the  filtered  range  profile  of  each  aspect  angle  at  each  image 
plane  pixel.  This  additional  coherent  process  will  make  the  reconstructed  image  more  robust 
to  the  Gaussian  noise.  If  the  object  consists  of  isotropic  point  scatterers,  i.e.,  the  scattered 
fields  of  these  point  scatterers  are  isotropically  distributed,  then  the  overall  coherent  gain  will 
be  NaNp  if  data  from  Na  aspect  angles  are  used  to  reconstruct  the  image.  In  practice  for  a 
metallic  object  for  example,  the  locations  and  scattering  strengths  of  the  scattering  centers  and 
end  points  of  edges  will  change  or  migrate  when  the  object  changes  aspect  (rotation)  during 
data  acquisition.  Therefore,  the  signal  (image  intensity  at  a  pixel)  to  noise  ratio  in  the  image 
space  cannot  be  expressed  in  the  above  simple  form. 

The  following  example  illustrates  the  effect  of  Gaussian  noise  on  the  reconstructed 
image. 


Figure.  3(a)  shows  the  mean  RCS  after  diaphanization  of  the  B-52  scale  model  versus 
aspect  angle.  We  used  the  computer  to  generate  Gaussian  noise  with  two  levels  and  N2 
(see  Hg.3(a))  and  add  those  to  the  measured  fields  to  get  two  sets  of  noisy  data.  The 
corresponding  S/N  levels  of  and  N2  are  about  3  dB  and  -2  dB  respectively  relative  to  the 
mean  RCS.  The  frequency  response  of  the  first  aspect  angle  (broadside  view)  without  and 
with  noise  levels  Nx  and  N2  are  shown  in  Figs.  3(b),  (c),  and  (d)  respectively.  The  range 
profiles  corresponding  to  the  data  in  Figs.  3(b),  (c),  and  (d)  are  shown  in  Figs.  3(e),  (f),  and 
(g)  respectively.  If  we  stack  the  range  profiles  of  the  first  64  view  angles  together  in  a  per¬ 
spective  representation,  the  range  profiles  without  and  with  the  two  levels  of  Gaussian  noise 
will  appear  as  shown  in  Figs.  3(h),  (i),  and  (j).  The  images  reconstructed  from  these  three 
sets  of  range  profiles  are  shown  in  Figs.  3(k),  (m)  and  (n)  respectively.  The  results  presented 
in  Fig.  3  show  clearly  that  Gaussian  noise  affects  the  frequency  response  most,  the  range 
profile  less,  and  the  reconstructed  image  least  because  of  coherent  or  superposition  gain. 

V.  FUNDAMENTAL  CONCEPTS  FOR  DISTORTING  IMAGES 

From  the  above  examples,  it  can  be  seen  that  low  S/N  or  strong  Gaussian  noise  is  not  a 
sufficient  criterion  for  totally  obscuring  the  images.  In  order  to  develope  techniques  for  mak¬ 
ing  the  image  unintelligible  or  obscured,  one  must  resort  to  understanding  of  the  scattering 
mechanism  and  the  reconstruction  algorithm  involved. 

Radiation  can  originate  from  several  locations  on  a  scatterer.  These  include  excitation 
regions,  impedance  loads,  sharp  bends,  and  open  ends  [8].  The  finite  passband  data  and  finite 
range  resolution  make  the  range  profile  accentuate  the  discontinuities  in  reflectivity  function, 
and  the  image  can  be  interpreted  as  a  collection  of  the  discontinuities  of  reflectivity  function 
[5].  The  reconstruction  algorithm  followed  by  include  steps:  obtaining  filtered  range  profiles 
followed  by  back-projection.  Both  steps  as  pointed  out  earlier  are  coherent  summation 
processes.  Based  on  these  observations  we  will  propose  some  fundamental  concepts  or  rules 
for  obscuring  the  image. 

1.  Creation  of  artificial  discontinuities:  The  image  is  a  collection  of  the  discontinuities  of 

the  reflectivity  function.  For  a  set  of  given  objects,  their  images  can  be  obtained  and 


stored  to  act  as  the  data  base  for  subsequent  recognition  purpose.  If  unexpected  discon* 
tinuities  were  created,  the  reconstructed  images  of  the  target  will  be  different  from  their 
original  ones.  In  that  case,  visual  recognition  will  be  more  difficult  or  may  fail. 

2.  Creation  of  multiple  reflections:  Locations  of  peaks  in  the  range  profile  due  to  multiple 
reflections  correspond  to  the  differential  ray  path  between  the  multiple  scattering  points 
on  the  object  Such  peaks  correspond  to  as  physical  detail  on  the  object  and  help  there¬ 
fore  to  distort  and  obscure  information  about  object  If  the  object  is  shaped  so  that  multi¬ 
ple  reflections  are  deliberately  produced  to  be  the  dominant  contributions  to  the  scattered 
fields,  then  the  image  might  be  obscured  or  distorted. 

3.  Manipulation  of  local  reflectivity:  Here  we  can  make  the  reflectivity  at  a  given  point  a 
function  of  aspect  angle  and  frequency  or  randomly  modulate  it  in  time.  The  contribu¬ 
tions  to  a  pixel  from  different  filtered  range  profiles  obtained  at  different  aspect  angles 
are  added  coherently.  If  the  reflectivities  of  a  point  obtained  at  different  aspect  angles 
have  the  same  phase,  that  image  point  will  be  intensified  after  back-projection.  If  the 
phase  of  reflectivity  at  that  point  is  randomly  varied  for  each  frequency  and  aspect 
angle,  the  range  profile  might  produce  random  peaks  and  the  peaks  in  the  range  profile 
corresponding  to  those  reflection  points  might  become  smaller  because  their  reflectivities 
at  each  frequency  is  not  in  phase.  After  back-projection,  the  intensity  at  that  image  point 
might  become  smaller  and  the  image  might  be  contaminated  by  random  noise. 

Based  on  the  above  ideas,  we  will  describe  next  a  study  of  some  diaphanization  tech¬ 
niques  which  might  be  useful  in  reducing  the  RCS  or  obscure  the  images  and  make  them 
intelligible. 

VI.  EFFECT  OF  DIAPHANIZATION  TECHNIQUES  ON  RCS  AND  IMAGING 

Some  techniques  for  reducing  RCS  have  been  extensively  analyzed  in  [7],  which  include 
covering  the  body  with  absorbing  layers,  shaping  the  body,  and  impedance  loading.  In  the  fol¬ 
lowing,  we  will  use  microwave  diversity  imaging  as  a  new  tool  to  experimentally  and  numeri¬ 
cally  study  the  effect  of  these  techniques  on  RCS  and  imaging. 


6.1  Covering  With  Absorbers 

Microwave  absorbers  can  be  divided  into  four  groups  [9]:  magnetic  material  absorbers, 
resonant  absorbers,  broad-band  absorbers,  and  surface  current  absorbers.  To  reduce  the  RCS 
over  a  broad  spectral  bandwidth,  broad-band  absorbers  are  preferred.  However,  broad-band 
absorbers  available  today  are  not  suitable  for  aerodynamic  or  spacedynamic  applications 
because  they  usually  consist  of  light  foam  material.  Nevertheless,  broadband  absorbers  can  be 
a  useful  tool  in  experimental  study  the  effect  of  diaphanization  on  RCS  and  imaging  in  the 
context  of  the  work  reported  here. 

Shown  in  Fig.4(a)  and  (b)  are  the  arrangement  for  measuring  the  scattered  field  of  a 
square  conducting  plate  (40  x  40  cm)  and  the  mean  RCS  using  co-polarized  waves  before  and 
after  diaphanization  plotted  against  rotation  angle  $  with  zero  degrees  taken  as  the  broadside 
direction.  Both  faces  of  the  plate  were  covered  with  an  Emerson  and  Cumming  Eccosorb 
AN72;  the  spectral  sampling  range  of  the  received  fields  consists  as  before  of  201  frequency 
points  spanning  the  spectral  window  (6  to  17.5  GHz).  The  mean  RCS  after  diaphanization  is 
seen  to  have  been  reduced  for  low  rotation  angles,  and  the  amount  of  reduction  decreases  as 
the  rotation  angle  increases.  However,  the  mean  RCS  after  diaphanization  are  greater  than 
those  before  diaphanization  when  the  rotation  angle  <J>  is  greater  than  60°.  This  simple  experi¬ 
ment  illustrates  that  the  use  of  broad-band  absorber  does  not  guarantee  the  mean  RCS  reduc¬ 
tion  for  all  aspect  angle. 

The  reconstructed  projective  images  of  the  plate  before  and  after  diaphanization  are 
shown  in  Fig.4(c)  and  (d)  respectively.  The  dynamic  display  ranges  of  the  two  images  shown 
are  different  with  that  in  (d)  being  lower  than  in  (c).  From  the  above  example  together  with 
the  example  shown  in  Fig.2,  one  can  see  that  diaphanization  using  a  broad-band  absorber  to 
cover  hot  spots  is  generally  not  very  effective  in  obscuring  the  image. 

6.2  Target  Shaping 

As  one  of  the  methods  of  RCS  reduction,  the  aim  of  target  shaping  is  to  redirect  the 
incident  radiation  away  from  the  probing  system.  As  a  method  of  distorting  the  image,  target 
shaping  can  also  be  a  technique  of  generating  multiple  reflections  or  creating  unexpected 


discontinuities  or  reflections.  Multiple  reflections  will  distort  and  degrade  die  image  and  unex¬ 
pected  reflection  can  decoy  the  image. 

A  good  example  of  target  shaping  is  the  B-l  aircraft,  which  is  designed  to  have  little 
specular  reflection  over  practically  encountered  angular  window.  A  1:100  scale  model  B-l 
was  mounted  on  a  rotating  pedestal  with  inclination  angle  20°  as  shown  in  Fig.5(a).  The 
reconstructed  images  with  an  angular  window  extending  from  head-on  to  tail-on  (-90°,  90°)  of 
the  co-polarized  and  cross-polarized  echoes  are  shown  in  Fig.  5(b)  and  5(c)  respectively. 

An  important  feature  of  an  image  is  its  edges.  The  microwave  image  of  an  edge  is 
either  two  points  or  a  line  depending  on  the  angular  window  available  [5].  If  the  circumfer¬ 
ence  of  the  edge  is  covered  or  painted  with  resistive  materials  with  serrated  edges,  the  range 
profiles  will  be  distorted  and  will  contain  more  peaks  due  to  arising  reflections  from  added 
discontinuities.  It  is  expected  that  the  reconstructed  image  will  be  different  than  from  that 
obtained  without  artificial  discontinuities.  To  illustrate  this  a  conducting  plate 
(40  cm  x  40  cm )  is  covered  with  ceramic  substrates  as  shown  in  the  pictorial  view  of  Fig. 
6(a).  The  dielectric  constant  of  the  substrate  is  about  10  and  the  thickness  is  25  mils .  The 
plate  was  mounted  on  the  pedestal  with  an  inclination  angle  0  of  about  30°  as  shown  in 
Fig.6(b).  The  averaged  back-scattered  fields  of  the  plate  without  and  with  the  ceramic  sub¬ 
strate  covering  are  shown  in  Fig.  6(c).  The  reconstructed  images  of  the  plate  with  and  without 
the  substrate  using  angular  widows  (-30°,  60°)  and  (5°,  60°)  are  shown  in  Figs.  6(d)  to  (g) 
respectively.  The  bright  portion  in  the  center  part  of  the  plate  image  is  due  to  secondary 
reflections  from  the  support  holder  in  the  presence  of  the  plate  and  cannot  be  subtracted 
effectively  during  calibration.  Figure  6(d)  does  not  look  much  different  from  Ftg.6(e).  The 
negligible  difference  is  a  result  of  the  brightness  of  the  edges,  to  which  the  fields  of  those 
angles  close  to  0°  (broad-side)  contribute.  The  difference  between  Fig.  6(0  and  Fig.  6(g)  is 
more  discernible.  Both  of  these  images  were  reconstructed  from  an  angular  window  excluding 
those  views  close  to  the  broadside.  The  image  of  the  plate  has  been  decoyed  successfully  by 
the  presence  of  the  substrate  covering  for  certain  view  angles.  Visual  analysis  of  the  image 
might  accordingly  be  confused  and  produce  a  wrong  identification. 
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Uatetanaaely,  obacuring  the  image  in  this  fashion  causes  an  increase  of  the  RCS  due  to 

the  earn  reflections. 


f  3  impoflamea  Loading 

The  effect  of  linear  imp-danr*  loading  on  back-scattering  from  a  conducting  object  has 
been  cxaamtvdy  *******  An  N-port  loaded  acatterer  can  be  represented  by  a  general  network 
diagram  [101.  The  back-scattered  field  could  be  maxi  mi  Ted,  minimized  or  nullified  by  lumped 
i-pAwx  iniAn  The  required  losing,  however,  is  highly  dependent  on  frequency,  loading 
P nriiiM.  mmnk»r  of  uniting  points  and  aspect  angles  [2].  Therefore,  fixed  linear  impedance 
loading  is  not  an  effective  means  for  reducing  RCS  over  a  broad  spectral  range.  However,  to 
sady  the  effect  of  linear  loafing  on  the  image,  a  straight  wire  scatter er  with  linear  impedance 
is  mod  as  the  object,  and  the  moment  method  is  used  to  numerically  analyze  the 
desired  properties.  The  back-scattered  field  of  a  thin  wire  scatterer  with  and  without 
i— Um^mg  wen  numerically  calmlatnd  Shown  in  Fig.  7(a)  is  the  geometry  of  a  wire 
seamier  with  seven  symmetrically  spaced  impedance  loading  points.  The  locations  of  the 
w^tiag  points  are  shown  in  the  figure.  The  overall  length  is  30  on,  and  the  radius  to  length 
ratio  is  1/100.  The  spectral  window  used  is  from  6  to  16  GHz.  In  short,  the  length  of  the  wire 
in  terms  of  wavelength  ranges  from  6  to  16.  The  impedance  at  each  loading  point  is  confined 
to  30  Q.  The  p— above  me  arbitrarily  chosen.  The  complex  frequency  response  of 
the  wire  without  and  with  ;Tnpr«r— -  loading  me  shown  in  the  Fourier  slices  [3]  of  Fig.  7(b) 
and  7(e).  The  radial  distance  of  a  given  point  in  there  plots  represents  the  operating  fre¬ 
quency.  while  the  polar  angle  repwaena  the  angle  of  incidence  6  in  Fig.  7(a).  The  brightness 
of  each  point  is  proportional  to  the  amplitude  of  the  frequency  response.  The  reconstructed 
images  of  dm  wire  icanerer  without  and  with  loading  using  angular  windows  (20*  S  6  i  90°) 
ami  (20*  $  •  i  84*)  with  •  defined  in  Fig.  7(a)  are  shown  in  Figs.  7(c).  (d)  and  (f).  (g).  It  is 
unarming  to  note  that  a  ring  with  a  radius  equal  to  half  the  length  of  the  wire  appears  ut  the 
reconstructed  imr«  shown  in  Ftp.  7(c)  and  (d).  This  is  due  to  the  waves  traveling  from  one 
cad  poire  to  the  other  end  point  of  the  wtre(ll].  The  image  of  the  wire  with  impedance  load¬ 
ing  (see  fig.  7(d))  differs  notably  from  that  without  loading  (are  Fig  7(g)).  The  image  is  seen 


to  be  blurred.  The  loading  points  cause  discontinuities  that  vendiaie  strongly.  The  reflection 
coefficient  of  a  certain  loading  point  is  different  for  each  frequency  and  each  aqpect  angle, 
and  each  loading  point  will  cause  additional  traveling  waves  or  cause  additional  reflections. 
This  makes  the  image  of  a  loading  point  not  only  a  point  bat  instead  a  tweeping  curve. 


The  phenomenon  of  traveling  waves  along  the  wire  is  very  pronounced  because  the 
currents  are  confined  to  flow  hi  one  direction  only,  and  the  physical  optics  can  not  be  applied 
in  this  structure.  The  image  of  a  wire  scaoerer  is  therefore  not  expected  to  reflect  its  shape. 
For  a  general  2-d  or  3-d  object,  the  effect  of  traveling  waves  may  not  be  so  pronounced.  The 
effect  of  linear  loading  on  the  image  of  a  2-d  or  3-d  object  needs  more  experimental  and 
numerical  study. 


m  CONCLUSION 


We  have  employed  microwave  diversity  imaging  to  treat  the  problem  of  RCS  manage¬ 
ment  RCS  management  is  baaed  on  understanding  of  dm  scattering  merhaaitme  of  the  object 
and  the  reconstruction  algorithm  in  microwave  diversity  f**g  A  new  term  "diaphauiza- 
hon,"  defined  as  all  techniques  of  reducing  the  RCS  and  techniques  of  obscuring  the  image, 
has  been  introduced.  RCS  management  is  traaaad  not  only  from  the  detection  point  of  view 
but  also  from  the  imaging  perspective.  The  use  of  microwave  diversity  imaging  as  a  tool 
opens  a  new  direction  for  future  rematch  in  RCS  ' agetnent  stadias  We  have  also  thown 

that  microwave  diversity  imaging  is  quite  robust  to  the  Gaum tan  notes.  We  also  have  pro¬ 
posed  some  fundamental  concapts  for  obecuring  an  image.  Soma  tachniquae  of  reducing  the 
RCS  have  been  applied  to  experimentally  nady  their  effect  oa  RCS  and  imaging  It  is  found 
images  are  usually  obscured  at  the  price  of  increasing  the  RCS.  How  to  reduce  RCS  and 
obscure  images  umultaneousiy  is  under  investigation. 


RCS  have  been 
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FIGURE  CAPTION 


Fig  1.  Simplified  sketch  of  a  microwave  diversity  imaging  arrangement 

Hg.2  (a).  Sketch  of  the  B-52  test  object  and  the  angular  window  utiliaed.  (b).  Recon¬ 
structed  image  of  B-52  (c).  Sketch  of  the  paction  of  the  test  object  (darkened  lines) 
covered  with  broad-band  absorber,  (d).  The  mean  RCS  versus  aspect  angle  before  and 
after  diaphanization.  (e).  The  reconstructed  image  after  diaphanisarion  displayed  with 
the  same  dynamic  range  as  the  undiaphanieed  image  in  (b).  (f).  Reconstructed  image 
displayed  with  a  lower  dynamic  range. 


Fig.3  (a).  Averaged  scattered  fields  of  the  scale  model  B-S2  after  diaphani ration  versus 
aspect  angle  and  the  generated  Gaussian  noises  with  level  N  t  and  N*  Measured  fre¬ 
quency  response  of  die  first  view  angle  (broadside)  (b).  without  added  noise,  (c).  with 
noise  level  added,  and  (d)  with  noise  level  N2  added,  (e).  The  range  profile  of 
data  in  (b);  (f).  The  range  profile  of  data  in  (cX  (g).  the  range  profile  of  data  in  (d); 

A  perspective  stack  range  profiles  of  the  first  64  view  angles  (h)  without  adding  noise; 
(i)  with  noise  level  N{;  and  (j)  with  noise  level  N2.  The  reconstructed  image  (k) 
without  adding  noise;  (m)  with  noise  level  Af|;  and  (n)  with  noise  level  N2. 


Fig.4  (a).  Arrangement  for  the  measurement  of  a  square  conducting  plate,  (b).  Aver¬ 
age  scattered  fields  of  a  (40  x  40  cm)  conducting  plate  before  and  after  diaphanization. 
Reconstructed  images  of  the  plate  (c)  before  and  (d)  after  diaphani  ration. 


Fig.3  (a).  Pictorial  View  of  scale  model  of  B-l  bomber,  (b).  Reconstructed  image 
using  co-polarized  waves,  (c).  Reconstructed  image  using  cross-polarized  waves. 


Fig.6  (a).  A  conducting  plate  covered  with  ceramic  substrate,  (b).  Measurement 
configuration,  (c).  Mean  scattered  fields  of  the  plate  without  (solid  line)  and  with  sub-  I 

urates  (dash-dotted  line),  (d).  Reconstructed  image  of  the  conducting  plate  with  angu¬ 
lar  window  (-30s,  60s).  (e).  Reconstructed  image  of  the  ceramic -covered  plate  with 
angular  window  (-30s,  60s).  (f).  Reconstructed  image  of  the  conducting  plate  with 
angular  window  (5s,  60s).  (g).  Reconstructed  image  of  the  ceramic -covered  plate  with 
angular  window  (5s,  60s). 

Fig.7  (a).  Geometry  of  a  straight  wire  with  impedance  loading,  (b).  Real  pan  and 
imaginary  pan  of  the  frequency  response  of  the  straight  wire  without  impedance  load¬ 
ing  (c).  Reconstructed  image  of  the  wire  with  angular  window  (20s.  90s)  containing 
broadside  view.  (d).  Reconstructed  image  of  the  wire  with  angular  window  (20s.  84s). 
excluding  broadside  view  (e).  Reconstructed  image  of  the  loaded  wire  with  angular 
window  (20°,  90s).  (f).  Reconstructed  image  of  the  loaded  wire  with  angular  window 
(20s,  84s). 
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Fig.  2  (continued) 

(d) .  The  mean  RCS  venue  upect  angle  before  and  after  diaphaniration 

(e) .  The  reconstructed  image  after  diaphanixation  displayed  with  the 

time  dynamic  range  aa  the  undiaphaniaed  image  in  (b). 

(f) .  Reconstructed  image  displayed  with  a  lower  dynamic  range. 
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Fig.3  (a).  Averaged  scattered  fields  of  the  scale  model  B-52  after  diaphani- 
zation  versus  aspect  angle  and  the  generated  Gaussian  noises  with 
level  N\  and  AT 2' 

(to  be  continued) 
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Fig.  3  (continued) 

Measured  frequency  response  of  the  first  view  angle  (broadside) 
(b).  without  added  noise,  (c).  with  noise  level  N,  added, 
and  (d)  with  noise  level  N2  added. 

(to  be  continued) 
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Fig.  3  (continued) 

(e).  The  range  profile  of  data  in  (b);  (f).  The  range  profile  of  data  in  (c); 
(g).  the  range  profile  of  data  in  (d); 

(to  be  continued) 
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Fig.6  (a).  A  conducting  plate  covered  with  ceramic  substrate, 
(b).  Measurement  configuration. 

(to  be  continued) 
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Fig.  6  (continued) 

(c) .  Mean  scattered  fields  of  the  plate  without  (solid  line)  and  with 

substrates  (dash-dotted  line). 

(d) .  Reconstructed  image  of  the  conducting  plate  with  angular  window  (-30°,  60°). 

(e) .  Reconstructed  image  of  the  ceramic-covered  plate  with  angular  window  (-30°,  60°). 

(f) .  Reconstructed  image  of  the  conducting  plate  with  angular  window  (5°,  60°). 

(g) .  Reconstructed  image  of  the  ceramic-covered  plate  with  angular  window  (5°,  60°). 
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Kg.7  (a).  Geometry  of  a  straight  wire  with  impedance  loading, 
(to  be  continued) 
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Fig.  7  (continued) 

(b) .  Real  part  and  imaginary  part  of  the  frequency  response  of  the 

straight  wire  without  impedance  loading 

(c) .  Reconstructed  image  of  the  wire  with  angular  window  (20°,  90°) 

containing  broadside  view. 

(d) .  Reconstructed  image  of  the  wire  with  angular  window  (20°,  84°). 

excluding  broadside  view 

(e) .  Reconstructed  image  of  the  loaded  wire  with  angular  window  (20°,  90°). 

(f) .  Reconstructed  image  of  the  loaded  wire  with  angular  window  (20°,  84°). 
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A  aaaauraaant  aye tee  oapabla  of  acquiring  data  covering  tha  j 

aiorowmve/ml  1 1 Imeter-vave  (ww/aaw)  ration  of  tha  spectrum  la  a  valuable  tool 
in  aany  araaa  of  rsesaroh.  For  example,  in  tha  atudy  or  inverae  scattering, 
tha  aoattarad  fields  froa  a  teat  object  are  dependent  on  the  wavelength  of 
tha  illuaination,  therefore  the  availability  of  aoattered  data  over  a  wide 
bandwidth  would  be  laportant  in  gaining  Insight  into  the  scattering 
phenomenon  and  in  characterising  the  object.  In  the  areas  of  coherent 
laaglng,  the  wide  bandwidth  iaplieo  better  range  resolution  which  leads  to 
hifi  resolution  taagta  that  are  suited  for  recognition  by  the  human  eye-brain 
ays  tea  aa  in  aiorowave  and  allliaater  wave  diversity  imaging  where  angular, 
apeotral,  and  polarisation  degrees  of  freedom  are  combined  to  cost- 
effectively  form  images  of  scattering  objects  with  near  optical  resolution. 

It  la  well  known  that  combinations  of  different  polarization  states  of 
tha  incident  field  and  of  the  antennas  receiving  the  scattered  field  provide 
different  information  about  the  scatterer.  A  total  of  four  polarization 
state  combinations,  corresponding  to  the  four  elements  of  the  scattering 
matrix,  is  therefore  needed  to  fully  characterize  scattering  by  an  object. 


provided  that  other  factors  such  as  antenna  positions  remain  fixed.  Thus 
being  able  to  change  the  polarization  states  of  the  incident  field  and  of  the 
antennas  Measuring  the  scattered  field  is  important  in  the  study  of 
depolarisation  effects  and  target  identification  utilizing  polarisation 
descriptors.  Zt  is  also  worth  noting  that  polarisation  effects  are  sore 
pronounced  in  aultistatic  and  bistatlo  Measurements  than  aonostatlc 
Measureasnt. 

Besides  the  bandwidth  and  the  polarisation  oonalde rat Iona,  the  scattered 
fields  are  also  dependent  on  the  locations  and  the  directions  of  the 
transMittlng  and  the  receiving  (T/R)  antennas.  A  system  that  allows  change 
in  the  locations  and  the  orientations  of  the  T/R  antennas  is  useful  in 
aperture  synthesis  and  in  bistatlc  and  aultistatic  scattering  and  laaglng 
studies.  The  system  can  also  be  used  to  study  the  scattering  effects  in 
different  radiation  zones,  e.g.  near-field  versus  far  Tield  coherent  laaglng. 

In  summary,  a  systea  that  provides  wide  bandwidth,  variable  polarization 
states  of  the  incident  field  and  the  measured  field,  and  relocatable 
positions  of  the  transmitting  and  receiving  antennas  relative  to  the  object 
is  a  highly  desirable  tool  for  the  study  of  uw/mmw  scattering,  Imaging,  and 
target  representation.  The  aforementioned  features  are  also  referred  to  as 
wavelength,  polarization,  and  positional  or  angular  (aspect  related) 
diversities  respectively.  In  the  following  sections,  a  measurement  system 
designed  to  meet  these  desirable  features  is  described.  The  system  is  an 
upgraded  version  of  the  original  9-2018)  GHz  system  configuration  residing  at 
the  Electro-Optics  &  Microwave-Optics  Laboratory  of  the  Moore  School  of 
Electrical  Engineering  at  the  University  of  Pennsylvania  which  was  utilized 
in  our  earlier  microwave  diversity  imaging  studies.  The  upgraded  system  is 
intended  to  extend  operational  capabilities  to  60  GHz  and  to  enable 


aultlstatlo  acatterlng  esaaurwnta  employing  the  sum  operation  principles 
as  in  the  original  systaa  but  employing  different  gear. 


II.  ORIGINAL  SYSTEM  CONFIGURATION 

Shown  in  Fig.  1  is  a  schematic  diagram  of  the  present  system  operating 
from  2  to  18  GHz  (S,  C,  X,  and  Ku  bands  in  the  radar  letter-band 
nomenclature).  Two  microwave  sources,  a  sweeper  oscillator  (HP-8620C)  and  a 
sweep  synthesizer  (HP-8340A),  are  available;  either  one  can  be  used  as  the 
illuminating  source.  The  coherent  detector  is  a  network  analyzer  (HP-841  OB 
or  HP-8510A)  with  a  reflection/transmission  test  unit  (HP-8743  or  HP-8513). 
The  combination  of  HP-8620C/841  OB  can  operate  from  2  to  18  GHz,  while  HP- 
8340A/85120A  combination  allows  operation  from  45  MHz  to  26.5  GHz.  The 
illuminating  signal  is  amplified  by  a  solid-state  amplifier  (Litton  M5401-00) 
which  furnishes  1  watt  in  the  2-6  GHz  range,  and  by  a  TUT  amplifier  (Varian 
VZM6991K3)  with  maximum  10  watts  output  in  8-18  GHz  range  before  being  fed  to 
the  transmitting  antenna.  The  TWTA,  however,  is  usable  in  6-8  GHz  with  lower 
gain.  Thus  the  operating  bandwidth  is  limited  by  the  power  amplifier  to  the 
2-18  GHZ  range.  Switching  between  amplifiers  is  performed  manually. 

A  portion  of  the  transmitted  signal  is  extracted  for  use  as  the 
reference  signal  and  for  the  source  leveling  circuit  in  the  HP-8620C/8410B 
system.  For  the  HP-8340A/8510A  combination,  a  fraction  of  the  output  of  the 
source  before  amplification  is  used  directly  as  the  reference. 

A  parabolic  reflector  antenna  (AEL  APN  102C  18")  with  interchangeable 
right-hand-circular  and  left-hand-circular  polarization  (RHCP/LHCP)  feeds 
(Transco  9C27400)  is  used  in  the  2-1 8  GHz  to  produce  circularly  polarized 
illumination.  The  receiving  antenna  is  a  dual  polarization  horn  (EM  A6100) 
with  a  90  degree  hybrid  for  circular  polarization  measurement.  To  obtain 
measurements  of  horizontal  and  vertical  (H/V)  polarizations,  a  pair  of  dual 
polarization  horns  (EM  A6100)  is  used  without  the  hybrid.  The  circular 


polarisation  atatas  of  the  RHCP/LHCP  measurements  ara  switched  annually, 
while  thoaa  of  tha  H/V  polarisation  aaaauraaants  can  be  switched 
electronically. 

Tha  neoeived  fiald  la  aaplifiad  by  low  noise  solid-state  amplifiers,  an 
Avantek  MT-R035M  in  2-6  OHs  and  an  Avantek  AWT-18039M  in  8-18  GHz,  with  both 
having  about  90  db  gain.  Tha  pra-aapllf iars  ara  switchad  aanually.  Tha 
aaplifiad  signal  is  than  fad  to  tha  tast  channel  of  the 
reflection/ transmission  teat  sat  operating  in  transmission  node,  and  detected 
by  tha  network  analyser. 

Tha  teat  range  is  a  rectangular  anecholc  chamber  with  dimensions  50* 
long,  20'  wide,  and  20'  high.  Tha  tast  object  is  ao unted  on  a  pedestal 
looated  at  tha  oanter  of  tha  ohanber.  Tha  T/R  antennas  are  located  at  one 
and  of  tha  ohanber,  with  slight  angular  separation  between  them  for  insertion 
of  electromagnet ic  wave  absorbing  aaterials  to  minimize  mutual  coupling.  The 
main  beams  of  tha  T/R  antennas  ara  aanually  steered  toward  the  test  object 
for  batter  slgnal-to-nolse  ratio.  Aside  from  tha  slight  bistatlc  angle 
between  tha  T/R,  the  arrangement  is  essentially  monostatic.  The  pedestal  is 
rotated  by  a  stepping  motor  to  change  the  aspect  angles  of  the  test  object 
relative  to  the  T/R  which  are  stationary. 

III.  UPGRADED  SYSTEM  CONFIGURATION 

The  new  system  to  be  described  next  is  an  upgraded  version  of  the 
original  system  described  in  the  proceeding  section  with  several  new  salient 
features  added.  Many  of  the  components  used  in  the  original  system  are 
Incorporated  into  the  new  system.  Added  features  of  the  new  system  include: 

1.  Wide  bandwidth  coverage  ranging  from  2  GHz  to  60  GHz.  Thus  the 
operating  bandwidth  is  nearly  tripled,  covering  the  microwave  region  and  a 


lower  segment  of  the  millimeter-wave  region.  In  conventional  radar  letter- 
band  nomenclature,  the  system  coverge  is  from  S  band  up  to  U  band.  Small 
frequency  gaps  may  exist  between  bands  due  to  hardware  limitations. 

2.  Host  of  the  switching  between  components  of  different  bands  are 
controlled  electronically.  That  part  that  cannot  be  switched  remotely  is  due 
to  unavailable  hardware.  But  the  system  may  become  fully  automated  in  the 
future  as  broadband  switches  become  available. 

3*  An  x,y,z,0,4  positioner  that  enables  positioning  a  broadband 
receiver  module  to  nearly  any  coordinate  location  in  the  anechoic  chamber  and 
enables  orienting  the  receiver  module  towards  the  stationary  target  greatly 
enhances  the  versatility  of  the  new  system.  Bistatic  and  multistatic 
measurements  for  polarization  and  aperture  synthesis  studies  that  were 
impossible  or  difficult  earlier  can  be  performed  readily  with  the  new  system. 
A  sketch  of  the  anechoic  chamber  depicting  the  positioner  and  showing  the 
scattering  object  located  on  a  separate  elevation  over  azimuth  positioner  to 
change  the  target  aspect  relative  to  the  receiving  module  and  a  stationary 
transmitter  antenna  is  given  in  Fig.  2. 

Figure  3  depicts  the  basic  blocks  of  the  upgraded  system.  A  central 
consideration  in  the  design  of  the  upgraded  system  is  the  realization  of  a 
(2-60)  GHz  receiver  module  that  is  translatable  for  carrying  out  bistatlc  and 
multistatic  scattering  measurements.  The  source  generates  highly  accurate 
and  stable  RF  signals  from  2  to  60  GHz  in  either  sweep  mode  or  steppe* 
frequency  CW  mode.  The  transmitter  assembly  comprises  a  ban*  of  pewe- 
amplifiers  and  the  switches  needed  to  guide  the  signal  through  ar  ip- » ■  - 
signal  path  in  accordance  with  its  frequency  as  detailed 
the  amplified  signal  is  coupled  out  and  downconverted  *  .  •-  • 
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A— Anechoic  chamber 

T-Transmitting  antenna  (illuminator) 

O-Test  object  (scatterer)  mounted  chi  its 
own  election  over  azimuth  positioner 

R-Positioner  mounted  receiver  module 

P-Positioner  frame 

S-Window  for  cable  connections  to  external 
instrumentation 

Fig.  2.  Sketch  of  anechoic  chamber  showing  receiver  module  R  mounted  on 
x,y,z,0,$  positioner  and  scattering  object  0  mounted  on  an  elevation 
over  azimuth  positioner  in  front  of  a  stationary  illuminating 
antenna  T. 
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reference  signal  for  the  coherent  detector  (HP-851  OA  network  analyzer).  A 
local  oscillator,  phase-locked  to  a  10  MHz  or  5  MHz  external  reference 
derived  from  the  time  base  reference  oscillator  of  the  source,  furnishes  26.5 
GHz  for  the  down-conversion.  The  major  portion  of  the  amplified  signal 
radiates  through  the  transmitting  antenna  cluster,  which  consists  of  high 
gain  antennas  of  different  bands.  Further  details  of  the  transmitter 
assembly  are  given  in  Fig.  5* 

The  scattered  field  is  captured  by  the  receiving  antenna  cluster,  the 
front  end  of  the  receiver  assembly,  shown  schematically  in  Figs.  6  and  7 
which  consists  of  several  antennas  of  different  bands.  The  received  signal 
is  passed  through  a  bank  of  switches  to  effect  polarization  diversity 
reception  and  a  bank  of  low  noise  amplifiers  to  attain  sufficient  signal 
level  for  downconversion.  Again;  a  local  oscillator,  also  phse-locked  to  the 
10  MHz  or  5  MHz  signal  derived  from  the  time  base  reference  oscillator  of  the 
source,  provides  the  26.5  GHz  reference  to  downconvert  the  signal  to  below  18 
GHz.  The  downconverted  signal  is  then  passed  through  a  second  bank  of 
amplifiers  to  Increase  the  power  level  before  being  fed  to  the  test  channel 
of  the  coherent  detector  (HP  851  OA  network  analyzer)  via  a  low-loss  flexible 
coaxial  cable. 

The  ooherent  detector  is  in  essence  the  HP-851  OA  network  analyzer  with 
the  HP-8513  reflection/ransmission  test  set.  The  reference  channel  input, 
derived  from  the  transmitted  signal,  and  the  test  channel  input,  which  is  the 
downconverted  returned  signal,  are  again  downconverted  twice,  first  to  20  MHz 
in  the  HP-8513,  then  to  100  KHz  inside  the  HP-851 0A,  before  being  sampled  and 
fed  to  the  amplitude-phase  detector  portion  of  the  HO  851 0A.  The  network 
analyzer  is  always  phase-locked  to  the  reference  channel  input. 


There  are  two  position  control  interfaces:  one  is  for  the  receiver 
nodule  (receiver  assembly  With  the  receiving  antenna  cluster),  and  the  other 
for  the  test  object  aspect  control.  The  receiver  nodule  can  be  moved  to 
almost  any  location  in  the  anechoic  chamber,  and  be  oriented  so  that  the 
antenna  cluster  is  aiming  toward  the  test  object.  The  test  object  aspect 
control  is  realized  by  mounting  it  on  an  elevation  over  azimuth  positioner. 
All  these  units,  Including  the  source,  the  transmitter/receiver  assemblies, 
the  coherent  detector,  and  the  position  controls,  are  under  the  coordination 
of  a  computer/controller  shown  in  the  block  diagram  of  Fig.  3*  The 
computer/controller  controls  these  units  via  the  IEEE-488  bus  and  parallel 
digital  I/O  lines.  Next  we'll  describe  each  block  in  Fig.  3  in  more  detail. 

IV.  THE  SOURCE 

A  simplified  version  of  the  source  in  shown  in  Fig.  8.  It  is  comprised 
of  two  parts:  a  HP-8340A  sweeper  synthesizer  and  a  Watkins- Johnson  WJ- 1204-42 
frequency  extender,  (see  also  Figs.  4  and  5).  This  combination  delivers 
signals  with  greater  than  0  dBm  unleveled  power  and  synthesized  frequencies 
over  0.1-60  GHz  range.  The  HP-8340A  covers  frequency  range  10  MHz  to  26.5 
GHz  with  minimum  leveled  power  output  of  about  5  dBm  over  the  entire 
frequency  range.  Leveled  power  output  is  somewhat  higher  for  lower 
frequency.  For  the  WJ-1 204-42  frequency  extender,  0  dBm  input  power  is 
sufficient.  The  synthesizer  operates  in  either  CW  mode,  stepped  frequency 
mode,  or  sweep  mode.  The  synthesizer  is  controlled  by  the  IEEE-488  GPIB.  A 
highly  stable  Internal  time  base  at  10  MHz  of  the  synthesizer  is  also  used  as 
the  phase-locked  signal  for  the  local  oscillators  in  the  transmitter/receiver 


assemblies. 
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Simplified  schematic  of  the  source  in  the  transmitter 
assembly  of  Fig.  4. 


The  WJ-1 204-42  frequency  extender  is  in  essence  a  frequency  multiplier. 
It  passes  the  RF  0.1-26.5  GHz  signal  from  the  source  to  one  of  the  three 
output  ports  without  much  attenuation.  For  higher  frequencies,  the  RF  signal 
is  amplified,  then  switched  to  either  a  doubler  or  a  tripler,  with  the  former 
providing  frequencies  in  the  26.5-40  Qiz  range,  and  the  latter  in  the  40-60 
GHz  range.  The  minimum  unleveled  output  power  at  these  ports  is  -3  dBm  at 
24.0  GHz.  The  selection  of  any  of  the  three  output  ports  is  by  digital  I/O, 

3  bits  at  TTL  level. 

V.  TRANSMITTER  ASSEMBLY 

A  simplified  schematic  of  the  transmitter  assembly  is  shown  in  Fig.  9. 
The  transmitter  assembly  is  a  complex  of  power  amplifiers  and  switches.  It 
also  provides  the  reference  signal  for  the  coherent  detector  as  shown  earlier 
in  Figs.  4  and  5.  The  output  from  the  WJ-1 204-42  extender  is  fed  into  a  bank 
of  5  amplifiers  through  proper  switching.  After  the  amplification,  the  major 
part  of  the  signal  is  directed  to  the  transmitting  antenna  cluster  for 
illuminating  the  scattering  object,  while  a  fraction  of  it  is  extracted  to 
provide  the  reference  signal.  The  extracted  signal  is  down converted  to  below 
18  GHz  depending  on  band  by  a  local  reference  at  26.5  3iz,  which  is  generated 
by  a  local  oscillator,  phase-locked  to  the  time  base  of  the  source.  The 
downconverted  signal  is  to  be  used  as  the  reference  for  the  coherent 
detector. 

The  antenna  cluster  consists  of  a  broadband  parabolic  reflector  (2-40 

GHz,  3’  diameter).  Other  antennas  covering  the  40-60  GHz  band  are  yet  to  be 
» 

purchased.  There  are  two  spiral  feeds  for  the  parabolic  dish,  corresponding 

•Unanticipated  increases  in  the  cost  of  several  components  of  the  system 
prevented  the  purchase  of  these  antennas  and  a  26.5  GHz  phase-locked  local 
oscillator.  We  expect  these  items  to  be  purchased  from  future  grants  in 
order  to  extend  operation  to  60  GHz. 
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to  RHCP  and  LHCP  respectively.  The  switching  problem  between  the  output  of 


the  amplifier  bank  to  the  transmitting  antenna  cluster  prevents  at  present 


the  complete  automation  of  the  system.  The  recent  commercial  introduction  of 


an  electronically  controllable  switch  covering  the  2  to  40  GHz  band  by 


Tektronix  will  alleviate  this  problem.  Another  problem  is  the  switching  of 


RHCP/LHCP  spiral  feeds.  At  present  this  can  only  be  done  manually.  A 


Cassegrain  reflector  may  allow  mechanical  rotation  of  the  primary  feed3  if 


full  automation  is  desirable.  Even  though  semi -automated  the  present  (2-60) 


GHz  system  will  be  of  great  utility. 


VI.  RECEIVER  ASSEMBLY 


The  receiver  assembly  or  module  is  encased  in  a  box  (see  Figs.  6  and  7), 


and  will  be  attached  to  the  end  of  the  retractable  z-coordinate  boom  of  the 


x,y,z,e,<J>  positioner.  In  addition  to  the  translational  movement,  the 


receiver  module  can  be  steered  angularly  to  aim  the  antenna  main  beams 


towards  the  test  object.  Shown  in  Fig.  10  is  a  simplified  schematic  of  the 


receiver  assembly.  The  front  end  of  the  receiver  assembly  is  a  cluster  of 


seven  horn  antennas;  one  dual-polarizer  horn  for  the  2-18  GHz  band,  one  pair 


of  horns  of  opposite  circular  polarization  for  the  18-26.5  GHz  band,  and 


another  similar  pair  for  26.5-40  GHz,  and  finally  a  pair  of  standard  gain 


horns  for  40-60  GHz  where  at  present  only  linear  polarization  measurements 


can  be  conducted.  Immediately  following  the  antennas  are  polarization 


switches  to  conrol  the  polarization  component  of  the  scattered  field  to  be 


received.  The  received  signal  is  then  amplified  to  get  sufficient  power  for 


downconversion  to  within  18  GHz.  A  local  oscillator,  phase-locked  to  the 


reference  oscillator  source  time  base,  provides  the  26.5  GHz  reference  signal 
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Fig.  10.  Simplified  version  of  receiver  assembly  shown  in  Fig. 


for  the  mixers.  After  the  downconversion,  the  signal  is  passed  to  the  second 
stage  amplifier  bank  before  being  relayed  by  a  flexible  cable  to  the  coherent 
detector.  The  receiver  module  is  therefore  translatable  and  angularly 
steerable  to  enable  carrying  out  multistatic  measurements. 

VII.  COHERENT  DETECTOR 

The  coherent  detetor  is  the  HP-851  OA  network  analyzer  system  with  the 
HP-8513  reflect lon/transmiss ion  test  set.  For  our  purpose,  the  measurement 
.parameter  is  set  to  measure  S21*  with  reference  channel  input  derived  from 

the  transmitter  assembly  and  the  test  channel  input  derived  from  the  receiver 
assembly.  A  simplified  version  of  the  receiver  assembly  given  in  Figs.  6  and 
7  is  shown  in  Fig.  10  in  block  diagram.  A  simplified  block  diagram  of  the 
coherent  detector  (network  analyzer)  is  shown  in  Fig.  11.  The  test  set  acts 
as  a  first  stage  IF  converter,  downconverting  all  input  RF  signals  to  20  MHz, 
which  is  also  the  frequency  of  the  reference  time  base  within  the  HP-8510A. 
There  is  a  restriction  on  all  signal  power  at  the  first  harmonic  mixers.  Too 
much  power  will  damage  the  harmonic  mixers  permanently,  and  too  little  will 
cause  the  phase-locked  loop  to  lose  track.  The  acceptable  power  to  the 
mixers  should  not  exceed  -10  dBm;  damage  level  is  16  dBm.  Reference  power  to 
the  mixer  less  than  -50  dBm  will  not  be  adequate  for  phase-locking  of  the 
Internal  local  oscillator.  All  the  power  is  referred  to  the  mixer  inputs; 
there  is  additional  28  dB  loss  between  the  reference  input  and  the  mixer. 
Additional  attenuation  pads  may  have  to  be  inserted  to  ensure  good 
measurement  accuracy. 

The  first  IFs  at  20  MHz  enter  the  HP-851 0A  for  second  stage 
downconversion;  this  time  down  to  100  KHz.  The  second  IFs  are  then  sampled 
and  processed  to  be  placed  in  proper  format  for  feeding  to  the  digital 


Fig.  11.  Simplified  schematic  of  the  coherent  detector. 


phase/amplitude  detector  of  the  HP-851  OA.  The  data  are  read  out  via  the 
IEEE-1188  GPIB,  which  is  also  the  controlling  bus  for  the  HP-851  OA. 


For  our  application,  the  system  bus  connecting  the  HP-851  OA  and  the 
source  HP-8340A  must  be  localized  in  order  to  control  the  HP-8340A 
independently.  Due  to  this  localization,  some  utilities  encoded  in  the 
operating  system  of  the  HP-851 OA  cannot  be  used.  We'll  address  this  problem 
in  a  later  section. 

VIII.  x,y,z,9,»  POSITIONER 

The  positioner  is  designed  to  allow  multistatic  scattering 
measurement.  To  ameliorate  the  formidable  task  of  having  the  capability  of 
moving  the  receiver  assembly  and  the  transmitter  assembly  to  arbitrary 
locations,  we  opted  to  provide  mobility  for  the  receiver  assembly  alone.  The 
bulk  of  the  transmitter  assembly  itself  prohibits  such  mobility.  The 
transmitter  assembly  is  therefore  stationary  but  additional  degrees  of 
freedom  are  achieved  by  varying  the  aspect  angle  of  the  scattering  object. 
This  is  accomplished  by  mounting  the  test  object  on  a  computer  controlled 
elevation  over  azimuth  positioner.  Thus  the  position  and  direction  control 
is  achieved  by  the  combination  of  moving  and  steering  the  receiver  assembly, 
and  rotating  the  test  object.  A  sketch  of  the  x,y,z,9,$  positioner  within 
the  anechoic  chamber  showing  its  relation  to  other  system  components  was 
given  earlier  (see  Fig.  2).  Two  parallel  rails  run  along  the  length  of  the 
anechoic  chamber  (x  direction),  each  48'  in  length.  On  these  parallel  rails 
moves  the  transversal  rails  across  the  width  of  the  chamber  (y  direction), 
18'  in  length.  A  retractable  boom  (z  direction),  to  the  end  of  which  is 
attached  the  receiver  assembly,  rides  on  the  transversal  rail,  with  fully 
extended  length  of  18*.  There  are  two  additional  stepping  motors  used  to 


orient  the  elevational  and  the  azimuthal  angles  (0  and  $)  of  the  receiver 
assembly  such  that  the  antennas  can  be  made  to  point  toward  the  test  object. 
Thus  we  have  a  total  of  seven  stepping  motors  to  be  driven  independently. 

The  rails  are  chromium-plated  stainless  steel  shaft  supported 
intermittently  by  aluminun  rail  supporting  blocks.  The  rails  are  mounted  on 
I-beams  made  of  reinforced  fiber  glass,  chosen  to  minimize  radar  cross- 
section  and  minimize  reflections  within  the  chamber.  The  parallel  rails  are 
attached  to  the  wall  of  the  anechoic  chamber,  with  supporting  blocks  made  of 
fiber  glass.  The  transversal  rails  move  on  the  parallel  rails  by  ball¬ 
bearing  pillow-block  riding  on  stainless  shafts.  Similarly,  the  boom  is 
attached  to  a  platform  moving  on  the  transversal  rails  with  ball-bearing 
pillow-blocks  underneath.  The  z-coordinate  collapsable  boom  consists  of 
three  sections  of  fiber  glass  tubing  in  different  sizes,  each  section  is  6’ 
in  length.  The  three  telescoping  sections  contract  within  each  other  and 
expand  under  stepper  motor  control. 

To  move  the  positioner,  a  reference  coordinate  frame  and  its  origin  must 
be  defined  with  respect  to  the  anechoic  chamber.  The  coordinate  is  chosen  to 
coincide  with  the  chamber  frame,  and  the  origin  to  be  the  center  of  the 
chamber  or  the  center  of  rotation  of  the  scattering  object  0  mounted  on  the 
elevation  over  azimuth  positioner  shown  in  Fig.  2.  Whenever  the  system  is 
powered  up,  or  before  it  is  being  turned  off,  the  x,y,z,e,$  positioner  must 
be  relocated  to  that  origin.  Therefore  there  must  be  a  power-up  procedure  or 
power-off  procedure  to  initiate  the  positioner.  After  that,  the  movement  of 
the  positioner  and  the  steering  of  the  receiving  antenna  are  referred  to  the 
origin.  Stepper  motor  pulse  counting  or  mechanical  or  optical  sensors  would 
be  incorporated  in  the  positioner  to  feedback  a  signal  to  monitor  the 
position  and  orientation  of  the  receiver  module.  For  example,  by  centering 


the  receiver  module  at  a  point  representing  the  center  of  rotation  of  the 
scattering  object  mounted  on  the  elevation  over  azimuth  positioner  an  origin 
for  a  spherical  coordinate  system  can  be  defined.  The  position  vector  of  the 
receiver  module  after  displacement  can  then  be  determined  by  computing  the 
relative  x,yf  and  z  displacements  by  counting  the  number  of  calibrated  steps 
of  the  x,y,z  stepper  motors  executed  to  get  to  that  point.  The  values  of  x,y 
and  z  determined  in  this  fashion  enable  calculating  the  radial  coordinate  r  - 
2  2  2  1  /2 

(x  +y  +z  )  and  the  angular  spherical  coordinates  6'  and  ♦  *  which  enable 
defining  the  orientation  0  and  $  of  the  receiver  module  required  to  point  it 
towards  the  object. 

IX.  SOME  DESIGN  CONSIDERATIONS 

Several  aspects  of  the  measurement  system  design  are  currently  being 
addressed.  The  disassociation  of  the  source  HP-8340A  from  the  control  of  the 
HP-851  OA  network  analyzer  is  necessitated  by  the  fact  that  the  new  system 
must  accommodate  measurement  at  frequencies  higher  than  the  capabilities  of 
the  HP-83^0A/8510A.  In  normal  applications,  the  HP-851  OA,  running  under  the 
software  provided  by  Hewlett-Packard  Co.,  acts  as  the  controller  of  the  HP- 
8340A.  The  HP-851 OA  is  the  only  instrument  to  be  controlled  by  the  external 
controller.  In  the  new  system,  the  frequency  control  of  the  HP-SS^OA  is 
relegated  from  the  HP-851  OA  internal  controller  to  the  external  controller. 
Except  for  the  frequency  stimulus  functions,  the  remaining  utilities  offered 
by  the  HP-851 OA  can  still  be  used.  The  same  situation  is  encountered  in  the 
HP  and  Hughes  millimeter-wave  network  analyzer  systems,  since  these  systems 
also  use  HP-8510  as  the  key  building  block.  The  operations  involving  the 
frequency  stimulus  must  be  programmed  by  the  user.  More  information  may  be 
needed  from  Hewlett-Packard  Co.  to  write  such  a  program. 


In  designing  a  network  analyzer,  efforts  are  made  to  maintain  the 
attenuation  and  phase  balance  of  the  test  channel  and  the  reference  channel. 
Consider  our  system  as  a  network  analyzer  operating  in  transmission  mode  with 
the  transmitting  section  as  the  reference  channel  and  the  receiving  section 
as  the  test  channel.  It  is  ideal  that  the  path  length  and  the  attenuation  of 
these  two  channels  are  identical.  The  inequality  of  path  length  introduces 
linear  phase  difference.  Different  attenuation  causes  wrong  magnitude  ratio 
reading.  Should  any  nonlinearity  exist,  e.g.,  phase  instability  or  excess 
power  in  any  channel,  the  measurement  errors  will  be  noticeable.  We  expect 
programmable  power  level  feature  of  the  HP  83^  synthesized  sweeper  used 
will  enable  precalibration  of  power  levels  through  the  system  at  different 
frequencies  to  alleviate  or  eliminate  any  nonlinearity  errors. 

Every  component  has  a  power  handling  capability  such  that,  within  a 
certain  power  range,  the  component  is  operating  in  its  linear  region. 
Therefore,  it  is  essential  to  ensure  such  linear  operation  for  each 
component;  not  only  to  remove  nonlinear  distortion,  but  also  to  use  the 
component  at  its  best  performance.  To  satisfy  this  consideration, 
attenuators  are  Inserted  at  various  stages  inside  the  assemblies.  The  price 
paid  for  achieving  this  is  that  some  power  is  wasted  and  the  lengths  of 
signal  paths  becomes  unequal.  But  the  gain  in  performance  sanctions  such 
trade-in.  The  phase  difference  caused  by  Inequality  of  signal  paths  can  be 
eliminated  by  a  calibration  procedure  which  measures  the  total  phase 
difference  of  a  through  connection  and  stores  it  as  a  calibration  set,  l.e., 
by  measuring  and  the  frequency  response  of  the  measurement  system  Itself  as 
was  done  in  the  original  measurement  system.  The  measurements  made  afterward 
are  corrected  against  the  calibration  data  set  to  compensate  the  phase 


difference.  The  same  calibration  procedure  can,  to  a  certain  degree,  correct 
the  residual  nonlinearity  and  power  level  fluctuations  within  the  system. 

Although  the  design  of  the  system  has  been  for  continuous  frequency 
coverage  from  (2-60)  GHz,  a  few  small  frequency  gaps  may  exist  within  the 
overall  bandwidth  because  of  hardware  limitations.  The  missing  of  these 
frequencies  may  not  have  significant  effects  on  various  applications. 

The  second  stage  amplification  and  the  frequency  downconversion  in  the 
receiver  assembly  are  indispensable  due  to  the  considerable  length  of  the 
flexible  coaxial  cable  connecting  the  receiver  assembly  to  the  coherent 
detector.  Flexibility  is  needed  since  the  receiver  assembly  must  be 
translatable  for  multistatic  operation.  To  enable  multi  static  measurements 
by  positioning  the  receiver  module  within  the  front  half  of  the  chamber  a  35' 
coaxial  cable  is  required.  A  low-loss  phase  stable  coaxial  cable  must  be 
used.  An  ideal  attractive  future  alternative  is  to  use  a  microwave  optical 
fiber  link  when  the  bandwidth  of  such  links  is  made  to  approach  Id  GHz. 

The  actual  data  acquisition  locations  produced  by  the  x,y,z,0,$ 
positioner  are  critical  in  multistatic  measurements,  especially  when 
frequencies  in  the  millimeter  wave  range  are  utilized.  In  terms  of  imaging, 
knowledge  of  the  exact  positions  where  the  data  are  taken  is  essential  to 
guarantee  good  image  reconstruction.  Because  of  the  considerable  weight  of 
the  transversal  rails  and  the  boom,  slight  bending  of  the  rails  is 
inevitable.  We  estimate  that  there  will  be  0.6"  bending  in  the  z  direction 
at  most.  It  would  be  better  if  this  bending  can  be  compensated,  say  by 
retracting  or  extending  the  boom.  It  may  become  necessary  to  calibrate  the 
positioner  for  all  locations  in  the  chamber  if  the  position  deviations  turn 
out  to  be  appreciably  detrimental.  Another  aspect  worth  considering  is  the 
wobbling  of  the  boom  as  it  changes  positions.  Adequate  time  must  be  allowed 


until  the  boom  stabilizes  before  taking  measurement.  This  would  be  taken 
Into  account  in  writing  data  acquisition  programs. 

Despite  the  fact  that  materials  chosen  are  intended  to  minimize  the 
reflection  from  the  positioner  and  the  receiver  assembly,  the  effect  from 
these  reflections  can  be  further  minimized  by  two  ways.  The  first  method  is 
to  cover  reflecting  parts  of  the  positioner,  Including  the  rails,  the  boom, 
and  the  receiver  assembly  box  with  appropriate  broadband  absorber  material. 
The  secod  method  is  to  measure  and  store  the  chamber  response  for  every 
positioner  location  on  an  intended  scanning  locus,  and  subtract  this  response 
from  later  measurements.  This  method  may  be  impractical  in  that:  firstly, 
the  same  measurement  sequence  has  to  be  done  twice,  one  for  room  clutter, 
then  for  test  object;  secondly,  the  accuracy,  or  rather  the  inaccuracy,  of 
the  stepping  motors  may  prevent  the  positioning  repeatability  of  the  two 
successive  measurements.  The  latter  may  nullify  the  very  intention  for 
clutter  cancellation.  A  combination  of  both  approaches  may  offer  the  best 
optimal  solution. 

X.  DATA  ACQUISITION  PROGFWtgNG 

This  section  outlines  the  programming  for  data  acquisition.  As  we 
discussed  earlier,  the  x,y,z,8,4  positioner  would  place  Itself  at  the  origin 
of  the  chamber  coordinate,  the  center  of  the  chamber,  before  and  after  data 
acquisition.  This  may  be  done  by  a  start-up  program  whose  sole  finction  is 
to  move  the  positioner  until  it  finds  the  chamber  origin.  This  routine 
should  be  invoked  whenever  the  system  is  started  and  before  it  is  turned  off. 

We  concentrate  now  on  the  programming  for  data  acquisition,  assuming 
that  the  start-up  routine  has  been  executed  and  the  positioner  is  at  the 
origin.  Shown  in  Fig.  12  is  the  program  block  diagram  of  a  typical  data 


acquisition  session.  Also  shown  below  each  block  are  the  subroutines  needed 
to  perform  the  function  of  that  block;  The  principle  of  modularity  in 
programming  is  closely  observed  so  that  each  user  can  tailor  his  own  program 
easily  using  these  routines  as  building  blocks.  A  description  of  these 
subroutines  is  given  next: 

SOURCE 

The  source  subroutine  performs  dialogue  to  accept  frequency  definitions, 
e.g.,  start  and  stop  frequencies,  number  of  frequency  samples,  power 
level  at  each  frequency,  etc.,  from  the  user.  From  these  user-specified 
parameters,  the  routine  initializes  the  source  and  the  coherent 
detector,  and  generates  the  switching  map  for  each  of  the  cross-band 
switches  if  the  frequency  range  involves  cross-band  measurement.  Also 
the  desired  polarization  combination  is  selected. 

FREQ 

The  frequency  subroutine  uses  the  parameters  generated  by  SOURCE  to 
perform  frequency  setting  of  a  single  sweep  of  the  specified  frequency 
range.  It  not  only  controls  the  frequency  setting  and  timing,  it  is 
also  responsible  for  setting  up  the  proper  signal  path  for  cross-band 
measurement  using  the  switching  map  generated  by  SOURCE. 

MOVE 

The  move  subroutine  controls  the  relative  movement  of  the  x,y,z,0,4> 
positioner  and  the  orientation  0,<>  of  the  receiver  module  through  the 
control  of  different  stepping  motors.  The  inputs  to  this  routine  are 
the  distances  in  x,y,z  directions  and  the  elevation  (0)  and  azimuthal 


($)  angles  of  the  receiver  module.  All  parameters  are  relative,  not 
related  to  the  absolute  coordinates  of  the  chamber.  Therefore,  these 
parameters  are  pre-calculated. 


DATA 


The  data  subroutine  reads  the  measurement  of  one  stepped  frequency  sweep 
from  the  network  analyzer  and  delivers  the  data  to  the  auxiliary  storage 
devices.  The  data  format  must  be  conformed  to  the  other  data  processing 
programs. 


LOCUS 


The  locus  subroutine  calculates  the  scanning  locus  of  a  data  acquisition 
session.  It  is  highly  dependent  on  the  specific  locus  the  user  intends 
to  synthesize.  It  is  difficult  to  write  a  subroutine  that  gives  options 
for  all  possible  loci,  although  some  regular  loci  may  be  implemented. 
The  key  parameters  to  be  generated  in  this  routine  are  the  relative 
changes  in  positions  x,y,z  and  angles  8,<fr  given  a  locus  definition.  The 
locus  definition  includes  the  starting  position  for  scanning,  the  number 
of  positions,  the  locus  equations,  etc. 


ROTATE 


The  rotate  subroutine  controls  the  rotation  of  the  test  object  mounted 
on  the  pedestal  in  order  to  change  its  aspect  relative  to  the  stationary 
illuminating  antenna  and  the  mobile  receiver  module.  This  enables 
complete  multistatic  measurement  to  be  carried  out. 


-31- 


& 


I 


3 


4 


The  aforementioned  subroutines  are  the  basic  building  blocks  for  any 
data  acquisition  program.  We  describe  next  the  steps  in  a  typical  data 
acquisition  session,  illustrated  in  Fig.  12. 

The  first  step  is  the  source  initialization,  the  SOURCE  subroutine  is 
the  key  routine  to  be  called.  The  functions  to  be  performed  at  this  stage 
are: 

*  F!requency  sweep  format  definition:  start  and  stop  frequencies,  sweep 
mode,  sampling  number. 

*  HP-83J10A  initialization:  power  level,  sweeping  time,  averaging 
factor. 

*  HP-851  OA  initialization:  measurement  parameter  S  21  *  system 

localization,  display  formats. 

*  Switching  map  calculation  for  setting  the  signal  path. 

*  Selection  of  polarization  states. 

After  system  initialization,  a  calibration  measurement  is  performed  for 
later  data  correction.  Two  subroutines  are  used  in  this  procedure: 

*  Move  the  receiver  assembly  to  the  calibration  location  by  the  use  of 
MOVE  routine.  A  good  choice  of  the  calibration  location  is  the 
location  closest  to  the  transmitter  assembly  as  if  in  monostatic 
arrangement.  Also  the  receiver  assembly  would  be  pointing  towards 
the  origin  of  the  coordinate  (e.g.,  center  of  rotation  of  object  on 
elevation  over  azimuth  positioner). 

*  A  reference  target,  a  metallic  sphere  or  a  metallic  cylinder,  is 
centered  at  the  origin  of  the  coordinates. 

*  A  single  frequency  sweep  of  data  is  taken  and  stored  in  the  HP-851  OA 
or  the  computer  for  later  correction.  Subroutines  used  are  FREQ  and 
DATA. 
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After  the  calibration  data  measurement,  the  user  enters  the  locus 
definition  by  the  use  of  LOCUS  subroutine  as  follows: 

*  Scanning  locus  definition:  starting  coordinate,  locus  equation,  and 
number  of  scanning  positions. 

*  Scanning  position  calculation.  Coordinate  increments  of  successive 
scanning  locations  are  calculated. 

*  Scanning  angle  calculation.  Angular  increments  needed  to  steer  the 
receiver  module  toward  the  object  at  successive  scanning  locations 
are  calculated. 

*  Define  the  test  object  orientation  for  each  scanning  location.  This 
step  is  needed  only  when  more  complex  loci  is  desired. 

After  calculating  all  the  parameters  needed,  now  the  positioner  is  moved 
to  the  initial  scanning  location: 

*  Move  the  positioner  to  the  starting  scanning  location  by  the  use  of 
MOVE. 

*  Rotate  the  test  object  to  the  initial  orientation  by  the  use  of 
ROTATE. 

The  system  now  is  ready  to  collect  the  data.  The  data  acquisition 
consists  of  the  following  steps: 

*  Collect  one  frequency  sweep  opf  data  by  the  use  of  FREQ. 

*  Get  the  data  from  the  system  and  store  them  by  the  use  of  DATA. 

*  Move  to  the  next  scanning  location  and  orientation  of  the  receiver 
module  by  the  use  of  MOVE. 

*  Change  the  object  orientation  as  desired  by  the  use  of  ROTATE. 

*  Repeat  the  same  procedure. 

After  the  end  of  data  acquisiton,  the  positioner  and  the  object  are 
reset  to  their  initial  location  and  orientations. 


Note  that  in  Fig.  12,  there  is  an  intermediate  step  of  clutter 
measurement.  This  step  is  essentially  the  same  as  data  acquisition,  except 
that  now  there  is  no  object  placed  on  the  pedestal.  The  measured  quantity  is 
the  aspect  dependent  room  clutter.  After  clutter  measurement,  the  positioner 
is  moved  to  the  initial  scanning  location. 

XI.  CONCLUSIONS 

The  new  upgraded  system  possessing  wavelength,  polarization,  and 
position  diversity  measurement  capability  in  the  (2-60)  GHz  range  is 
described.  The  equipment  and  components  are  delineated.  Many  aspects  of  the 
design  problems  and  the  viable  solutions  to  these  problems  are  presented. 
The  programming  considerations  are  addressed,  with  a  typical  data  acquisition 
session  as  an  example.  Most  of  the  components  are  in  hand  and  in  a  state  of 
partial  assembly.  The  construction  and  assembly  of  the  x,y,z,0,<j>  positioner 
is  about  SOt  completed.  The  implementation  of  the  system  is  a  challenging 
task.  However  once  fully  operational,  the  system  will  be  a  unique  tool  in 
the  study  of  multistatic  scattering,  polarization  effects,  and  radar  imaging 
modalities  in  the  (2-60)  GHz  range.  The  system  will  also  be  invaluable  in 
the  generation  of  high  fidelity  target  representations  and  signatures  in  the 
study  and  development  of  automated  recognition  schemes  based  on  models  of 
neural  networks. 

The  components  of  the  system  have  been  chosen  with  multifunctionality  in 
mind.  For  example  the  components  of  the  frequency  agile  phase-locked 
reference  oscillator  in  Fig.  4  can  be  reconfigured  as  shown  in  Fig.  13  to 
furnish  a  high-speed  synthesized  source  for  illuminating  a  target  for  rapid 
data  acquisition  that  is  one  to  two  order  of  magnitude  faster  than  with  the 
hp  8340A  based  system. 
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Fig.  13.  Fast  switching  (.7-40)  GHs  synthesized  transmitter 
module. 


APPENDIX  IV 


OPTICAL  ANALOGS  OF  TWO-DIMENSIONAL  NEURAL 
NETWORKS  AND  THEIR  APPLICATION  IN  RECOGNITION  OF  RADAR  TARGETS 

N.H.  Far hat,  S.  Mlyahara  and  K.S.  Lea 
University  of  Pennsylvania 
The'  Electro-Optics  and  Microwave-Optics  Laboratory 
200  S.  33rd  Street 
Philadelphia,  PA  19104-6390 

ABSTRACT 

Optical  analogs  of  2-D  distribution  of  Idealized  neurons  (2-D 
neural  net)  baaed  on  partitioning  of  the  resulting  4-D  connectivity 
matrix  are  discussed.  These  are  desirable  because  of  compatibility 
with  2-D  feature  spaces  and  ability  to  realize  denser  networks.  An 
exsnple  of  their  use  with  sinogram  classifiers  derived  from  realis¬ 
tic  radar  data  of  scale  models  of  three  aerospace  objects  taken  as  learn¬ 
ing  set  Is  given.  Super-resolved  recognition  from  partial  informa¬ 
tion  that  can  be  as  low  as  20X  of  the  sinogram  data  Is  demonstrated 
together  with  a  capacity  for  error  correction  and  generalization. 

INTRODUCTION 

Neural  net  models  and  their  analogs  furnish  a  new  approach  to 
signal  processing  that  Is  collective,  robust,  and  fault  tolerant. 

Optical  Implementations  of  neural  nets2*2  are  attractive  because  of 
the  Inherent  parallelism  and  massive  Interconnection  capabilities 
provided  by  optics  and  because  of  emergent  optical  technologies  that 
promise  high  resolution  and  high  speed  programmable  spatial  light 
modulators  (SIMs)  and  arrays  of  optical  bistability  devices  (optical 
decision  making  elements)  that  can  facilitate  the  Implementation  and 
study  of  large  networka.  Optical  implementation  of  a  one-dimension¬ 
al  network  of  32  neurons  exhibiting  robust  content -addressability 
and  associative  recall  haa  already  been  demonstrated  to  illustrate 
the  above  advantages.2  Extension  to  two-dimensional  arrangements 
are  of  Interest  because  these  are  suitable  for  processing  of  2-D 
Image  data  or  Image  classifiers  directly  and  offer  a  way  for  optical 
implementation  of  large  networka . 4 

In  this  paper  we  will  discuss  content  addressable  memory  (CAM) 
architectures  based  on  partitioning  of  the  four  dimensional 

memory  or  Interconnection  matrix  encountered  in  the  storage  of  2-D 
entitles.  A  specific  architecture  and  implementation  based  on  the 
use  of  partitioned  unipolar  binary  (u.b.)  memory  matrix  and  the  use 
of  adaptive  thresholding  in  the  feedback  loop  are  described.  The 
use  of  u.b.  memory  masks  greatly  simplifies  optical  Implementations 

3  4 

and  facilitates  the  realization  of  larger  networks  '(10  -10  neurons)  . 
Numerical  simulations  showing  the  use  of  such  2-D  networks  in  the 
recognition  of  dilute  point-like  objects  that  arise  in  radar  and 
other  similar  remote  sensing  Imaging  applications  are  described. 

Dilute  objects  pose  e  problem  for  CAM  storage  because  of  the  small 


Hawing  distance  betv««n  than.  Her*  w«  show  chac  coding  in  Che  form 
of  a  binogAM  ctaSi-i^izA.  of  eha  dilute  object  can  remove  this  limita¬ 
tion  permitting  recognition  from  partial  versiona  of  the  stored 
entities.  The  advantage  of  this  capability  in  super-resolved  recog¬ 
nition  of  radar  targets  is  discussed  in  the  context  of  a  new  type  of 
radar  diversity  imaging,  studied  extensively  in  our  laboratory,  that 
is  capable  of  providing  sinogram  information  compatible  with  2-D  CAM 
storage  and  Interrogation.  Super-resolved  automated  recognition  of 
scale  noddle  of  three  aero-space  objects  from  partial  information  as 
la#  as  20Z  of  a  learned  entity  is  shown  employing  hetero-associative 
storage  where  the  outcome  is  a  word  label  describing  the  recog¬ 
nised  object.  Capacity  for  error  correction  and  generalisation  were 
also  observed. 


TWO-DIMENSIONAL  NEURAL  NETS 


content  addressable  or 
i  set  of  M  2-D  bipolar 


Storage  and  readout  of  2-D  entities  in  a 
associative  memory  is  described  next.  Given 

.  fm\ 

binary  pattens  or  entitles  v^  m»l,2...M  each  of  HxN  elements  rep¬ 
resented  by  a  matrix  of  rank  N,  these  can  be  stored  in. a  maimer  that 
is  a  direct  extension  of  the  1-D  case  as  follows:  For  each  element  of 
a  matrix  a  new  NXN  matrix  is  formed  by  multiplying  the  value  of  the 
element  by  all  elements  of  the  matrix  including  itself  taking  the 
self  product  as  zero.  The  outcome  is  a  new  set  of  N*  binary  bipolar 
matrices  each  of  rank  N.  A  formal  description  of  this  operation  is. 
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(1) 


which  is  a  four  dimensional  matrix.  An  overall  or  composite  synaptic 
or  connectivity  memory  matrix  is  formed  then  by  adding  all  4-D 
,(m) 


matrices  T 
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This  symmetric  4-D  matrix  has  elements  that  vary  in  value  between  -M 
to  M  also  in  steps  of  two  as  for  the  1-D  neural  net  case  and  which 
assume  values  of  *4*1  and  -1  (and  zeros  for  the  self  product  elements) 
when  the  matrix  is  clipped  or  binarized  as  is  usually  preferable  for 
optical  Implementations.  Two  dimensional  unipolar  binary  entities 


b^j  are  frequently  of  practical  Importance, 
in  the  usual  way  into  bipolar  binary  matrices  through  v 


These  can  be  transformed 
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which  are  then  used  to  fora  the  4-D  connectivity  matrix  or  memory  as 
described.  Also,  as  in  the  1-D  neural  net  case,  the  prompting  entity 

k(m> 


can  be  unipolar  binary  b 


ij 


which  would  simplify  further  optical 


implementations  in  incoherent  light. 

Architectures  for  optical  implementation  of  2-D  natural  nats 
muat  contand  with  tha  task  of  realizing  a  4-D  memory  matrix.  Hara 
a  schama  is  prasantad  that  is  basad  on  tha  partitioning  of  tha  4-D 
memory  matrix  into  an  array  of  2-D  matrices  of  rank  N. 

Nearest  neighbor  search  of  tha  memory  matrix,  for  a  given  entity 

(no) 

is  dona  by  forming  tha  estimate, 

Sy*0  »  ^  Tljk£  b^Jo)  ....  i.l.k.l  -  1,2,... N  (3) 


followed  by  thresholding  to  obtain  a  new  u.b.  matrix  which  is  used 
fao) 

to  replace  "iT  in  eq.  (3)  and  tha  procedure  is  repeated  until  tha 
resulting  matrix  converges  to  tha  stored  entity  closest  to  tha  ini- 
tinting  matrix  b^  .  Tha  operation  in  aq.  (3)  can  be  interpreted 
as  first  partitioning  of  tha  4-D  matrix  into  an  array  of  2-D 

subantrlces  of  rank  N:  T11V>,  T1?vf . T^, 

T2HkiJ  hva*  TN2ki*  *  *  *  ’  TNNki  M  *«P*c««*  schematically  in  Pig. 

1(a)  where  tha  partition  submatrlcas  are  arranged  in  a  2-D  array. 

This  first  step  is  followed  by  multiplication  of  by  each  of  tha 

partition  submatrlcas,  on  an  element  by  element  basis,  and  summing 
tha  products  for  each  submatrix  to  obtain  tha  first  estimate  b^0^ . 

Tha  tensor  multiplications  and  summation  operations  called  for  in 
aq.  (3)  are  carried  out  in  Fig.  1(a)  by  placing  a  spatially  inte¬ 
grating  photodetector  (PD)  behind  each  submatrix  of  the  partitioned 


Pig.  1.  Optical  analog  of  2-D  neural  net.  (a)  Architecture  based 
on  partitioning  of  connectivity  matrix,  (b)  Opto-elec- 
tronlc  embodiment. 


memory  aask  which,  is  assumed  for  the  tlaa  boiag  to  bo  realized  by 
pixel  transmittance  aodulotlon  in  on  ldool  tranoporoncy  capablo  of 

assuming  nogaclva  eranaalctanca  valuaa.  Tha  input  entity  b£“^  la 
aaauaod  to  bo  displayed  on  a  suitable  LED  array.  The  LED  diaplay  of 
la  multiplied  by  the  Ideal  tranaalttance  of  each  of  the  parti¬ 
tion  aubaatrlcae  by  iasglng  the  diaplay  on  each  of  these  with  exact 
registration  of  pixels  by  aeons  of  a  lenslet  array  as  depicted  in 
Fig.  1(b).  The  output  of  each  PD,  proportional  to  one  of  the  com¬ 
ponents  of  eq.  3,  Is  thresholded ,  sap Ilf led,  and  fed  back  to  drive 
an  associated  LED.  The  (1, j)-th  LED  Is  paired  with  the  (i,j)-th  PD. 
This  completes  the  interconnection  of  the  2-D  array  of  NxN  neurons 
In  the  above  architecture  where  each  neuron  coansmlcates  its  state 
to  all  other  neurons  through  a  prescribed  four  dimensional  synaptic 
or  asaory  aatrlx  in  which  information  about  M  2-D  binary  aatrlces  of 
rank  N  (entitles)  have  been  stored  distributive!? .  The  number  of 
2-D  entitles  that  can  be  stored  In  this  fashion  iiM*  n2/81oN,  which 
follows  directly  from  the  storage  capacity  formula  for  the  1-D  neural 
net  case  by  replacing  >  by  M*. 

The  added  complexity  aaaoclated  with  having  to  realise  a  bi¬ 


polar  transmittance  in  the  partitioned  T 


IjkJL 


ik  of  Fig.  1 


can  be  avoided  by  using  unipolar  tranaalttance.  This  can  lead  how¬ 
ever  to  some  degradation  In  performance.  A.  systematic  numerical 
simulation  study*  of  a  neural  net  CAM  in  which  statistical  evalua¬ 
tion  of  the  performance  of  the  CAM  for  various  types  of  asaory  aasks 
(multivalued,  clipped  ternary,  clipped  u.b.)  and  thresholding  schemes 
(zero  threshold,  adaptive  threshold  where  energy  of  input  vector  is 
used  as  threshold,  adaptive  thresholding  and  relaxation)  was  carried 
out.  The  results  indicate  that  a  u.b.  aeaory  aask  can  be  used  with 
virtually  no  sacrifice  In  CAM  performance  when  the  adaptive  thresh¬ 
olding  and  relaxation  achene  is  applied.  The  scheme  assumes  an 
adaptive  threshold  Is  used  that  is  proportional  to  the  energy  (total 
light  intenalty)  of  the  input  entity  displayed  by  the  LED  array  at 
any  tlae.  In  the  scheme  of  Fig.  1(b)  this  can  be  realized  by  pro¬ 
jecting  an  image  of  the  input  pattern  directly  onto  an  additional  PD 
element.  The  PD  output  being  proportional  to  tha  total  intensity  of 
the  input  display  is  used  aa  a  variable  or  adaptive  threshold  in  a 
comparator  against  which  the  outputs  of  the  PD  elements  positioned 
behind  the  partitioned  components  of  the  memory  aask  are  com¬ 

pared.  The  outcomes,  now  bipolar,  are  attenuated  and  each  is  fad 
into  a  limiting  amplifier  with  delayed  feedback  (relaxation).  Each 
llalter /amplifier  output  Is  used  to  drive  the  LED  that  each  photo- 
detector  la  paired  with.  It  was  found5  that  this  schwas  yields  per- 
foraance  equivalent  to  that  of  aa  Ideal  CAM  with  multivalued  connec¬ 
tivity  aatrlx  and  zero  thresholding.  Mote  that  although  the  ini¬ 
tial  lzalng  2-D  entity  b^0^  la  unipolar  binary,  the  entitles  fed 

back  after  adaptive  thresholding  and  limited  amplification  to  drive 
the  LED  array  would  initially  be  analog  resulting  in  multivalued 
Iterates  and  intensity  displays.  However,  after  few  iterations  the 


outputs  become  binary  assuming  tha  extreme  valuas  of  tha  liaitar.  Tha 
ability  to  usa  u.b.  memory  aatricas  In  tha  fashion  dascribad  naans 
that  s lapis  black  and  white  photographic  transparanclas  or  binary 
SIMs  can  ba  usad  respectively  as  stationary  or  programmable  synaptic 
connectivity  masks  as  suggested  by  Fig.  1. 

SINOGRAM  CLASSIFIERS  AND  HETERQASSOCIATIVE  STORAGE 

Sinograms  are  object  rapraaantationa  encountered  In  tomography* * ^ 
They  are  also  useful  as  object  classifiers  specially  when  tha  objects 
are  point-like  and  dilute8.  Given  a  sat  of  2-D  dilute  objects  tha 
lannlag  distances  between  their  sinogram  classifiers  will  be  greater 
them  tha  Bamming  distances  batuaan  the  objects  themselves,  with  both 
sets  digitised  to  the  same  number  of  pixels,  making  it  easier  for  an 
associative  aanory  to  distinguish  between  the  sinograms8.  Sinogram 
classifiers  have  additional  advantages  that  enable  scale,  rotation, 
sad  shift  invariant  recognition  of  radar  target  which  can  not  be 
detailed  here  because  of  limited  space.  A  sinogram  la  a  cartesian 
plot  of  the  polar  projections  of  object  detail.  For  example  referring 
to  Fig.  2(a)  which  represents  a  dilute  object  consisting  of  16  points 
on  a  32x32  pixel  grid,  the  distance  that  the  projection  of  each  point 
makes  on  the  y  axis  as  meeeured  from  the  origin  whan  the  object  is 
rotated  about  the  origin  traces  a  sinusoidal  pattern  when  plotted 
against  rotation  angle  as  shown  in  Fig.  2(b).  Figure  2(c)  is  a 


rotation  $  results  in  a  sinogram  classifier  or  feature  space  of  the 
Car get  which  characterizes  it  at  any  fixed  elevation  angle  8.  The 
top  row  of  Fig.  3  shows  three  digitized  sinogram  classifiers  of 
scale  models  of  three  aerospace  targets  plotted  on  a  32x32  pixel 


grid.  These  art  treated  as  a  learning  sat  and  stored  heesro-associa- 
tively  rather  than  autoassociatively  be  replacing  in  eq.  (1)  by 

/  v  JCX 

k, 1*1, 2. .,.32;  a-1,2,3  where  r^'  represents  abbreviated  word 

labels  shom  In  the  botton  row  of  Fig.  3  with  which  the  three  test 
objects  are  to  be  associated. 


Fig.  3.  Hetero-esaoclatlve  storage.  Digitized  sinograms  (top)  and 
associated  word  labels  (bottom). 


RESULTS 

Representative  results  of  numerical  simulation  of  exercising  the 
hateroassoclatlva  memory  matrix  with  complete  and  partial  versions  of 
one  of  the  stored  entitles  in  which  the  fraction  n  of  correct  bits  or 
pixels  in  the  partial  versions  ranged  between  1  and  .1  are  presented 
in  Fig.  4.  Reliable  recognition  was  found  to  occur  after  one  iteration 
for  all  entitles  stored  down  to  n  •  .2.  For  n  •  .1  or  less  successful 
recall  of  correct  labels  was  found  to  depend  on  the  angular  location  of 
the  partial  data  the  memory  is  presented  with  as  illustrated  in  the  two 


u  •  100X  7f  •  601  77-  20X  77-  10X  77 -  10X 


Fig.  4.  Example  of  recognition  from' partial  information.  Complete 

sad  partial  sinograms  of  data  set  3  used  as  input  (top) ,  and 
final  memory  state-recognized  label  (bottom) . 


right-most  examples  In  Fig.  4.  Here  the  memory  could  not  label  the 
partial  input  correctly  but  converged  Instead  Onto  a  label  chat  It 
did  not  learn  before*  This  appears  to  be  a  generalization  (mixture) 
of  the  three  entitles  stored  earlier.  This  is  quite  analogous  to  the 
generalization  capability  of  tha  brain.  Note  the  generalization  is 
contrast  reversed  as  vs  know  that  stable  states  of  a  memory  with 
symmetric  connectivity  matrix  are  not  only  the  entitles  stored  but 
also  their  compliments. 

CONCLUSIONS 

Architectures  for  optical  implementation  of  2-0  neural  nets 
baaed  on  partitioning  of  the  4-0  connectivity  matrix  are  shown  to  be 
suitable  for  use  with  2-0  object  classifiers  or  feature  spaces.  An 
example  of  their  utility  in  super-resolved  recognition  (labeling)  of 
radar  targets  characterized  by  sinogram  classifiers  is  presented. 

The  results  show  that  neural  net  models  and  their  opto-electronlc 
analogs  furnish  a  new  viable  approach  to  signal  processing  land  auto¬ 
mated  recognition  that  is  both  robust  and  fault  tolerant. 
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A  NEW  ITERATIVE  ALGORITHM  FOR  EXTRAPOLATION  OF 
DATA  AVAILABLE  IN  MULTIPLE  RESTRICTED  REGIONS 
WITH  APPLICATION  TO  RADAR  IMAGING 

HnahJyh  Li,  Nabil  H.  Farhat  and  Yuhxyan  Shan 

TTm  Mom  Irhnni  nf  rknrir d  Enfinwrim 
Uahenky  o t  Pmylwii 
PbOmUpUa,  PA.  19104 


L  INTRODUCTION 

It  is  well  known  that  die  resolution  of  microwave  diversity  imaging  systems  [1]  depends 
on  the  spectral  and  angular  (aspect  related)  windows  and  on  the  use  of  polarization 

diversity  to  enhance  the  amount  of  information  about  a  object  collected  by  a  highly 

thinned  recording  aperture.  To  obtain  die  range  information  of  the  target,  one  can  use  a  pulsed 
signal  analyzed  in  die  time  domain  and  map  thus  the  range  profile  of  the  target  as  function  of 
aspect  angle  or  use  a  broad-band  CW  signal  analyzed  in  frequency  domain  to  yield  its  fre¬ 
quency  response.  The  range  resolution  is  inversely  proportional  to  the  bandwidth  coverage  of 
the  measurement  system.  The  latter  approach  is  ordinarily  n»nmr  to  implement,  in  practical 
situations,  however,  due  to  limitation  of  the  measurement  system  or  restriction  of  bandwidth 
allocation,  the  observed  data  can  lie  in  multiple  restricted  spectral  regions  which  we  call  pass 
bands.  Several  methods  of  extrapolating  the  measured  data  beyond  the  observed  regions  have 
been  proposed  and  tested  CL3.4]  in  an  attempt  to  achieve  the  full  resolution  of  the  unrestricted 
spectral  range,  when  a  prior  knowledge  of  the  maximum  dimension  of  the  object  exists,  and 
an  iterative  procedure  is  applied. 

To  increase  die  resolution  obtained  from  spectral  data  of  such  limited  extent,  techniques 
of  nonlinear  power  spectrum  estimation  have  been  used  with  notable  success  [5].  These 
include  autoregression  (AR),  linear  prediction  (LP),  and  maximum  entropy  method  (MEM). 
For  a  stationary  Gaussian  process,  the  above  methods  can  be  shown  to  be  equivalent  [5). 


Although  the  spectra  estimated  by  MEM  or  AR  can  be  very  sharp  and  well  resolved,  this 
may  not  be  to  advantage  in  a  microwave  imaging  system.  If  the  data  are  not  sampled  densely 
enough  in  the  spectral  domain  as  the  sharp  well  resolved  component  may  be  missed,  and  the 
results  may  not  reflect  faithfully  the  actual  spectral  amplitudes.  Besides,  image  reconstruction 
from  microwave  diversity  imaging  systems  involves  coherent  superposition  of  the  data  in  the 
spectra,  or  range  profiles,  of  die  scatterer  (obtained  at  different  looks  or  aspect  angle),  where 
these  are  estimated  from  partial  data  available  in  segmented  bands  [1].  If  the  estimated  ampli¬ 
tudes  of  the  range  profiles  obtained  by  MEM  or  AR  depart  from  the  desired  values  because  of 
undersampling,  image  degradation  will  result  Therefore,  to  overcome  die  dense  sampling 
requirement  it  may  be  preferable  to  extrapolate  the  data  available  in  the  various  pass  bands 
into  the  vacant  bands  before  the  spectra  or  range  profile  are  formed,  and  image  reconstruction 
is  hwHm*iIwii 

If  one  wishes  to  predict  the  next  value  in  a  sequence  on  the  basis  of  previous  samples, 
one  can  use  the  linear  prediction  model  which  consists  of  finding  the  parameters  that  minim. 
ize  die  summation  of  prediction  error  energies  which  will  be  defined  in  the  later  section.  To 
extrapolate  the  data  beyond  the  observed  region,  an  intuitive  way  is  to  predict  die  exterior 
data  by  using  the  same  parameters  obtained  by  the  linear  prediction  model.  The  most  popular 
approach  for  linear  prediction  parameters  estimation  with  N  data  samples  is  the  Burg  algo¬ 
rithm  [5,6].  For  a  given  number  of  data  samples  in  a  given  observation  interval,  in  order  to 
separate  the  discrete  spectra  (in  this  paper,  spectrum  is  defined  as  the  Fourier  transform  of  the 
observed  data),the  required  model  order  in  the  linear  prediction  method  increases  as  the 
separation  of  spectra  decreases,  Le,  it  is  easier  to  model  the  data  sequence  for  spectra  with 
larger  separation  than  those  with  closer  separation,  which  translates  into  well  separated 
scattering  centers.  In  addition,  for  a  given  model  order  and  given  number  of  sampling  points, 
it  is  easier  to  distinguish  the  two  close  spectra  components  (scattering  centers)  by  a  data  set 
with  longer  observation  interval  than  that  by  a  data  set  with  shorter  observation  interval.  It 
was  also  suggested  that  the  model  order  should  not  exceed  half  of  the  number  of  data  points 


for  short  data  segment  because  otherwise  die  linear  prediction  spectral  estimate  will  exhibit 
spurious  peaks  [5].  From  the  above  observations,  one  can  conclude  that  it  would  be  more 
difficult  to  resolve  two  closer  point  targets  (Fourier  transform  of  the  observed  data  in  fre¬ 
quency  domain)  with  short  data  band.  If  all  die  observed  data  within  multiple  restricted 
regions  can  be  fully  utilized,  better  resolution  can  be  expected. 

It  is  known  that  the  linear  prediction  method  is  especially  suited  for  those  cases  when  die 
spectra  are  in  discrete  type.  Under  high  frequency  condition,  the  scattered  fields  of  a  complex 
target  can  be  ascribed  to  a  few  discrete  scattering  centers  and  edge  diffractions.  It  will  be 
shown  that  under  die  high  frequency  approximation  the  locations  of  the  scattering  centers  and 
their  scattering  strengths  are  independent  of  the  operating  frequency  for  a  given 
transmittertieceiver  pair.  These  phenomena  provide  the  motivation  to  investigate  a  new 
method  for  extrapolation  of  scattered  field  data  available  in  multiple  restricted  frequency 

« - a- 

MIKU. 

In  this  paper,  a  new  iterative  algorithm  which  uses  the  Burg  algorithm  to  find  the  linear 
prediction  parameters  and  an  iterative  method  to  modify  the  prediction  parameters  is  proposed 
and  tested  with  both  simulated  and  realistic  measured  data  generated  in  our  anechoic  chamber 
experimental  microwave  imaging  and  measurement  facilities.  With  this  algorithm,  one  can 
obtain  acceptable  extrapolation  beyond  the  observed  region  if  the  spectra  are  in  discrete  forms 
and  die  separatum  of  the  spectra  are  not  too  dose.  Both  simulation  and  experimental  results 
are  presented  to  demonstrate  as  an  example  the  effectiveness  of  the  method  in  microwave 
diversity  radar  imaging. 

IL  THE  NEW  ITERATIVE  ALGORITHM 

The  most  popular  approach  to  linear  prediction  parameter  estimation  with  N  data  samples 
{  xi , . . . ,  Xff  }  was  introduced  by  Burg  [6].  The  linear  prediction  parameters  are  obtained  by 
mhrinrizmg  the  sum  of  the  forward  and  backward  prediction  error  energies  , 
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sobject  to  die  constraint  that  the  prediction  parameters  satisfy  a  recursion  relationship  [5], 
where  is  die  forward  prediction  error  with  model  order  p  and  is  given  by 

m  ^£tapkxn-k*  (2) 

*■0 

and  bp*  is  the  backward  prediction  error  with  model  order  p  and  is  given  by 

V*  £«**•-*♦*•  0) 

km  0 

are  called  linear  prediction  parameters,  and  asterisk  denotes  the  complex  conjugate  opera¬ 
tion. 

If  one  is  going  to  extrapolate  from  the  available  data  beyond  die  observed  region,  a 
straight  forward  way  is  to  use  die  estimated  prediction  parameters  and  the  measured  data 
by  the  following  equations, 

*N+jm  ft,^pkxN+j~k  J>  0;  (4) 

*-j*  ;>0,  (3) 

km  1 

where  the  A  denotes  the  estimated  value. 

If  die  data  available  are  confined  to  multiple  separate  spectral  regions  or  pass-bands  of 
equal  width  as  illustrated  in  Fig.  1,  and  one  tries  to  extrapolate  from  the  observed  data  to  the 
vacant  bands,  an  intuitive  method  is  to  divide  the  inner  vacant  band  into  two  parts  of  equal 
width  and  to  extrapolate  into  the  left  part  by  using  the  prediction  parameters  obtained  from 
the  data  set  of  region  I  and  extrapolate  into  the  right  part  by  using  the  model  parameters 
obtained  from  the  data  set  of  region  H 

If  the  data  sequence  can  be  correctly  expressed  by  the  prediction  parameters,  then  the 
extrapolation  error,  which  is  defined  as  the  absolute  value  of  the  complex  difference  between 
the  theoretical  values  and  extrapolated  values,  would  be  very  small.  However,  if  the  prediction 
parameters  cannot  model  the  sequence  correctly,  the  error  of  extrapolation  may  accumulate. 
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We  have  found  that  the  longer  the  available  data  string,  the  larger  die  model  order,  the  linear 


prediction  model  which  characterizes  the  data  sequence  would  be  more  accurate  especially  in 
the  presence  of  noise.  However,  the  model  order  should  not  exceed  half  of  the  number  of  the 
samples  because  the  estimated  spectrum  will  produce  spurious  peaks  [5]. 


In  order  to  utilize  the  information  available  in  different  regions,  a  new  iterative  algorithm 


using  the  Burg  algorithm  to  estimate  the  prediction  parameters  and  a  new  iterative  method  is 
proposed.  The  procedure  illustrated  in  Fig.  2  is  as  follows: 


1.  Divide  the  inner  vacant  band  into  two  parts  of  equal  width.  Extrapolate  into  the  left 


part  by  using  die  prediction  parameters  obtained  from  the  data  set  of  region  I  and 


extrapolate  into  die  right  part  by  using  the  prediction  parameters  obtained  from  the 
data  set  of  region  IL  If  the  bands  are  not  equal  in  width,  unequal  division  of  die 
vacant  intervening  bands  may  be  appropriate. 


2.  With  the  "vacant  bandV  data  together  with  the  observed  data,  use  the  Burg  algo¬ 


rithm  to  find  a  new  set  of  prediction  parameters. 


3.  Using  this  set  of  prediction  parameters  and  die  data  of  region  I,  extrapolate  into  the 
left  part  of  vacant  bands,  and  using  the  same  set  of  prediction  parameters  and  the 
data  in  region  n,  extrapolate  into  the  right  part  of  the  vacant  bands. 


4.  Using  this  set  of  parameters  together  with  the  extrapolated  data,  estimate  the  data  in 


the  observation  region  I  and  IL  Calculate  the  error  energy  between  the  measured 


data  and  the  estimated  data  in  the  observation  regions.  The  error  energy  is  denoted 


by  £t  and  is  given  by: 


i  i 


where  x*  are  the  measured  data,  are  the  forward  estimation  of  x, ,  i  \  are  the 


backward  estimation  of  xj ,  is  the  forward  prediction  error,  and  bt  is  the  backward 

prediction  error. 
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5.  With  the  measured  data  together  with  the  estimated  "vacant  band'Y  data,  use  the 
Burg  algorithm  to  find  a  new  set  of  prediction  parameters.  From  the  measured  data 
and  this  new  set  of  prediction  parameters,  extrapolate  die  "vacant  band'Y  data  as 
described  in  step  3. 

6.  Use  the  same  procedure  of  step  4  to  calculate  the  new  error  energy  of  the  pass- 
bands,  call  it  Ej. 

7.  Compare  £t  with  Bj,  if  £2  i*  smaller  than  Elt  replace  the  error  energy  £x  by  £2  , 
repeat  step  5. 

8.  If  £2  is  greater  than  Elf  stop  the  iteration,  and  take  the  extrapolated  data  of  the  pre¬ 
vious  loop  as  the  final  result 

In  step  1,  if  die  width  at  a  single  band  (band  I  and/or  band  II)  is  not  large  enough,  the 
extrapolation  errors  produced  by  the  prediction  parameters  obtained  from  single  pass-band 
data  may  be  very  large,  in  that  case,  we  can  set  the  data  in  the  vacant  bands  to  zero. 

The  above  iterative  method  can  be  easily  applied  to  the  case  where  only  one  single  data 
band  is  available  to  extrapolate  outside  it  The  procedures  are  almost  the  same  except  that 
only  one  data  band  sequence  is  used  to  extrapolate  to  the  exterior  bands  and  to  calculate  the 
extrapolation  errors. 

HL  SCATTERING  PROPERTIES  OF  A  METALLIC  OBJECT 

In  this  section  we  shall  show  that  under  the  high  frequency  approximation  the  scattered 
fields  of  a  metallic  object  can  be  expressed  as  superposition  of  scattered  fields  of  discrete 
scattering  centers.  These  phenomena  allow  us  to  apply  the  proposed  extrapolation  algorithm 
to  radar  imaging. 

For  a  object  large  compared  with  wavelength,  the  scattering  mechanism  can  be 

divided  into  the  following  components  [7]: 


1.  Specular  scattering  points 

2.  Scattering  from  surface  discontinuity:  edges,  context,  dpt  etc. 

3.  Scattering  from  surface  derivative  discontinuities 


4.  Cscepmg  waves 

5.  Traveling  waves 

6.  Scattering  from  coocave  repoos 

7.  Multiple  scattering  points 

For  most  situations,  tbe  major  contributions  to  the  scattered  waves  are  ascribed  to  die 
specular  scattering  points  and  edge  diffractions. 

i.onnocr  ft  mwiiiic  oojoci  wiw  oo  ft  rocaipo  procirn  m  munuyinpa  uy  ft  punc  wito 
as  shown  in  Fig.  3.  The  distance  between  tbe  rotation  center  O  and  the  transmitter  and 
receiver  are  Rt  and  R,  respectively,  and  the  unit  vectors  in  the  direction  of  transmitter  and 
receiver  are  It  and  4  respectively.  Under  the  physical  optics  and  Bom  approximations,  the 


£(*)  -  L2* <** ) x  «  (*>  dS ' , 


(7) 


where  &  is  the  wave  number,  Sa  the  flluminated  region,  A  ( f )  the  unit  normal  vector  at  the 
surface  point  7**,  and  ft0  (k)  the  incident  magnetic  field  at  the  rotation  center. 

As  k  approaches  infinity,  the  asymptotic  expression  of  the  above  equation  can  be 
obtained  by  applying  the  stationary  phase  method  [8]  to  eq.  (7).  The  result  is. 


(«> 


where  1**  >  are  the  vectors  such  that 


0 » 


(9) 


and 


(10) 
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r-r, 

where  *r— -7  is  the  derivative  with  respect  to  the  surface  curvature.  The  points  correspond* 

ing  to  the  solutiocs  of  eq.  ( 9 )  are  called  stationary  points  or  equi-phase  points  or  the  scatter* 
inf  centers;  the  term  4  ( "f*j )  x  ff0  (k)  t^lSj  ■  aj  is  called  the  scattering  strength  for  die  partic¬ 
ular  scattrring  center  at  "fj  .  It  is  seen  that  the  locations  of  the  scattering  centers  depend  on 
the  directions  of  4 »  4  **  well  as  the  shape  of  the  metallic  surface.  The  scattering  strength 
depends  on  the  local  properties  of  the  scattering  centers.  The  above  analysis  illustrates  that  the 
object  function  we  would  be  dealing  with  in  high  frequency  radar  imaging  are  of  discrete 
form  consisting  of  point  scattering  centers. 

If  the  received  scattered  Adds  have  been  calibrated  with  a  reference  target  [1],  the 
corrected  Adds  (£)  can  be  expressed  as 

^  «*K*“*) 'r'j]  •  (11) 
J 

The  Fourier  transform  of  eq.  (11)  will  give  the  range  profile  and  scattering  strength  of 
the  scattering  centers. 

IV.  RESULTS 

In  this  section,  die  performance  of  the  proposed  new  algorithm  using  both  simulated  and 
realistic  data  will  be  evaluated.  First,  assume  for  simplicity  an  object  consisting  of  n  point 
scatteren  located  at  (r9  +  yj)  is  illuminated  by  a  plane  wave,  where  r9  is  the  distance 
between  the  tranamilterirecdver  and  a  reference  point  of  the  object  and  y;-  is  the  differential 
range  of  the/th  scatterer  (range  relative  to  rj.  Under  far  field  condition  and  ignoring  multiple 
scattering,  and  considering  for  simplicity  a  version  of  eq.  (11),  the  corrected  scalar 

field  can  be  expressed  as 

E 's  (*)  -  Z  <*j  •J2k'h 
J 


(12) 


In  the  following  simulation,  the  theoretical  values  of  E (k)  are  calculated  in  200 
equally  spaced  frequency  steps  covering  the  frequency  range  /i  »  6  GHz  to  fa  *  16  GHz, 
with  an  assumed  signal  to  noise  ratio  set  to  40  dB.  These  values  anticipate  the  realistic  experi¬ 
mental  data  in  testing  the 

Assume  die  available  (computed)  data  are  in  the  following  pass-band  (/30  ,/M )  and 
(fa  >fa  )•  We  want  to  extrapolate  the  data  to  the  vacant  bands  (fi  ,fo ),  (fn  ,fa  ),  and 
(/m  00  )•  The  range  resolution  obtained  by  the  DFT  method  using  the  whole  bandwidth 

(fi  ,fa  )  is  about  L5  cm.  The  resolution  using  a  single  frequency  band  is  about  5 Jem.  The 
resolution  using  both  frequency  bands  is  about  2.0  cm,  however,  very  high  side-lobe  level  will 
be  produced.  We  consider  a  scatterer  with  seven  point  scatteren,  the  location  and  scattering 
strength  for  each  point  scattenrare  ( *  -  30  cm,  ax  *  0J),  (  r2  *  -  20  cm,  <*2*0.5), 
(  r3  ■  -  10  cm,  oj  »  (X5),  (  r4  *  -  2  cm,  44  *  1),  (  r5  *  10  cm,  45  *  0 .25),  ( r6  »  20  cm, 

4s  *  0 25),  (  r?  *  -  30  cm,  07  *  0.25).  The  values  of  the  field  at  each  sampled  frequency  f} 

* 

are  calculated  using  eq.  (12). 

Define  the  extrapolation  enor  at  frequency  fj  as 

(13) 

whree  E  \  (fj)  ii  the  extrapolated  nine  at  each  fj .  The  extrapolation  errors  for  different 
algorithms  are  compared  and  shown  in  Fig.  4(a).  The  bold  solid  curve  is  the  amplitude  of  the 
theoretically  computed  fields  E  \  (fj ),  the  thin  solid  curves  are  the  extrapolation  error  after 
100  iterations  using  die  algorithm  proposed  in  [3],  die  dashed  curves  are  obtained  by  using 
die  Burg  algorithm  to  find  the  prediction  parameters  from  the  respective  pass-band,  and  using 
this  set  of  parameters  together  with  data  in  each  pass-band  to  extrapolate  to  the  outside 
regions  (bands  m  and  IV).  The  dotted  line  curves  are  obtained  using  this  new  algorithm  with 
one  iteration  and  with  model  order  25.  The  algorithm  proposed  in  [3]  involves  basically  appli¬ 
cation  of  the  Gerchberg  algorithm  to  data  in  the  multiple  restricted  regions.  However,  no 
numerical  or  experimental  results  are  given  in  that  paper.  It  is  clear  from  the  results  obtained 


-  10- 


here  Out  the  algorithm  in  [3]  seems  not  to  be  effective  in  the  case  considered  as  the  errors 
can  exceed  die  amplitude  of  die  theoretical  fields.  Extrapolation  from  single  pass-band  are  not 
good  in  this  example,  because  the  model  order  is  not  sufficient  to  model  the  riant  series  in  the 
presence  of  noise.  The  proposed  new  method  after  one  iteration  is  seen  to  produce  small  error. 

The  Fourier  transform  (FT)  of  the  all-band  data  (Le,  data  in  region  I  to  IV),  pass-band 
data  only,  pass-band  0ns  extrapolated  data  using  the  algorithm  in  [3]  and  pass-band  plus 
extrapolated  data  with  the  new  proposed  method  are  shown  in  Fig.  4(b)  to  4(e)  respectively. 
Note  that  the  FT  of  spectral  data  yields  range  profile  of  the  scattering  object  It  is  clear  that 
FT  using  pass-band  data  only  (Fig.  4(c))  has  very  high  side-lobe  structure,  the  FT  of  the 
extrapolated  data  using  the  algorithm  in  [3]  (Fig.  4(e))  is  totally  different  from  the  original  of 
Fig.  4(b).  The  result  obtained  by  Fourier  transforming  die  data  generated  by  the  proposed 
algorithm  is  shown  in  Fig.  4(d)  which  exhibits  excellent  agreement  with  the  all-band  result  of 
Fig.  4(b).  The  magnitudes  of  the  peaks  in  Figs.  4(b)  and  4(d)  depart  from  die  original 
assigned,  values  because  of  zero  padding  used  in  the  fast  Fourier  transform  (FFT)  algorithm. 
Thu  lack  of  fidelity  in  scattering  strength  reconstruction  does  not  have  a  discernible  degrading 
effect  on  the  quality  of  image  reconstructed  as  will  be  illustrated  below,  but  is  important  and 
must  be  dealt  with  when  quantitative  analysis  of  scattering  strengths  is  needed. 

If  the  frequency  coverage  is  increased  to  (/j  *  6  GHz./gx)  *  20  GHz)  with  the  number 
of  sampling  points  being  fixed  to  200  and  the  pass-bands  are  kept  at  (/30  ,  fn  )  and  (fm  , 
/i7o),  the  computed  fields  and  the  extrapolation  errors  would  be  as  shown  in  Fig.  5  .  It  is  seen 
that  the  extrapolation  error  indicated  by  the  dashed  line  becomes  smaller.  If  the  frequency 
coverage  is  decreased  to  (/1  *  6  GHz,/^  *  12  GHz),  the  results  would  be  as  shown  in  Fig. 
&  It  is  seen  that  die  extrapolation  errors  indicated  by  the  dashed  and  dotted  curves  are  now 
both  high.  The  FFT  of  the  whole  band  data,  pass-band  data  only,  and  the  extrapolated  plus 
pass-band  data  using  this  method  are  shown  in  Figs.  6(b)  to  6(d)  respectively.  The  results  in 
Figs.  5  and  6  indicate  the  desirability  of  using  segmented  spectral  data  spanning  wider  spec¬ 
tral  ranges. 
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The  performance  of  the  algorithm  using  realistic  data  is  also  evaluated.  The  test  object,  a 
metalized  100:1  scale  model  of  a  B-52  aircraft  with  79  cm  wing  span  and  68  cm  long 
foaelage  was  mounted  on  a  computer-controlled  elevation-over-azimuth  positioner  situated  in 
an  anechoic  chamber  environment.  201  equal  frequency  steps  covering  the  /t  ■  6.1  to 
fax  m  17  J  GHz  range  were  used  to  obtain  the  frequency  response  of  the  object  as  described 
in  [1].  The  target  is  positioned  for  a  fixed  elevation  angle  of  30*  while  the  trimuih  angle 
was  altered  between  0*  and  90*  in  steps  of  0.7*  for  a  total  of  128  angular  looks. 

The  pass-band  is  first  defined  ss  (/»  ./jq)  and  (fm  ,/l70).  The  measured  values  and  the 
extrapolated  errors  of  the  broad-tide  look  which  is  90*  from  die  head-on  lode  are  shown  in 
Fig.  7(a).  The  solid  line  curve  is  the  amplitude  of  die  range-phase  corrected  field  (see  [1]). 
The  dashed  curve  wpwsmts  the  extrapolation  enor  resulting  from  extrapolating  from  each 
single  band  (bands  I,  II)  with  model  order  25  ss  described  in  step  1  of  die  proposed  algo¬ 
rithm.  The  dotted  line  curves  are  obtained  using  the  new  algorithm  with  1  iteration  and 
model  order  25.  The  extrapolation  enor  for  measurement  is  defined  in  a  manner  similar  to  the 
definition  of  enor  in  numerical  shnularirm  as  the  of  the  difference  between  the 

corrected  measured  Adds  and  extrapolated  Adda.  The  Fourier  transfonn  from  the  whole  band 
data,  the  pass-band  data  only,  and  the  pass-band  together  with  extrapolated  data  are  shown  in 
Figs.  7(b),  7(c)  and  7(d)  respectively.  Fourier  transform  of  die  corrected  scattered  fields  will 
give  the  range  profile  of  the  target  in  that  view,  hi  this  figure,  it  is  seen  that  the  extrapolation 
errors  do  not  improve  after  one  iteration.  The  reason  can  be  explained  from  the  plot  of  the 
range  profile  shown  in  Fig.  7(b).  hi  this  view  direction,  the  major  contributions  to  the  scat¬ 
tered  fields  are  due  to  fuselage  and  primarily  those  engines  and  fuel  tank  which  are  on  the 
illuminated  side.  Specular  scattering  from  these  points  are  well  separated  in  time  or  distance 
and  their  number  is  small.  Hence  the  linear  prediction  parameters  obtained  from  single  pass- 
bend  are  sufficient  to  model  the  data  sequence.  The  extrapolation  errors  are  not  as  small  as 
those  obtained  by  simulations.  The  reason  of  tins  is  that  the  applicability  of  linear  prediction 
model  to  the  extrapolation  of  scattered  fields  of  a  metallic  object  is  based  on  the  high 


frequency  approximation.  In  the  measurement  data,  however,  polarization  effects,  edge 
diffraction,  multiple  scattering  and  the  Mure  to  satisfy  the  high  frequency  approximation  in 
the  lower  region  of  the  frequency  band  utilized  in  the  measurement  will  degrade  the  perfor¬ 
mance  of  die  algorithm. 

The  reconstructed  images  of  die  test  object  using  data  collected  in  an  angular  windows 
of  90°  extending  from  head-on  to  broadside  in  128  looks  (see  [1]  for  detail)  and  different  fre¬ 
quency  bands  are  shown  in  Fig.  8.  The  transmitting  antenna  is  right-hand  circularly  polarized 
and  the  receiving  antenna  is  left-hand  circularly  polarized,  which  constitutes  by  the  convention 
given  in  [7]  a  co-polarized  transmitting/receiving  system.  Figure  8(a)  is  obtained  by  using  die 
whole  band  data;  Fig.  8(b)  is  obtained  by  using  die  pass-band  data  alone.  Figures  8(c)  and 
8(d)  are  obtained  by  extrapolating  without  iteration  and  after  one  iteration  respectively.  The 
model  order  used  is  M  *  25  in  both  cases. 

If  the  pass-band  is  defined  as  5  ,/130),  the  reconstructed  images  obtained  by  using  the 
pass-band  data  alone  and  by  extrapolation  without  iteration  and  after  1  iteration  are  shown  in 
Figs.  8(e),  8(f),  and  8(g)  respectively.  The  model  order  used  is  also  Af  *  25. 

ft  is  seen  that  the  image  quality  of  Figs.  8(c)  and  8(e)  are  as  good  as  that  of  Fig.  8(a). 
These  results  show  the  effectiveness  of  the  application  of  the  proposed  algorithm  to  radar 
imaging  from  segmented  data  bands. 

V.  CONCLUSION 

A  new  method  employing  the  Burg  algorithm  and  an  iterative  procedure  to  extrapolate 
observed  data  beyond  restricted  regions  of  observation  has  been  proposed  and  tested.  Simula¬ 
tion  and  experimental  results  prove  the  effectiveness  of  this  proposed  method.  The  algorithm 
is  especially  effective  when  the  spectra  of  the  collected  data  (the  object  range  profile  in  this 
case)  are  in  discrete  form.  Possible  applications  of  this  new  method  can  be  found  in  diverse 
flykft  whenever  the  data  is  available  in  restricted  bands.  For  example,  in  multiple  band 
microwave  imaging  system,  the  quality  of  the  image  obtained  by  extrapolating  from  a  much 


smaller  bandwidth  can  be  as  good  as  that  obtained  by  data  in  the  full  bandwidth.  The  cost  of 
die  imaging  system  can  hence  be  reduced  drastically  as  the  cost  of  the  required  gear  can  be 
much  lower  than  die  cost  of  the  gear  to  cover  the  full  bandwidth  and  restrictions  on  use  of 
frequency  bands  can  be  accommodated. 

• 

ACKNOWLEDGEMENT 

This  research  was  supported  by  the  Air  Force  Office  of  Scientific  Research,  Air  Force 
system  Command,  under  grant  AFOSR-81-0240F  and  by  the  Army  Research  Office  under 
contract  DAAG  29-83-K-0 1 20-P02. 


REFERENCE 


[1]  N.H  Farhat,  CL.  Werner  and  Til  Chu,  "Prospects  for  Three-Dimensional  Projective  and 
Tomographic  Imaging  Radar  Networks,"  Radio  Science,  VoL  19,  No.  5,  pp.  1347-1355, 
1984. 

[2]  R.W.  Gerchberg,  "Super-Resolution  Through  Error  Energy  Reduction,"  Optica  Acta,  VoL 
21,  No.  9,  pp.  708-720,  1974. 

[3]  YoshHri  Yamakoshi  and  Takuso  Sato,  "Iterative  Image  Restoration  From  Data  Available 
in  Multiple  Restricted  Regions,"  Applied  Optics,  VoL  21,  No.  24,  pp.4473-4480,  1982. 

[4]  C.Q.  Lan,  ILK.  Xu  and  Glen  Wade,  "Limited  Angle  Diffraction  Tomography  and  Its 
Application  to  Planar  Scanning  Systems,"  IEEE  Trans,  on  Sonics  and  Ultrasonics,  VoL 
SU-32,  No.  1,  pp.  9-16, 198S. 

[5]  Steven  M.  Kay  and  Stanley  Lawrence  Marple,  "Spectrum  Analysis  -  A  Modem  Perspec¬ 
tive,"  Proc.  of  the  IEEE,  VoL  69,  pp.1380-1419,  November  1981. 

[6]  Burg,  John  P.,  "A  New  Analysis  Technique  for  Time  Series  Data”,  NATO  Advanced 
Study  Institute  on  Signal  Processing  with  Emphasis  on  Underwater  Acoustics,  pp.  12-23, 
August  1968. 

[7]  G.  T.  Rude,  es  aL,  Radar  Cross  Section  Handbook,  edited  by  G.  T.  Ruck,  Chapter  8,  Ple¬ 
num  Press,  New  York  1970. 

[8]  Jon  Mathews  and  RX.  Walker,  Mathematical  Methods  of  Physics,  2nd  Ed.,  W.A.  Benja¬ 
min,  Inc.,  Menlo  Path,  California,  1970. 


FIGURE  CAPTIONS 


Fig.  1  Available  data  in  multiple  regions.  Passbandftharicri  region) 
suxxoanded  by  vacant  bands. 

Fig.  2  Schematic  diagram  o f  the  proposed  new  iterative  extrapolation  method. 

Fig.  3  Geometry  of  the  scattering  measurement  system. 

Fig.  4  (a)  Magnitude  of  theoretical  fields  and  comparison  of  extrapolation 
errors  of  different  methods, /i  *  6  GHz,  fa  ■  16  GHz. 
magnitude  of  theoretical  fields, 
extrapolation  error  from  a  single  passband,  no  iteration, 
extrapolation  errors  from  hew  iterative  algorithm. 

_ extrapolation  errors  from  algorithm  proposed  in  [3]. 

(b)  FFT  of  the  whole  band  data. 

(c)  FFT  of  the  passband  data. 

(d)  FFT  of  the  passband  and  extrapolated  data  with  1  iteration. 

(e)  FFT  of  the  passband  and  extrapolated  data  using  algorithm  proposed  in  [3]. 


Fig.  5  Magnitude  of  theoretical  fields  and  comparison  of  extrapolation 
errors  with  and  without  iteration,  ft  ■  6  GHz,/^  *  20  GHz. 

—  magnitude  of  theoretical  fields. 

■ . -  extrapolation  error  from  respective  passband,  no  iteration. 

_ extrapolation  errors  from  new  iterative  algorithm. 


Fig.  6  (a)  Magnitude  of  theoretical  fields  and  comparison  of  extrapolation 
errors  with  and  without  iteration,  fx  *  6  GHz,  =  12  GHz. 

■ . ■mi  mggnifnA»  of  theoretical  fields. 

— — —  extrapolation  error  from  respective  passband,  no  iteration, 
extrapolation  errors  from  new  iterative  algorithm. 

(b)  FFT  of  the  whole  band  data. 

(c)  FFT  of  die  passband  data. 

(d)  FFT  of  the  passband  and  extrapolated  data  with  1  iteration. 


Fig.  7  (a)  Magnitude  of  die  measured  fields  and  comparison  of  extrapolation 
errors  without  and  with  1  iteration 
■■■—i  magnitude  of  theoretical  fields. 

- extrapolation  error  from  respective  passband,  no  iteration. 

extrapolation  errors  from  new  iterative  algorithm. 

(b)  FFT  of  the  whole  band  data. 

(c)  FFT  of  the  paisband  data. 

(d)  FFT  of  the  passband  and  extrapolated  data  with  1  iteration. 
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Hg.  8  Reconstructed  images  of  the  metalized  scale  B-52  aircraft 

an  angular  window  of  90°  extending  from  head-on  to  broad-side  in  128 
looks  and  different  spectral  coverage. 

Reconstructions  from: 

(a)  entire  bandwidth  (fx  ,fm\ 

(b)  passband  (/» (fm  ,/170). 

(c)  passband  (/*  (fm  ,/170) 

and  extrapolation  data  (extrapolated  data  into  empty  bands)  without  iteration. 

(d)  passband  (/*,  ,/go),  (fm  ,/170) 

and  extrapolation  data  with  1  iteration. 

(e)  passband  (fa  ,/&$. 

(f)  passband  (fa  ./^  and  extrapolation  data  without  iteration. 

(g)  passband  (fa  ,/130)  and  extrapolation  data  with  1  iteration. 
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ARCHITECTURES  FOR  OPTO-ELECTRONIC  ANALOGS  OF  SELF-ORGANIZING  NEURAL  NETWORKS 
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Abstract 


Architectures  for  partitioning  opto-electronic  analogs  of  neural  nets 
into  Input/output  and  internal  units  to  enable  self-organization  and 
learning  where  a  net  can  form  its  own  internal  representations  of  the 
"environment"  are  described. 


1.  INTRODUCTION:  In  our  proceeding  work  on  optical  analogs  of  neural  nets, 
[1], [2],  the  nets  described  were  programmed  to  do  a  specific  computational 
task,  namely  a  nearest  neighbor  search  by  finding  the  stored  entity  that  is 
closest  to  the  address  in  the  Hamming  sense.  As  such  the  net  acted  as  a 
content  addressable  associative  memory.  The  programming  was  done  by 
computing  first  the  interconnect lvity  matrix  using  an  outer-product  recipe 
given  the  entities  we  wished  the  net  to  store  and  become  familiar  with 
followed  by  setting  the  weights  of  synaptic  Interconnections  or  links 
between  neurons  accordingly. 

In  this  paper  we  are  concerned  with  architectures  for  opto-electronic 
implementation  of  neural  nets  that  are  able  to  program  or  organize 
themselves  under  supervised  conditions,  l.e.,  of  nets  that  are  capable  of 
(a)  computing  the  interconnectivity  matrix  for  the  associations  they  are 
to  learn,  and  (b)  of  changing  the  weights  of  the  links  between  their  neurons 
accordingly.  Such  self-organizing  networks  have  therefore  the  ability  to 
form  and  store  their  own  internal  representations  of  the  entities  or 
associations  they  are  presented  with. 

Multi-layered  self-programming  nets  have  been  described  recently  [3]~ 
[5]  where  the  net  is  partitioned  into  three  groups.  Two  are  groups  of 
visible  or  external  Input/output  units  or  neurons  that  interface  with  the 
outside  world  l.e.,  with  the  net  environment.  The  third  is  a  group  of 
hidden  or  Internal  units  that  separates  the  input  and  output  units  and 
participates  in  the  process  of  forming  internal  representations  of  the 
associations  the  net  is  presented  with,  as  for  example  by  "clamping"  or 
fixing  the  states  of  the  input  and  output  neurons  to  the  desired 
associations  and  letting  the  net  run  through  its  learning  algorithm  to 
arrive  ultimately  at  a  specific  set  of  synaptic  weights  or  links  between  the 
neurons  that  capture  the  underlying  structure  of  all  the  associations 
presented  to  the  net.  The  hidden  units  or  neurons  prevent  the  input  and 
output  units  from  communicating  with  each  other  directly.  In  other  words  no 
neuron  or  unit  in  the  input  group  is  linked  directly  to  a  neuron  in  the 
output  group  and  vice-versa.  Any  such  communication  must  be  carried  out  via 
the  hidden  units.  Neurons  within  the  input  group  can  communicate  with  each 
other  and  with  hidden  units  and  the  same  is  true  for  neurons  in  the  output 
group.  Neurons  in  the  hidden  group  can  not  communicate  with  each  other. 
They  can  only  communicate  with  neurons  in  the  input  and  output  groups  as 
stated  earlier. 


Two  adaptive  learning  procedures  In  such  partitioned  nets  have 
attracted  considerable  attention.  One  is  stochastic  Involving  a  simulated 
annealing  process  [6],C7]  and  the  other  Is  deterministic  involving  an  error 
baok'-propogatlon  process  [<*].  There  Is  general  agreement  however;  that 
because  of  their  Iterative  nature,  serial  digital  computation  of  the  links 
with  these  algorithms  Is  very  time  consuming.  A  faster  means  for  carrying 
out  the  required  computations  is  needed.  Never-the-less  the  work  mentioned 
represents  a  milestone  in  that  It  opens  the  way  for  powerful  collective 
computations  in  multilayered  neural  nets  and  in  that  It  dispels  earlier 
reservations  C8]  about  the  capabilities  of  early  models  of  neural  nets  such 
as  the  Perceptron  [9]  when  the  partitioning  concept  Is  Introduced.  What  Is 
most  significant  and  noteworthy,  in  our  opinion.  Is  the  ability  to  now 
define  buffered  input  and  output  groups  with  unequal  number  of  neurons  In  a 
net  which  was  not  possible  with  earlier  nets  where  all  neurons  participate 
in  defining  the  initial  (input)  and  final  (output)  states  of  the  net. 


2.  ANALOG  IMPLEMENTATIONS:  Optics  and  opto-electronlc  architectures  and 
techniques  can  play  an  Important  role  in  the  study  and  implementation  of 
self-programming  networks  and  in  speeding-up  the  execution  of  learning 
algorithms.  We  have  done  some  exploratory  work  in  this  regard  to  see  how 
the  neurons  In  an  opto-electronlo  analog  of  a  neural  net  can  be  partitioned 
Into  groups  with  specific  Interconnection  patterns.  Here,  for  example,  a 
method  for  partitioning  an  opto-eleotronlc  analog  of  a  neural  net  into 
Input,  output,  and  Internal  units  with  the  selective  communication  pattern 
desoribed  earlier  to  enable,  stochastic  learning,  i.e.,  carrying  out  a 
simulated  annealing  learning  algorithm  in  the  context  of  a  Boltzmann  machine 
formalism  is  described,  (see  Fig.  1(a)).  The  arrangement  shown  In  Fig.  1(a) 
derives  from  the  neural  network  analogs  we  described  earlier  [2].  The 


Fig.  1.  Partitioning  concept  (a)  and  method  for  rapid  determination  of  the 
net's  energy  E. 

network,  consisting  of  say  N  neurons,  is  partitioned  into  three  groups.  Two 
groups,  V1  and  V2, represent  visible  or  exterior  units  that  can  be  used  as 

Input  and  output  units  respectively.  The  third  group  H  are  hidden  or 
internal  units.  The  partition  is  such  that  N1^N24,N^>N  where  subscripts 

1,2,3  on  N  refer  to  the  number  of  neurons  in  the  VirV2  and  H  groups 

respectively.  The  Interconnectivity  matrix,  designated  here  as  W^,  is 

partitioned  into  nine  submatrices,  A,B,C,D,E,  and  F  plus  three  zero  matrices 
shown  as  blackened  or  opague  regions  of  the  W^  mask.  The  LEO  array 


represents  the  state  of  the  neurons,  assumed  to  be  unipolar  binary  LED  on  - 
neuron  firing,  LED  off  -  neuron  not-firing.  The  W^j  mask  represents  the 

strengths  of  interconnection  between  neurons  in  a  manner  similar  to  earlier 
arrangements  [2].  Light  from  the  LEDs  is  smeared  vertically  over  the 

mask  with  the  aid  of  an  anamorphic  lens  system  (not  shown  in  Fig.  1(a))  and 
light  emerging  from  rows  of  the  mask  is  focused  with  the  aid  of  another 
anamorphic  lens  system  (also  not  shown)  onto  elements  of  the  photodetector 
(PD)  array.  Also  we  assume  the  same  scheme  utilized  in  [2]  for  realizing 
bipolar  values  of  in  incoherent  light  is  adopted  here,  namely  by 

separating  each  row  of  the  mask  into  two  subrows  and  assigning  positive 

values  to  one  subrow  and  negative  values  to  the  other,  then  focusing 

light  emerging  from  the  two  subrows  separately  onto  pairs  of  adjacent 
photosite  connected  in  opposition  in  the  ,  V2  and  H  segment  of  the 

photodetector  array.  Submatrix  A  with  N1  x  M1  elements,  provides  the 

interconnection  weights  of  units  or  neurons  within  group  V1 .  Submatrix  B 

with  Ng  x  N2  elements,  provides  the  interconnection  weights  of  units  within 

Vg.  Submatrices  C  (of  N1  x  elements)  and  D  (of  x  N1  elements)  provide 

the  interconnection  weights  between  units  of  V.  and  H  and  submatrices  E  (of 

N2  x  tfj  elements)  and  F  (of  x  N2)  provide  the  Interconnection  weights  of 

units  of  V2  and  H.  Units  in  V1  and  V2  can  not  communicate  with  each  other 

directly  because  locations  of  their  Interconnectivity  weights  in  the 

matrix  or  mask  are  blocked  out  (blackened  lower  left  and  top  right  portion 
of  W^j).  Similarly  units  within  H  do  not  communicate  with  each  other 

because  locations  of  their  interconnectivity  weights  in  the  mask  are 

also  blocked  out  (center  blackened  square  of  W^).  The  LED  element  6  is 

always  on  to  provide  a  fixed  or  adaptive  threshold  level  to  all  other  units 
by  contributing  to  the  light  focused  onto  only  negative  photosites  of  the 
photodetector  (PD)  arrays. 

By  using  a  computer  controlled  nonvolatile  spatial  light  modulator  to 
implement  the  mask  in  Fig.  1(a)  and  including  a  computer/controller  as 

shown  the  scheme  can  be  made  self-programming  with  ability  to  modify  the 
weights  of  synaptic  links  between  its  neurons  to  form  Internal 
representations  of  the  associations  or  patterns  presented  to  it.  This  is 
done  by  fixing  or  clamping  the  states  of  the  V,  (input)  and  V2  (output) 

groups  to  each  of  the  associations  we  want  the  net  to  learn  and  by  repeated 
application  of  the  simulated  annealing  procedure  with  Boltzmann,  or  other, 
stochastic  state  update  rule  and  collection  of  statistics  on  the  states  of 
the  neurons  at  the  end  of  each  run  when  the  net  reaches  thermodynamic 
equilibrium. 

For  each  clamping  of  the  V1  and  V2  units  to  one  of  the  associations, 

annealing  is  applied,  starting  from  an  arbitrary  Wy,  with  switching  states 

of  units  in  H  until  thermodynamic  equilibrium  is  reached.  The  state  vector 
of  the  entire  net,  which  represents  a  state  of  global  energy  minimum.  Is 
then  stored  by  the  computer.  This  procedure  is  repeated  for  each 


association  several  times  recording  the  final  State  vectors  every  time.  The 
probabilities  of  finding  the  i-th  and  j-th  neurons  in  the  same  state  are 

then  obtained.  Next  with  the  output  units  V2  unclamped  to  let  them  free  run 

like  the  H  units  the  above  procedure  is  repeated  for  the  same  number  of 

r 

annealings  as  before  and  the  probabilities  P^  are  obtained.  The  weights 

« 

W^  are  then  incremented  by  AW^  •  n< P i j i j  ^  wh®r®  n  is  a  constant  that 

controls  the  speed  and  efficacy  of  learning.  Starting  from  the  new  W^  the 

above  procedure  is  repeated  until  a  steady  W^j  is  reached  at  which  time  the 

learning  procedure  is  complete.  Learning  by  simulated  annealing  requires 
calculating  the  energy  E  of  the  net  [3], [5].  A  simplified  version  of  a 
rapid  scheme  for  obtaining  E  opto-electronically  is  shown  in  Fig.  1(b).  A 
slight  variation  of  this  scheme  that  can  deal  with  the  bipolar  nature  of 

would  actually  be  utilized.  This  is  not  detailed  here  because  of  space 
limitation. 

3.  REMARKS:  The  partitioning  architecture  described  is  extendable  to 
multilayered  nets  of  more  than  three  layers  and  to  2-D  arrangement  of 
neurons.  Learning  algorithms  in  such  layered  nets  lead  to  multivalued  W^. 

Therefore  high-speed  computer  controlled  SLMs  with  graded  pixel  response  are 
called  for.  Methods  of  reducing  the  dynamic  range  of  W^  or  for  allowing 

the' use  of  W^  with  ternary  weights  are  however  under  study  to  enable  the 

use  of  commercially  available  nonvolatile  SLM  devices  that  are  mostly 
binary  e.g.,  Litton's  MOSLM. 
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APPENDIX  VII 


Phased  Array  Antenna  Pattern  Synthesis  By 

Simulated  Annealing 

N.H.  Farhat  and  B.  Bai 


ABSTRACT:  A  new  procedure  is  described  for  optimum  phased  array  syn¬ 
thesis.  The  synthesis  is  optimised  by  a  simulated  annealing  process  in  which  the 
energy  function  is  directly  related  to  the  far  field  intensity  of  a  phased  array. 
Numerical  simulation  results  are  presented.  A  possible  optical-digital  hybrid  im¬ 
plementation  that  can  perform  the  required  computation  at  higher  speed  than 
a  pure  digital  implementation  is  discussed. 


1  INTRODUCTION 

For  the  synthesis  of  an  antenna  array  with  uniformly  spaced  elements  it  is 
well  known  that,  if  the  current  distribution  function  is  not  restricted,  a  Dolph- 
Chebysher  distribution  function  over  the  antenna  gives  rise  to  an  optimum  pat¬ 
tern  which  has  the  lowest  sidelobe  level  for  a  specified  mainbeam  width  [1]. 
However,  if  for  the  purpose  of  easy  practical  implementation,  the  distribution 
function  is  restricted  to  some  specific  set,  other  methods  have  to  be  investigated 
and  used  for  optimum  synthesis.  The  simulated  annealing  method  presented 
here  is,  by  our  study,  one  of  the  choices  for  optimum  synthesis  of  phased  ar¬ 
rays  with  restricted  distribution  functions.  The  synthesis  is  optimal  in  the  sense 
that  the  lowest  sidelobe  level  is  achieved  while  the  specified  mainbeam  width  is 
wiMfitaiimi.  This  method  can  be  used  for  both  microwave  phased  arrays  and 
optical  arrays,  hi  our  study  so  far,  ws  have  been  mainly  concerned  with  optical 
arrays,  which  appear  to  be  technologically  feasible  with  present  electronic  and 

‘The  authors  are  with  tha  Electro-Optics  aad  Microwars  Optics  Laboratory,  tbs  Moors 
School  of  Electrical  Engineering,  University  of  Pennsylvania,  Philadelphia,  PA  19104. 


optical  technologic.  Hence,  the  parameters  assumed  in  the  simulations  below 
are  relevant  to  the  optical  case,  but  the  method  and  conclusions  apply  to  phased 
arrays  in  general. 

2  SIMULATED  ANNEALING  METHOD 

Metropolis  e t  al.  introduced  the  simulated  algorithm  for  calculat- 

ing  the  properties  of  any  system  of  interacting  individual  molecules  [2].  The 
algorithm  was  previously  applied  to  combinatorial  optimisation  problems,  in¬ 
cluding  physical  design  of  computers,  and  the  traveling  salesman  problem  [3]. 
The  method  can  be  extended  for  general  optimisation  problems.  For  the  system 
to  be  optimised,  an  "energy*  or  "cost”  function  E  is  first  established  and  a 
dynamic  variable  T,  the  "temperature*  of  the  system,  is  choeen  to  control  the 
precise  Starting  at  a  high  "temperature*,  the  system  is  slowly  cooled  down, 
until  the  system  "Creeses*  and  reaches  the  optimum  state  in  a  manner  similar  to 
annealing  of  a  crystal  during  growth  to  reach  a  near  perfect  structure.  At  each 
"temperature”,  a  change  in  the  system  is  made  according  to  a  certain  rule,  and 
then  the  "energy”  or  "cost”  change  of  the  system  A E  is  calculated.  If  A  E  <  0, 
the  system  alteration  is  retained  and  the  process  is  continued.  The  acceptance 
or  rejection  of  the  alteration  or  change  of  grain  of  the  system  when  A E  >  0  is 
treated  probabilistically.  Accordingly,  the  Boltzman  factor  /(A£)  =  exp(-j^) 
is  calculated,  where  K  is  a  constant  whose  dimension  depends  on  the  dimen¬ 
sions  of  A E  and  T.  Then  a  random  number  R  uniformly  distributed  in  the 
interval  [0,1)  is  choeen.  If  R  <  /(A15),  the  change  of  grain  is  retained;  on  the 
other  hand,  if  R  >  f{AE),  the  change  is  discarded,  that  is,  the  system  before 
change  is  used  for  the  next  step  of  the  process.  This  procedure  is  repeated  for 
each  "temperature"  until  the  system  is  optimised  by  arriving  at  a  global  energy 
minimum.  The  choices  of  K  and  the  initial  T  are  crucial  for  the  success  and 
speed  of  convergence  of  simulated  annealing  process.  Because  of  the  probabilis¬ 
tic  Boltsman  selection  rule  of  the  A  E  >  0  case,  the  process  can  always  get  out  of 
a  local  minimum  of  the  "energy"  or  "cost”  function  in  which  it  could  get  trapped 
and  proceed  to  the  desired  global  minimum.  This  makes  simulated  annealing 
different  from  iterative  improvement  procedure  [3j-[6j. 


3  PHASED  ARRAY  SYNTHESIS 


For  «  2-dimensional  phased  array  with  (2A/+1)  x  (21V+1)  identical  subapertnres, 
the  far  field  is  the  product  of  its  subaperture  factor  and  its  array  factor.  The 
array  factor  can  be  written  as, 

hi  M 

Q{f*J »)*  £  £  Wr>«.«p{-j2x(»ni4/,  +  nBf,)\  (1) 

mm-Hnm-N 

where,  A  and  B  are  the  spacing  between  elements  in  an  orthogonal  array  coor¬ 
dinate;  f*  m  ft  33  xfc*  x  and  y  being  the  transverse  far  field  coordinates, 
and  Rq  is  the  distance  between  point  (x,  y)  and  the  array  origin. 

Since  the  subaperture  factor  is  fixed  for  a  given  subaperture  function  and 
it  varies  much  more  slowly  than  the  array  factor,  the  effect  of  the  subaperture 
factor  is  insignificant  and  hence  ignored  in  the  present  synthesis.  The  array 
factor  Q(/.,/f)  will  be  studied  by  changing  the  weighting  factor  W**,  for  all 
possible  m  and  n  to  achieve  the  optimum  radiation  pattern. 

The  "energy"  or  “cost"  function  in  simulated  annealing  can  be  established  in 
many  different  ways  for  phased  array  synthesis.  Since  our  primary  work  is  done 
for  optical  and  infrared  phased  arrays  and  this  kind  of  array  has  a  relatively 
large  size  compared  to  the  wavelength  A  used  (A  on  the  order  of  10~*  m),  the 
beam  width  is  very  small  (on  the  order  of  10“ 1  degrees).  It  is  important  to. 
achieve  the  lowest  aidelobe  level.  For  this  purpose,  one  obvious  way  to  choose 
the  ‘‘energy*  or  “cost"  function  for  phased  array  synthesis  would  be  the  energy 
outside  a  specified  main  lobe.  When  this  energy  function  is  minimized,  under  the 
assumption  that  the  total  energy  remains  constant,  it  could  be  expected  that 
the  energy  would  become  concentrated  in  the  main  beam.  Consequently,  the 
relative  aidelobe  level  could  be  minimized.  However,  from  simulations  we  have 
run,  it  turns  out  that  this  energy  constraint  cannot  minimize  the  sidelobe  level, 
since  minimizing  the  energy  does  not  necessarily  constrain  the  peak  sidelobe 
level.  In  the  simulation  presented  here,  the  intensity  E  associated  with  the 
highest  sidelobe  level  is  used  as  the  "energy*  or  "cost"  function  for  simulated 
annealing.  Of  course,  if  necessary,  a  constraint  about  the  beamwidth  could  be 
included  in  the  "energy"  or  "cost"  function.  When  the  "energy"  is  minimized, 
the  relative  sidelobe  level  is  now  also  minimized.  The  beamwidth  at  the  final 
run  will  be  taken  as  the  optimum  beamwidth  associated  with  the  final 


weights.  The  relationship  of  beamwidth  and  sidelobe  level  is  then  optimum  in 
the  sense  of  Dolph-Chebyshev,  that  is,  the  sidelobe  level  is  the  lowest  for  the 
given  beamwidth. 

4  SIMULATION 

Since  the  simulated  annealing  process  is  the  same  for  a  2-dimensional  array  as 
for  a  1-dimensional  array  and  considering  the  limited  computation  power  of  the 
MICRO  PDP  — 11  computer  available  for  our  simulation,  a  1-dimensional  array 
is  used  in  our  study.  The  l-dimensional  array  is  also  assumed  to  be  a  continuous 
one  with  many  desired  pixels.  Each  of  the  pixels  acts  as  a  subaperture.  The 
simulation  starts  with  a  uniform  distribution  (all  subapertures  with  1  or  -1). 
The  far  field  pattern  is  shown  in  Fig.l(a)  for  this  uniform  distribution  case  and 
has  the  usual  —  13.3(dJ?)  sidelobe  level.  The  distribution  function  is  restricted 
to  the  set  of  1,  ->1,  which  means  real  transmittance  with  binary  phase  (phase=  0 
or  it  )  modulation.  Then,  simulated  *w«*»n«g  is  carried  out  by  just  changing 
the  sign  of  each  subaperture  in  turn,  calculating  the  intensity  of  the  highest 
sidelobe  and  applying  the  algorithm  for  each  change.  At  each  “temperature”  for 
the  simulated  annealing  algorithm,  if  either  the  number  of  accepted  changes  on 
randomly  selected  elements  exceeds  £  of  the  total  subaperture  numbers  or  the 
number  of  attempts  exceeds  10  times  the  total  subaperture  numbers,  the  pro¬ 
cess  goes  to  the  next  lower  “temperature” .  If  the  desired  number  of  acceptances 
is  not  achieved  at  three  successive  “temperatures’1,  the  process  is  considered 
“frozen”  and  the  optimum  pattern  is  attained.  The  final  optimum  result  is 
given  in  Fig.l(b).  This  simulation  is  done  for  an  array  of  41  subapertures. 
The  subaperture  size  is  assumed  to  be  equal  to  the  spacing  A  between  subaper- 
tures  and  is  taken  to  be  61A  (figure  relevant  for  optical  arrays).  The  element 
distribution  function  that  gives  the  final  optimum  result  is  shown  in  table  1. 
This  function  is  seen  to  be  asymmetric.  From  formula  (1)  it  is  seen  that  two 
elements  with  opposite  index,  say  n  *  -10  and  n=  10  for  the  1-dimensional 
case,  have  different  phases  and  therefore  their  contributions  to  the  far  field  fac¬ 
tor  are  different.  In  general,  the  array  factor  given  in  formula  (1)  is  a  complex 
quantity.  Although  it  is  frequently  thought  that  the  optimum  array  factor  is 
real,  there  seems  to  be  no  reason  to  exclude  a  complex  quantity  as  an  optimum 


array  factor.  Even  if  the  array  factor  is  a  complex  quantity  and  asymmetric,  it 
can  be  shown  that  the  far  field  (power)  pattern  is  symmetric  for  arrays  with  real 
weights.  In  our  simulation,  the  weights  are  real  (1  or  —I)  and  the  pattern  in 
Flg.l(b)  is  therefore  symmetric.  The  distribution  function  in  table  1  for  the 
optimum  pattern  shown  in  Fig.l(b)  is  for  elements  —18,  —15,  and  —12  with 
weights  equal  to  -1;  if  we  make  elements  +18,  +15,  and  +12  also  with  weights 
equal  to  —1,  which  gives  a  symmetric  array  factor,  the  pattern  obtained  will  be 
that  shown  in  F!g.l(c)  and  is  seen  not  to  be  optimum.  From  the  simulation, 
it  is  found  that  the  result  of  the  simulated  annealing  procedure  is  relatively  ro¬ 
bust  for  different  initial  distribution  functions,  that  is,  it  eventually  gives  the 
same  optimum  pattern  for  different  initial  distribution  functions,  as  long  as  the 
initial  distribution  function  forms  a  pattern  with  a  single  main  beam.  The  final 
weighting  function  in  table  1  was  obtained  from  a  uniform  initial  distribution 
function  (all  weights^  1)  and  it  is  interesting  that  switching  of  weights  happens 
for  only  a  few  elements  (three).  The  function  in  table  1  gives  an  optimum 
pattern  only  for  broad  side  direction.  Like  the  Dolph-Chebyshev  optimum  func¬ 
tion,  the  optimum  function  given  by  simulated  annealing  is  different  for  different 
steering  angles,  but  an  optimal  weight  distribution  can  be  computed  for  each 
steering  angle.  From  the  simulation  result,  it  is  seen  that  the  optimum  far  field 
pattern  has  similar  features  to  the  pattern  given  by  the  Dolph-Chebyshev  dis¬ 
tribution  function.  In  the  Dolph-Chebyshev  pattern,  all  the  sidelobes  have  the 
same  level  for  a  specified  beamwidth.  A  numerical  example  in  [1]  shows  an  8- 
element  array  (element  separation  d  =  0.5A)  with  25.8(dB)  sidelobe  level  and 
40.8*  beamwidth.  The  optimum  pattern  given  by  our  simulation  shows  nearly 
equal  level  side  lobes  which  are  minimized  for  the  given  beamwidth. 

5  DISCUSSION 

Simulated  annealing  is  a  modification  of  the  iterative  improvement  algorithm  [4]. 
It  is  physically  more  meaningful  and  can  be  computed  more  systematically  than 
the  iterative  improvement  (4j.  Physically,  the  simulated  annealing  process  is 
analogous  to  the  cooling  of  atoms  in  crystal  growth:  careful  annealing  produces 
a  defect-free  crystal,  rapid  produces  a  defective  crystal  or  glass  [3]. 

The  probabilistic  treatment  with  the  probability  function  P(A£)  =  exp(—j£) 


provides  a  way  to  accept  the  unfavorable  changes  and  is  easy  to  compute.  From 
our  simulation,  it  has  been  found  that  the  simulated  algorithm  seems 

always  to  give  better  performance  than  the  iterative  improvement  algorithm. 

Since  simulated  annealing  is  a  modified  iterative  improvement  process  ,  it 
talus  a  relatively  long  time  to  do  an  optimisation  problem  just  as  iterative 
improvement  does  in  a  computer  calculation.  The  phased  array  synthesis  in 
our  simulation  runs  for  one  hour  or  so  for  an  array  of  41  elements  on  a  MICRO 
PDP-11  computer.  Finding  an  efficient  scheme  to  reduce  the  excessive  amounts 
of  computer  time  on  most  optimum  problems  has  always  been  of  concern  [5}-{7]. 
Otherwise,  if  enough  computation  power  is  available,  iterative  improvement  can 
be  run  from  random  starts  for  many  times  to  approach  the  optimum  state.  Fast 
opto-digital  computing  schemes  similar  to  those  described  in  [8]  may  also  be 
considered  for  phased  array  synthesis  by  simulated  it  is  understood 

that  the  for  field  is  the  Fourier  transform  of  the  array  distribution  function.  An 
optical  lens  can  be  used  for  computing  the  Fourier  transform  as  the  distribution 
function  is  inputed  to  the  front  focal  plane  of  the  lens  via,  for  example,  an  appro* 
priate  computer  driven  spatial  light  modulator  (SLM).  The  Fourier  transform 
in  the  back  focal  plane  can  be  recorded  and  fed  to  the  computer /controller  to 
make  the  simulated  annealing  decision.  The  outcome  is  fedback  to  the  SLM  to 
change  the  distribution  function  in  the  front  focal  plane.  The  hybrid  opto-digital 
scheme  will  do  the  Fourier  transform  instantly.  Li  this  fashion  the  computation 
associated  with  the  Fourier  transform  can  be  virtually  assuming  a 

high  speed  SLM  and  computer  interface  are  utilised.  An  optoelectronic  Bolts- 
man  machine  for  accelerating  the  selection  rule  has  also  been  proposed  earlier 
in  [8].  This  process  can  be  repeated  for  every  step  in  simulated  Mingling  Also, 
a  Cauchy  probability  selection  rule,  instead  of  the  Boltsman  selection  rule,  can 
be  used  to  speed  up  the  whole  process  further  (7]. 
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Fig.  1.  Simulated  annealing  result,  (a)  Far  field  pattern  for  the 
uniform  distribution;  (b)  Pattern  obtained  by  simulated 
annealing;  (c)  Pattern  obtained  from  a  symmetrized  version 
of  the  weights  distribution  obtained  by  simulated  annealing. 
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